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SECTION  1 


INTRODUCTION 

A  major  concern  with  the  increasing  use  of  composite  materials  and 
low  voltage  electronics  is  the  amount  of  electromagnetic  (EM)  coupling  to  the 
interior  of  an  aircraft  and  to  the  cables  and  electronic  devices  within  it. 

The  introduction  of  boron/epoxy,  graphite/ epoxy ,  and  Kevlar/epoxy  composite 
materials  as  structural  elements  in  modern  airframes  will  result  in  a  substan¬ 
tial  reduction  in  airframe  weight,  due  to  the  high  strength-to-weight  ratios 
of  these  materials.  The  use  of  these  new  composite  materials  has  raised  ques¬ 
tions  relative  to  the  aircraft  vulnerability  resulting  from  the  effects  of 
lightning,  high  power  radar,  nuclear  electromagnetic  pulse  (EMP) ,  and  precipi¬ 
tation  static.  The  problems  are  further  compounded  by  the  fact  that  these 
materials  are  relatively  easy  to  construct,  and  have  resulted  in  a  prolifera¬ 
tion  of  available  composite  materials. 

This  final  report  on  Office  of  Naval  Research  Contract  N00014-78-C- 
0673  describes  methods  for  determining  the  shielding  provided  by  an  aircraft's 
exterior  surface  and  the  coupling  of  the  interior  fields  to  cables  and  trans¬ 
mission  lines  within  aircraft  cavities.  This  data  is  used  to  determine  whether 
devices  commonly  found  on  aircraft  will  be  subject  to  upset  or  burnout. 

The  penetration  of  an  external  electromagnetic  field  into  the  in¬ 
terior  of  a  homogeneous  shell  enclosure  has  been  widely  studied  and  various 
formulations  can  be  found  in  the  open  literature.  Analytical  solutions  are 
available  for  the  canonical  geometries  of  twin  parallel  plates,  a  spherical 
shell,  and  an  infinitely  long  circular  shell.  The  utility  of  these  solutions 
is  manifested  in  a  transfer  function  relating  the  interior  field  at  a  point  to 
the  field  that  would  exist  there  in  the  absence  of  the  shell.  This  result  is 
usually  presented  in  the  frequency  domain  and,  for  low  frequencies,  obviates 
a  relatively  simple  relationship  between  the  interior  field  and  the  excitation 
field.  For  canonical  shell  geometries,  this  low  frequency  transfer  function 
is  written  in  terms  of  shell  wall  conductivity  and  thickness  and  shell  enclos¬ 
ure  volume— to-sur face  ratio.  Application  to  noncanonical  geometries  can  be 


made  as  long  as  volume- Co-surface  ratios  are  known.  As  the  frequency  content 
of  the  excitation  spectrum  becomes  large  enough  so  that  the  electrical  size  of 
the  shell  cross  section  becomes  resonant  (on  the  order  of  a  free-space  wave¬ 
length),  then  the  low  frequency  transfer  function  is  no  longer  adequate  to 
describe  the  penetrability  of  the  shell.  In  this  case  it  is  necessary  to 
resort  to  approximate  numerical  techniques. 

In  this  report,  several  models  of  the  shell  coupling  mechanism  are 
analyzed  with  frequency  regions  of  validity  from  dc  to  several  gigahertz 
depending  on  shell  cross  section  dimensions.  Two-dimensional  shell  enclos¬ 
ures  of  infinite  extent  in  one  dimension  are  considered  in  order  to  facilitate 
the  computer  program  solutions.  Results  from  these  theoretical  enclosures  are 
applicable  to  physically  realizable  three-dimensional  enclosures  which  are 
long  compared  to  their  cross  section  dimension  (i.e.,  some  airplane  fuselage 
and  wing  sections). 

In  Section  2,  various  integral  equation  formulations  are  outlined 
for  determining  the  induced  current  density  on  perfectly  conducting  two- 
dimensional  cylindrical  shells  having  an  arbitrary  cross  section.  Though 
no  penetration  occurs  if  the  shell  wall  is  a  perfect  conductor  (a  =  00 ) ,  the 
current  density  on  the  exterior  surface  caused  by  an  incident  field  is  much 
the  same  as  that  on  a  highly  conducting  (but  finite  a)  shell.  The  various 
integral  equations  are  solved  by  the  method  of  moments,  and  specific  matrix 
operators  are  defined  for  later  use.  A  user-oriented  computer  program  is 
given  in  Appendix  A  with  sample  input/output  data. 

In  Section  3,  two  approximate  shell  coupling  formulations  are  pre¬ 
sented  for  two-dimensional  shells  having  an  arbitrary  cross  section.  The 
matrix  operators  defined  in  Section  2  are  used  in  the  moment  method  solution 
of  the  resulting  integral  equations.  Comparison  with  the  exact  series  solu¬ 
tion  for  the  circular  cross  section  is  used  as  a  check.  A  user-oriented  com¬ 
puter  program  is  given  in  Appendix  B  with  sample  input/output  data. 


The  exact  series  solution  for  a  normally  incident  plane  wave  exciting 
a  shell  of  circular  cross  section  is  summarized  in  Section  4.  Simple  low- 
frequency  formulas  are  derived  for  the  interior  fields.  The  case  of  axial 
electric  current  line  source  excitation  is  also  analyzed  for  later  application 
to  the  case  of  near-strike  lightning. 

The  preceding  analysis  is  presented  only  in  the  frequency  domain. 
However,  if  the  excitation  spectrum  is  sufficiently  band  limited  and  the 
transfer  functions  for  the  canonical  geometries  are  valid  over  that  frequency 
range,  then  analytical  expressions  for  the  interior  field  may  be  derived  in 
the  time  domain.  This  is  done  in  Section  5  for  near-strike  lightning,  direct- 
strike  lightning,  and  a  nuclear  electromagnetic  pulse. 

These  techniques  may  be  integrated  to  provide  an  accurate  descrip¬ 
tion  of  the  penetration  fields  inside  a  homogeneous  two-dimensional  enclosure. 
The  effect  of  this  interior  field  on  circuits  situated  inside  the  enclosure  is 
of  primary  importance  and  an  interactive  computer  program  was  written  for  this 
purpose  which  utilizes  the  results  of  Section  5.  This  program  is  described  in 
Appendix  C. 
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SECTION  2 


PERFECTLY  CONDUCTING  CYLINDERS  OF  ARBITRARY 
CROSS  SECTIONAL  SHAPE 

Although  no  field  penetrates  an  enclosure  with  perfectly  conducting 
walls,  the  electric  current  density  induced  on  the  wall  exterior  due  to  an 
external  field  will  not  differ  greatly  from  that  in  the  case  of  walls  having 
a  finite  conductivity  of  5  :  1000.  In  fact  some  shell  coupling  formulations 
require  a  "short-circuit  current"  which  is  used  to  excite  an  equivalent  prob¬ 
lem  for  the  shell.  This  is  simply  the  current  flowing  on  the  outside  surface 
of  the  shell  when  the  walls  are  perfectly  conducting  (a  =  . 

The  purpose  of  this  section  is  to  present  the  E-field,  H-field,  and 
combired-f ield  integral  equation  formulations  for  perfectly  conducting  cylinders 
of  infinite  length  and  arbitrary  cross  section  illuminated  by  a  normally  inci¬ 
dent  plane  wave.  The  E-field  equation  is  obtained  by  requiring  that  the  total 
tangential  electric  field  be  zero  on  the  contour  C  defining  the  cylinder  cross 
section.  The  H-field  equation  is  obtained  by  requiring  that  the  total  tan¬ 
gential  component  of  magnetic  field  equal  zero  just  inside  C.  The  combined- 
field  equation  is  obtained  by  taking  a  linear  combination  of  the  E-field  and 
H-field  equations.  These  integral  equations  are  written  in  matrix  form  by 
using  a  method  of  moments  Galerkin  procedure.  The  unknown  electric  current  on 
C  is  then  solved  for  by  standard  matrix  methods.  The  exact  series  solution  is 
also  presented  for  comparison  purposes  when  C  is  a  circle.  Generalization  to 
oblique  incidence  is  also  outlined  but  not  programmed.  Computer  programs  are 
documented  in  Appendix  A  for  the  E-field,  combined-field,  and  exact  series 
solutions . 

2 . 1  INTRODUCTION 

The  E-field  and  H-field  formulations  for  this  problem  are  well  known 
[1,2]  and  some  E-field  computer  programs  have  been  documented  [3,4].  The  cross 
section  of  the  cylinder  is  defined  by  the  contour  C,  which  will  be  approximated 
by  straight  line  segments.  For  each  formulation,  an  integral  equation  is  written 
involving  an  equivalent  electric  current  which  replaces  the  conducting  contour  C. 


The  integral  equation  is  then  solved  for  the  electric  current  by  a  method  of 
moments  Galerkin  procedure  [2,5].  Once  this  electric  current  is  determined, 
quantities  such  as  the  scattered  far  field  pattern  and  radar  cross  section 
may  be  easily  computed. 

For  a  normally  incident  plane  wave,  as  discussed  in  the  next  sub¬ 
section,  the  total  field  may  be  expressed  as  the  superposition  of  a  TE  (trans¬ 
verse  electric  to  z)  part  and  a  TM  (transverse  magnetic  to  z)  part.  Since  it 
is  not  the  purpose  of  this  section  to  rigorously  derive  the  different  formula¬ 
tions,  they  are  presented  with  brevity  in  Subsections  2.2  through  2.4  where 
explicit  formulas  are  given  as  an  aid  in  understanding  the  programs.  Formulas 
for  the  scattered  field  pattern  are  given  in  Subsection  2.5.  For  comparison 
purposes,  one  may  check  the  programs  against  the  exact  series  solution  pre¬ 
sented  in  Subsection  2.6.  A  generalization  to  oblique  incidence  is  given  in 
Subsection  2.7  for  the  E-field  integral  equation.  The  special  problem  of  a 
longitudinal  impressed  current  excitation  is  considered  in  Subsection  2.8. 
Finally,  detailed  instructions  for  using  the  computer  programs  are  included 
in  Appendix  A. 

2.1.1  Excitation 

The  cylinder  is  assumed  infinite  in  the  z  direction  and  is  defined 
by  the  two-dimensional  contour  C  lying  in  the  x-y  plane.  The  shape  of  C  is 
independent  of  z.  For  simplicity,  the  cylinder  is  illuminated  by  a  normally 
incident  (kz  =  0)  uniform  plane  wave.  A  time  dependence  of  e~*u)t  is  implicit 
throughout.  This  excitation  gives  rise  to  a  scattered  field  which  is  also 
independent  of  z.  Thus  the  TE  case  (magnetic  field  parallel  to  z)  and  the  TM 
case  (electric  field  parallel  to  z)  may  be  treated  separately.  The  source  of 
the  scattered  field  is  postulated  to  be  an  electric  current  which  takes  the 
place  of  the  perfect  conductor  and  which  is  defined  on  T.  It  is  separated 
into  a  z  component  (TM  case)  and  a  transverse  component  (TE  case)  directed 
along  C.  For  the  more  general  excitation,  where  0  #  tt/2,  the  two  components 
of  electric  current  are  coupled  and  thus  both  polarizations  must  be  treated 
together  as  indicated  in  Subsection  2.7. 


la  terms  of  its  TE  and  TM  parts,  the  incident  field  may  be  expressed 


as 


E1  =  Eie  +  Eih 


(1) 


and 


H1  =  H16  +  Hlh 


(2) 


The  superscripts  e  and  h  denote  TM  and  TE,  respectively.  These  parts  are 
written  explicitly  as 
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(6) 


in  terms  of  the  coordinate  system  of  Fig.  2-1. 

In  the  above,  k  and  n  are  the  wave  number  and  impedance,  respectively, 

A 

of  the  space  surrounding  the  cylinder.  The  unit  vector  k  is  defined  by 


-k  =  x  cos 


,i  ~  ,  i 

$  +  v  sm  $ 


(7) 


where  X  is  the  angle  of  incidence  measured  counterclockwise  from  the  x  axis. 
The  vector  r_  is  from  the  origin  to  a  point  on  the  x-v  plane.  The  real  numbers 
a  and  b  are  chosen  so  that  aL  +  b  =1  and  choices  of  aeJ  and  beJ  determine 
the  polarization  of  the  incident  field,  which  is  elliptical  in  general.  For 
example,  a  choice  of  ae-'a  =  1  and  be-'®  =  0  gives  the  linearly  polarized  TM 
case.  A  choice  of  ae"'a  =  1//2  and  be^  =  j//2  gives  a  left-hand  circular 
polarization.  For  simplicity,  ae^a  and  be^  are  taken  to  be  equal  to  unity 
here. 


atal!  at — 


. 
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Original  Problem:  A; 
Wave  Normally  Incidei 


The  total  field  inside  C  is  zero  and  the  total  field  outside  C  is 
written  as  (E^  +  E^  ,  +  Ji  )  where  (E^,  H^)  is  the  field  which  exists  every¬ 

where  without  the  cylinder  present.  The  scattered  field  (ES ,  HS)  is  written 
in  terms  of  an  electric  current  J  as  [6]: 


ES  (J  )  = 


-j  n 


kA  +  £  V  7  •  A 


(8) 


HS  =  HS  (J  )  =  V  x  A 


(9) 


where  the  magnetic  vector  potential  A  is  given  by 


A 


1_ 


l 


77  /  J  (f)  H 


(2) 


(k  r  -  r’(t')l)  dt’ 


(10) 


The  symbols  £  and  H  denote  electric  and  magnetic  field  operators, 

respectively.  The  domain  of  integration  in  Equation  (10)  is  restricted  to  C, 

where  is  defined  in  terms  of  t',  the  arc  length  variable  along  C.  The 

vectors  jr  and  x'  denote  field  and  source  points,  respectively,  in  the  x-v 

,  (2) 

plane  and  denotes  the  Hankel  function  of  the  second  kind,  order  zero. 

2.1.2  Specification  of  Contour  C 

To  proceed  with  a  numerical  solution,  the  contour  C  is  approximated 
by  a  finite  number  (NO)  of  straight  line  segments  as  shown  in  Fig.  2-2.  This 
is  done  by  specifying  the  x-y  coordinates  of  the  end  points  of  each  segment 
starting  with  (x^,y^)  and  proceeding  clockwise  to  (*NC+1»  yfjC+l^‘  In  the  E-fiel-d 
formulation,  it  is  not  necessary  for  the  contour  to  be  closed.  A  closed  contour 
is  one  for  which  (x^,  y 

ever  for  the  H-field  formulation  and  hence  for  the  combined-field  formulation. 
Each  straight  line  segment  AC^  has  length  A^,  a  normal  unit  vector  n^,  and  a 
tangent  unit  vector  t^  for  integers  n  =  1,2,...,  NC.  These  unit  vectors  are 
related  by 

A  A  /V 

t  x  n  =  z  (11) 

— n  — n  — 


)  = 


(XNC+1’ 


yNC+l) ‘ 


This  requirement  must  be  met,  how- 
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Figure  2-2.  The  Contour  C  Approximated  by  NC 
Straight  Line  Segments 


2- 


The  parameter  t  is  introduced  to  represent  the  arc  length  along  C 
measured  from  the  point  (x^,  v^)  to  any  point  on  C.  Subscripted  values  of  t 
are  given  by  the  formula 


n-1 

cn  "  X  A  (12) 

i-1 

for  n  »  2,  3,  . . . ,  NC  and  with  t^  =  0. 

2.1.3  Definition  of  Expansion  Functions  and  Symmetric 
Product 


As  mentioned  earlier,  may  be  separated  into  z-directed  and 
transverse-directed  components.  This  is  written  as 


4  *  4  +  4  (13) 

Since  there  is  a  charge  associated  with  ,  it  is  desirable  that  its  represen¬ 
tation  in  terms  of  a  set  of  expansion  functions  be  differentiable.  There  is 
no  charge  associated  with  J^,  however,  but  Jz  does  become  unbounded  near  sharp 
edges  of  perfect  conductors.  With  this  in  mind,  we  define  a  set  of  triangle 
functions  as 


r  t  -  t  . 
f  m-i  ' 


-  -  w  L  i  ^  ^  L. 

t  -  t  ,  — tn-1  m-1  —  —  m 

m  m-1 


T 

— m 


(t> 


nH-l 


m+1 


t 

— m 


for  t  <  t  < 
m  — 


tm+l 


(1A) 


0 


for  t  elsewhere 
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and  a  sec  of  pulse  functions  as 


P  (t) 

— m 


1  z  for  t  <  c  <  t  , 
—  m  —  —  nri-l 


for  t  elsewhere 


(15) 


for  integers  m  *  1,  2,  . . . ,  NC  and  with  t  =  -A.„  and  t  5  r  . 

o  NC  —o  -NC 


The  electric  current  J  is  then  expanded  as 


J  (t) 
— c 


NC 

2  +  1 !  p_(t) 

n*l 


a  —a 


n  a 


(16) 


h  e 

where  I  and  I  are  complex  coefficients  to  be  determined  for  the  TE  and  TM 
n  n  r 

cases,  respectively.  For  the  TM  case,  J_c  -  0  and  for  the  TE  case,  =  0. 


In  the  Galerkin  procedure,  the  testing  functions  are  chosen  to  be 
identical  to  the  expansion  functions.  Hence,  to  carry  out  this  procedure,  a 
symmetric  product  is  defined  by 


<A,B>  = 


(17) 


with  A  and  B  defined  on  C. 


2.2  E-FIELD  FORMULATION 

The  E-field  integral  equation  is  obtained  by  setting  the  tangential 
component  of  the  total  electric  field  equal  to  zero  on  C.  This  is  written  as 


1/  ^  c  k.  i 

-  -  Eb  (J  )  -  -  Et  on  C 

n  — t  — c  n  -t 


(18) 
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where  the  extra  factor  of  k/rj  was  multiplied  through  for  later  convenience. 
The  operator  E  is  defined  by  Equations  (8)  and  (10)  and  the  subscript  t 
denotes  tangential  component  found  by  the  usual  — n  x  n  x  operation.  After 
expanding  in  terms  of  Equations  (14)  or  (15),  depending  on  the  polariza¬ 
tion  considered,  and  testing  Equation  (18)  with  the  same  functions  used  for 
expansion,  one  obtains  the  following  sets  of  matrix  equations: 


for  the  TM  case,  and 


[ze]  Te  =  vie 


[zh]  ih  -  vih 


for  the  TE  case.  The  vectors  I  and  I*5  contain  the  coefficients  of  expansion 
in  Equation  (16). 


2.2.1  Formulas  for  [Z] 

The  elements  of  the  matrices  [ Z& ]  and  [Zh)  are  given  by  the  follow¬ 
ing  formulas,  where  1  £  m  £  NC  and  1  _<  n  £  NC.  For  the  TM  case,  we  have 


Z  =  -  -  <  P  ,  E  (P  )  > 

mn  n  >  — m  — t  v-^n'  x 


-tf 


P  (t)  •  n  x  n 
~m  — m  — ni 


P  (f)  H  (k|r(t) 
— n  o  — 


-  r'(t')l)  dt'  dt 
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After  transforming  both  source  and  field  intervals  to  the  interval 
[-1,1]  and  letting  y  -  kA^,  we  obtain 


'y4r  /  /  "„<2) 

<1  *'-1 


y  ~  y 

^  u  t  -  ^  u'  t  +  R' 

2  — m  2  — n  — ra,n 


du'  du 

if  m  4  n 


mn 


(22) 


[r  /  [a  (r  (1  +  u))  +  a(r2  (1  -  4. 


du  if  m  *  n 


where  R'  is  k  times  the  vector  from  the  midpoint  of  AC  to  the  midpoint 
— m,n  r  n  r 

of  AC  .  The  function  ot  is  defined  bv 
m 


a(z ) 


z 

I 


H  (2)  (u)  du 
o 


(23) 


which  is  computed  using  Struve  functions  [7].  The  integrals  in  Equation  (22) 
are  readily  approximated  by  a  Gaussian  quadrature  integration  rule  [8,9].  For 
the  TE  case  we  have 


Zh  -  -  -  <  T  ,  ES  (x  )> 

mn  n  N  -m  —  t  — n 


-  -k.‘ 


L 

/ 


m+1 


(t)  •  n  x  n  *  [  1  +  ■=—  7  V  • 


(24) 


m-1 


/ 


n+1 


T  ( t ’ )  H  (2)  (klr(t)  -  r ' ( t ' ) j )  dt'  dt 
— n  o  —  — 


n-1 
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It  is 


where  the  unit  vector  n  resides  on  the  field  interval  AC  ,  U  AC  . 

—  m-i  m 

convenient  to  break  the  above  integral  up  into  four  parts.  Considering  the 
contribution  from  each  part  separately.  Equation  ( 24 )  is  rewritten  as 


SZh  (m  -  1,  n  -  1,  1,  1)  +  SZh  (m  -  1,  n,  1,  -1) 


+  SZh  (m,  n  -  1,  -1,  1)  +  SZh  (m,  n,  -1,  -1) 


where  the  function  SZ  is  defined  by 
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(2) 

where  H1  denotes  the  Hankel  function  of  second  kind,  order  one.  Note  that 
e  ^  h 

[Z  ]  and  [Z  ]  are  both  symmetric  matrices  so  that  one  need  only  compute  the 
upper  right  triangle  portion  of  each. 

2.2.2  Formulas  for  V1 

-»ie  -+ih 

The  elements  of  the  excitation  vectors  V  and  V  are  given  by  the 
following  formulas  where  1  <  m  <  NC.  For  the  TM  case  we  have 


/le  =  -  <  P  ,  Eie  > 

m  H  N  — m  — t 


(27) 


Y_  e 
m 


jk  *  ^  sin  -f  (k  •  tm) 


(k  •  t  ) 
2  TO 


f- 


where  R'  is  k  times  the  vector  from  the  origin  to  the  midpoint  of  AC  and  k  is 
— m  m  — 

defined  by  Equation  (7).  For  the  TE  case  we  have 


,ih 


k  ,  T  rih  \ 
n  v  ^m’  -t  > 


'm-1  m-1 


jk  •  R'  . 
-m-1 


j  b 


m-1 


jb 


m-1 


sin  b 


m-1 


m-1 


(28) 


Ytn  am  Jii  ’  & 
2  j  b  6 

-1  m 


sin  b  -1b 

m  J  m 
- ; -  -  e 


where  a  =  -n  •  k  ana  b  =  (y  /2)  k  *  t  - 
m  — m  —  mm—  — m 
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2.3  H- FI ELD  FORMULATION 


The  H-field  integral  equation  is  obtained  by  setting  the  tangential 
component  of  the  total  magnetic  field  equal  to  zero  just  inside  C.  This  is 
written  as 

-k  n  x  (J^)  =  k  n  x  H ^  on  C  (29) 

where  C  denotes  a  contour  just  on  the  -n  side  of  C.  The  factor  of  k  has  been 

multiplied  through  for  later  convenience.  The  magnetic  field  operator,  H  , 

is  defined  by  Equations  (9)  and  (10).  After  expanding  in  terms  of  Equa¬ 
tions  (13)  or  (14),  depending  on  the  polarization  considered,  and  testing 
Equation  (29)  with  the  same  functions  used  for  expansion,  one  obtains  the 
following  sets  of  matrix  equations: 

[Te]  I®  =  Iie  (30) 

for  the  TM  case,  and 

[Th]  Ih  =  Tlh  (31) 

-*h  _ 

for  the  TE  case.  The  vectors  I  and  I  again  contain  the  coefficients  of 

expansion  in  Equation  (16). 
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.3.1  Formulas  for  [T] 


0  V% 

The  elements  of  the  matrices  [T  ]  and  [T  ]  are  given  by  the  follow¬ 
ing  formulas,  where  1  <  m  <  NC  and  1  <  n  <  NC.  For  the  TM  case  we  have 


-k  <  P  ,  n  x  H  (P  ) > 
— m  —  —  — n 


=  k 


r  m+1 

/  * 


P  (t)  •  P  (t)  dt  -  yr 
“in  — n  4  j 


L 

/ 


m+1 


P  (t) 

— m 


(32) 


/ 


n+1 


(r  -  r') 


nm  X  4(C'}  x  -  r  1  I  H!(2>  (k|r  -  r  ’  j )  dt'  dt 


The  first  term  is  simply  the  Ampere's  law  contribution  to  the  inte¬ 
gral  when  the  field  point  is  on  C  .  Again,  after  some  algebra,  one  mav 
obtain 


r . 


if  m  = 


mn 


(33) 


1  1 


Y  Y  r  r  n  •  R 

m  n  /  /  — m  — m,n 

16j  K  K  i** 

-1  -1  — m,n 


(2) 


( ’ R '  ! )  du '  du 

— m,n 


if  m  ^  n 


In  the  above,  R'  is  given  bv 
— m,  n 


R' 

— ta,n 


R' 

— m,  n 


u  t 
— m 


t 

— n 


(34) 


t 
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where 


n  •  R' 

>(m,n)  =  =^—3^  H  2  ( j  R'  |) 


i  R '  , 

-m,n' 


— m,  n 


(38) 


2.3.2  Formulas  for  I 


■tx 


The  elements  of  the  excitation  vectors,  Iih  and  Iie,  are  given  by 
the  following  formulas  where  1  <  m  <  NC.  In  the  TM  case  we  have 


I16  =  k  P  ,  n  x  H1' 

m  x  — m  —  —  • 


(n 


A  V  /\ 

-  jk  •  R’  sin  (k  ♦  t  ) 

10  Y  e  - =SL 

—  m  y 


m 


(k  *  U 

— m  — m 


(39) 


For  the  TE  case  we  have 


r"h  =  k  <  T  ,  n  x  H1  > 
m  — ra  —  — 


m-l  e 


jk  •  R’  i 
—  —m-l 


jb 


m-l 


jb  .  sin  b 


m-l 


m-l 


m-l 


(40) 


y  '  & 

m  e 


jb. 


m 


sin  b 


-  e 


_jb 


2.4  COMBINED  FIELD  FORMULATION 

It  can  be  shown  [10]  that  the  E-field  or  H-field  equations  are  not 
sufficient  by  themselves  to  uniquely  determine  the  electric  current  distribu¬ 
tion,  J  .  That  is,  they  each  may  have  non-trivial  homogeneous  solutions  at 
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frequencies  which  correspond  to  internal  eigenfrequencies  of  the  closed  con¬ 
tour  C.  To  illustrate  this  correspondence,  consider  the  TE  interior  problem, 
where  the  total  electric  field  ji  inside  C  must  satisfy 

(V*  +  k^)  =  0  inside  C  (41) 

subject  to  the  boundary  condition 

E  =  0  on  C  (42) 

This  is  mathematically  identical  to  the  external  problems: 

0  H-field  formulation,  TM  case 

°  E-field  formulation,  TE  case 

Thus,  these  problems  have  the  same  eigenfrequencies.  Similarly,  for  the  TM 
interior  problem,  the  total  electric  field,  E^,  inside  C  satisfies 

(V*  +  k2)  E  =0  inside  C  (43) 

t  z 

subject  to  the  boundary  condition 

E  =  0  on  C  (44) 

z 

This  problem  is  mathematically  identical  to  the  external  problems 
°  E-field  formulation,  TM  case 

0  H-field  formulation,  TE  case 

Thus,  at  or  near  these  internal  resonant  frequencies,  the  E-field  and  H-field 
matrices  become  ill-behaved.  To  remedy  this  situation,  a  linear  combination 
of  the  E-  and  H-field  equations  is  formed: 

-k  n  x  HS  (J  )  -  6  -  (J  )  =  k  n  x  H1  +  S  £  (45) 

—  —  — c  n  -  t  — c  —  —  t 
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This  equation  is  referred  to  as  the  combined  field  formulation  and  it  can  be 
shown  [10]  that  Equation  (45)  has  a  unique  solution  for  for  any  contour  C 
at  all  frequencies  as  long  as  3  is  a  positive  real  number.  In  matrix  form. 
Equation  (45)  is  written  as 

[Te  +  3  Ze]  Ie  =  I1'+  3  V16  (46) 

for  the  TM  case,  and 

[Th  +  3  Zh3  Ih  =  I:h  +  s  Vlh  (47) 

for  the  TE  case.  The  formulas  for  these  matrix  elements  are  given  in  Sub¬ 
sections  2.2  and  2.3. 


2.5  FORMULAS  FOR  SCATTERING  CROSS  SECTION 

s  s 

Once  the  electric  current  is  found,  the  scattered  field  (E.  ,  H  ) 

is  readily  computed  from  Equations  (8)  through  (10).  In  the  far-field, 

[r(  >>  X,  there  are  two  quantities  of  interest  which  are  computed  from  J^. 

One  is  the  normalized  scattered  field  pattern.  This  is  simply  a  plot  of 
i  S  S  i  .  s  S  i 

E  /E  versus  <b  for  the  TM  case  and  H  /H  versus  i>  for  the  TE  case. 

1  z  z  max1  1  z  z  max1 

The  denominator  is  the  maximum  value  of  scattered  field.  The  second  quantity 
of  interest  is  the  scattering  cross  section.  For  the  TM  case,  this  is  defined 
by  the  equation  [2] 

S  2 
E*  (r , <t>) 

a(cj>)  =  11m  2ttt  — — -  (48) 

r-«°  E  1 

z 

Using  Equation  (3)  with  ae^a  *  1  and  specializing  H  (k.|r_  -  r 1  | )  to 
large  r_,  we  obtain 


(49) 
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Note  that  V™*1  and  V™6  are  identical  to  the  E-field  excitation  vectors  V**1  and 

+ie  i 

V  ,  resoectively ,  for  measurement  angle  <£  =  ()>. 


2.6  EXACT  SERIES  SOLUTION 

Here  the  cylinder  is  defined  by  a  contour  C  which  is  a  circle  of 
radius  a.  The  angle  of  incidence  will  be  chosen  at  ■  0.  Thus  for  the 
TM  case,  the  incident  field  at  a  point  r  is  expanded  as  [6] 


E 


i 

z 


n  e 


jkx 


ikr  cos  <J> 

*  h  e 


oo 

=  n  2  e  jn  J  (kr)  cos  n<5 
n  n  n 
n=0 


(56) 


Neumann's  number  e  is  defined  by 
n 


e 

n 


1  for  n  =  0 

2  for  n  >  0 


(57) 


For  the  TE  case,  the  incident  field  on  C  is  expanded  similarly  as 


jkx 

e 


OO 

=  2]  e  j  J  (kr)  cos  nO 

rv  n  n 

n=0 


(58) 


Writing  J  in  Equation  (13)  in  terms  of  its  components  on  C  we  have 


J  =  J  z  +  J,  o 
— c  z  —  ■:(>-*- 


(59) 
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where  J  and  J. 

z  <p 


may  be  expanded  in  a  Fourier  series  on  C  of  the  form 


2 

irka 


a  +  2  X  a  cos  nd> 
o  n  r 

n=l 


(60) 


and 


TTka 


c  +  2  V  c  cos  nd> 
o  n 

n-1 


(61) 


The  coefficients  are  obtained  by  enforcing  the  boundary  conditions 

i  S  iS 

on  the  tangential  components  of  E  +  E  and  H  +  H  at  r  =  a.  They  are  given 
by  [6,11] 


Hq(2)  <ka) 


H  (ka) 

o 


(62) 


,  .,n 

c-j). 


H  (2)  (ka) 
n 


H  <2>'  (ka) 
n 


A  A 


Note  that,  in  the  above  expansion  for  J,  ,  ^  =  -t_. 


Formulas  for  the  normalized  scattered  field  pattern  are  giver  by 


fa/\  =  \/2h 


b  +  2  V  b  cos  ncf> 
o  n 

n=l 


(63) 


for  the  TM  case  and 


/o7I  =  \fTR 


d  +  2  y  d  cos  n<!> 
o  n 

n*l 


(64) 
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for  Che  TE  case.  The  coefficients  are  given  by 


b 

o 


J 


H 

o 


(ka) 


(2) 


(ka) 


d 

o 


(ka) 

(ka) 


(65) 


b 

n 


(-1) 


J  (ka) 

H  (2)  (ka) 
n 


(-D 


■T  (ka) 

H  (2)’  (ka) 
n 


2.7  OBLIQUE  INCIDENCE 

A  plane  wave  which  is  incident  at  an  angle  <p^  from  the  x  axis  and 
91  from  the  z  axis  of  Fig.  2-1  is  written  in  the  form 

~  a 

jk(x  +  k  z)  •  r 
■  — t  z  — 

E  =  E  e  (66) 

—  — o 

where  the  argument  of  the  exponential  is  defined  by 


k  ■  2tt/ A  =  uj  vnr 

k  =  sin  9*  cos  ^  x  +  sin  9*  sin  v 

k  =  cos  9^ 

z 


(67) 


r_  =  x  x  +  y  2. 

The  z  dependence  of  the  incident  field  gives  rise  to  a  z-dependent  scattered 
field  and  it  is  no  longer  possible  to  decouple  the  z-directed  and  transverse- 
directed  components  of  J^.  The  E-field  integral  equation  is  now  written  as 

E^(r,z)  =»  -E*(r,z)  r  on  C  for  all  z  (68) 
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which  must  be  satisfied  everywhere  on  C.  We  define  the  Fourier  transform 
pairs 


E(r,z)  = 


\fli T 


00 

=  f  £(r. 

2tt  J 


jk,z 

k  )  e  dk 
z  z 


(69) 


and 


-jk,z 


€(r,  k  )  =  —4=  j  E(r,z)  e  dz  (70) 

2  V2tt  J 


Taking  the  transform  of  Equation  (68)  gives 


— t  ^r*  kz^  =  ~—t  ^r’  kz^  °n  C 

The  scattered  field  IS  is  caused  by  an  equivalent  electric  current  .J(r',  z’) 
and  is  written  as 


E  (r,z)  =  - 


jn  I  k  /  J  J_  G  dt'dz’  +  ^V  V  ■  f  I  J_ 

L  *x  K'  z' 


G  dt'dz' 


(72) 


where  G  is  the  Green's  function  defined  by 


G(r,z,r',z’)  = 


.  Wls  • £ 


2  2 

+  (z  -  z’) 


4tr 


(z  -  z')2 


^|r  -  r'|2  + 

The  V  •  operator  in  the  second  term  is  taken  inside  to  give 


(73) 


E  (r,z)  =  -jn 


("k  f  f  J  G  dt'dz’  +  £  V  /  (  V  •  J 

L  *C  z’ _  k'Jcl~  ~ 


G  dt'dz' 


(74) 


The  transform  of  Equation  (74)  is 


£S(r,k)  = 


V37 


if  h 


J ) 


C  z' 


-jk  z 

G  e  2  dz  I  dt ' dz ' 


-jk  z 
G  e  dz  /  dt 'dz 1 


(75) 


To  do  the  first  integral,  note  that  G  satisfies 


V2C  +  k2G  -  -5(r  -  r’) 


(76) 


This  is  rewritten  in  the  form 


7?G  +  ^-|  +  k2G  -  -6(R)  5(2  -  z’) 

32 


(77) 


where  R  *  (x  -  x' )  x  +  (y  -  y')  The  transform  of  Equation  (77)  with  respect 
to  (z  -  z')  is 


(72  +  k2)  G  =  -5  (R) 


(78) 


where 


G  = 


00 

V57  L 


-jk  (z  -  z’) 


G  e 


d(z  -  z’) 


(79) 


and 


2  2  2 
k  -  k  -  k 
t  z 


(80) 
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The  solution  to  Equation  (78)  is 


g  =  |jh0(2)  (kl*l> 


Thus  we  have 


-jkz’? 

G  e  dz  =  e  G 


The  second  integral  in  Equation  (75)  is  aided  by  the  identity 


-jk2z  -jk  z 1  ~  ~  ~ 

Ge  dz  =  e  G  -  j  k^  G  z 


Thus  we  have,  after  writirg  J_  =  J_c  +  £  Jz> 


•  'VS  kfclQ'*k 


-jkzz'  , 

J  )  e  G  dt'dz' 

z 


CD 

iff 

*  r  * 


“  i  k  z  *  A  A 

(7'  •  J)  e  7  G  -  j  k  G  z  dt’dz' 

—  —  —  t  z  — 


Now  let 


?U\k  ) 


CD 

-1—  (  JU’.z’ 

s/2^  y  „ 


-jkzz’ 

)  e  dz 


and 


J(t\z') 


(86) 


Expanding  che  divergence  of  J  gives 


T  •  I  (t*,z')  =  j  kz  J  (t’.z')  +  yp-  J  (t\z’) 


(87) 


Substituting  Equation  (87)  into  Equation  (84),  one  obtains 


£S(r,k  ) 


■  -J”  [k  J  ft a  +  i  j 


(t'.k’) 


+  j  kz  £  (t'.k^j  •  ?t  G  dt’j 


(88) 


+ 


3 

3t ' 


(-j  kz)  G 


dt' 


Thus  putting  Equation  (71)  into  matrix  form  would  vield  the  following: 


c. 

ZZ 

z  1 

2  t 

1 

N 

s 

I 

o , 

N 

- 1 

m'“tz 

r>: 

1 

r* 

1 _ 

V 

t 

L  J 

l  84) 


where  the  elements  of  the  submatrices,  Z,  are  found  from  Equation  |K8). 
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2.8  CURRENT  DISTRIBUTION  ON  CONDUCTOR  DUE  TO  IMPRESSED 
LONGITUDINAL  CURRENT 

In  che  preceding  sections,  formulations  were  presented  for  determin¬ 
ing  the  induced  surface  current  distribution  on  perfect  conducting  cylinders 
of  arbitrary  cross  section  and  infinite  length.  The  excitation  was  taken  to 
be  a  uniform  plane  wave.  Here  we  consider  a  different  type  of  excitation, 
namely,  that  of  a  steady-state  current  which  flows  axially  along  the  cylinder. 
If  the  cylinder  contour  C  is  not  circular,  the  current  density  on  C  is  dis¬ 
tributed  around  the  contour  due  to  inductive  effects.  It  is  this  current 
redistribution  which  is  solved  for  here  by  an  integral  equation  formulation. 

The  perfect  electric  conductor  C  is  infinite  in  the  z  direction  and 
carries  a  total  current  I  amp.  Let  the  surface  current  density  on  C  be 
denoted  by  J amp/m.  Then 

J  dc  =  I  (90) 

z 

where  dc  is  the  elemental  arc  length  on  C.  Since  I  is  independent  of  fre¬ 
quency,  there  is  no  electric  field.  Thus  the  magnetic  field  satisfies  the 
equations : 

V  •  _H  ■  0  everywhere  (91) 

IJ  on  C 
— z 

(92) 

0  in  R  and  R. 

—  01 

Since  C  is  a  perfect  conductor,  H=  0  in  the  internal  region,  R^,  shown  in 
Fig.  2-3.  For  the  region,  R  ,  H  may  be  written  as  [from  Equation  (92)] 


where  A  satisfies 


V  x  V  *  A  =  j  =  V (V  •  A)  -  7“  A 


(94) 


A  is  in  the  z  direction  only  and  independent  of  z  so 


7"A  =  -J 

z  z 


(95) 


The  boundary  condition  that  H  satisfies  is  [from  Equation  (91) 


n  •  H  =  0  on  C 


(96) 


This  is  rewritten  in  terms  of  A  as 

z 


3A 

_ z 

3c 


=  0  on  C 


(97) 


where  we  have  used  the  fact  that  n  x  t  =  z.  The  solution  to  Equation  (95)  may 
be  written  as 


=  -  2,  / 

*  n 


J  (c  ' )  In  |  r  -  r '  '  dc ' 


(98) 


where  r  is  a  point  in  Rq  and  r_'  =  _r'(c')  is  a  point  on  C.  Now  Equation  (97) 

implies  that  A^  =  constant  on  C.  Thus  the  integral  equation  that  J satis¬ 
fies  is 


-u 


J  (c1)  In  r  -  r'  dc'  =  K 

z  — 


(99) 


subject  to  the  constraint  of  Equation-  (90).  K  is  a  constant  which  depends  on 
the  geometry  of  C.  For  a  circle  of  radius  a. 


K  = 


2~t  a"  In  2a 


(.100) 
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Equation  (99)  may  be  solved  by  choosing  K  =  -1/2".  Proceeding  as  in  the  pre¬ 
vious  sections,  we  break  C  up  into  N  subsections.  Over  each  subsection,  J 

z 

is  assumed  constant.  This  is  equivalent  to  using  the  pulse  basis  defined 
by  Equation  (15).  Thus  we  have 


N 

J,  =  T  a.  P.(t)  (101) 

2  ■  i  J  1 

j=l  J  J 

where  ot.  are  unknown  coefficients.  This  is  substituted  into  Equation  (99)  to 
J 

obtain 


N 


(102) 


For  computational  simplicity,  a  point-matching  procedure  is  used  where  the 
impulse  functions 


11  if  £  is  on  A. 

(103) 

0  if  _r  is  elsewhere  on  C 

are  used  for  testing.  The  resulting  matrix  equation  is  given  by 


[Ml  a  =  K 


(109 ) 


where 
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ij 


A. 
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A. 


dt 


(105) 


2-30 


Equation  (105)  can  be  evaluated  analytically  and  is  given  by 
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if  i  =  j 


M.  . 
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(106) 


^  if  i  ^  j 

where  the  following  notation  is  used: 


r . 
— x 


vector  from  origin  to  midpoint  of  A 
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SECTION  3 


THIN  SHELLS  OF  ARBITRARY  CROSS  SECTION 
AND  FINITE  CONDUCTIVITY 


Two  approximate  methods  are  presented  here  for  computer-aided  analy¬ 
sis  of  field  penetration  into  cylindrical  shells.  The  first  method  effectively 
replaces  the  shell  by  an  equivalent  impedance  sheet  boundary  condition  which 
results  in  a  modified  E-field  integral  equation.  The  second  method  utilizes  a 
transmission  line  analysis  to  derive  a  surface  load  impedance  to  be  used  in  a 
loaded  body  E-field  integral  equation.  The  E-field  operator  developed  in 
Section  2  for  both  polarizations  is  used.  The  adjective  "thin",  as  used  here, 
means  the  shell  thickness  is  small  with  respect  to  a  wavelength  in  the  sur- 
rouding  medium  but  may  be  appreciable  with  respect  to  the  shell  material 
wavelength. 


3.1  INTRODUCTION 

The  primary  purpose  of  this  section  is  to  develop  some  approximate 
techniques  for  computing  the  electromagnetic  scattering  and  penetration  prop¬ 
erties  of  two-dimensional  shells  of  an  arbitrary  cross  section  and  having 
finite  conductivity.  Quantities  of  interest  are  thus  scattering  cross  section 
and  near  fields  inside  the  shell.  The  latter  are  characterized  by  the  "shield¬ 
ing  effectiveness"  of  the  shell  which  is  defined  here  as  [1] 


SE 


20  log 


(dB) 


(1) 


NS  S 

where  F  and  _F  are  fields  computed  at  a  point  without  and  with  the  shell 
present,  respectively. 


The  original  problem  is  shown  in  Fig.  3-1  where  a  plane  wave  illumi¬ 
nates  a  shell  of  thickness  d.  The  shell  is  made  up  of  material  with  consti¬ 
tutive  parameters  u  ,  £,  a  and  the  surrounding  material  is  free  space  (y  ,  e  ). 

o  0  0 


f 
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3 »  u  •  z. 
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Figure  3-1.  Original  Problem:  Plane  Wave  Illuminating 
a  Shell  of  Uniform  Thickness  d 

The  conductivity  0  may  also  be  a  function  of  position  in  the  shell.  The 
thickness  d  is  assumed  to  be  much  less  than  the  wavelength  of  free  space,  A  . 
Thus,  as  far  as  the  surrounding  medium  is  concerned,  the  shell  may  be  replaced 
by  a  single  contour,  C.  This  contour  is  further  approximated  by  a  finite 
number  of  straight  line  segments  AC^  for  i  =  1,  2,  ...  ,  NC.  This  is  shown 
in  Fig.  3-2.  The  original  properties  of  the  shell  are  accounted  for  by 
assigning  to  each  segment  ACk  a  value  of  d  and  O.  Each  line  segment  has 
length  A.  and  unit  vectors  t.  and  n.  such  that  t.  x  n.  =  z.  The  excitation 
consists  of  two  types  of  plane  waves,  each  to  be  considered  separately.  These 
are  the  TE  case  (z  component  of  magnetic  field  only)  and  the  TM  case  (z  com¬ 
ponent  of  electric  field  only). 

A  general  formulation  of  the  problem  in  Fig.  3-1  requires  the  use 
of  equivalent  electric  and  magnetic  currents  on  the  inner  and  outer  surfaces 
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Figure  3-2.  Approximation  of  Shell  by  Contour  C 

of  the  shell.  This  is  given  in  [2]  and  will  not  be  discussed  here.  Instead, 
some  approximate  formulations  will  be  developed  which  are  valid  for  certain 
types  of  shells.  As  a  starting  point,  the  shell  material  is  assumed  to  be  a 
fairly  good  conductor.  If  d  is  also  much  less  than  the  wavelength  in  the 
shell,  then  an  impedance  sheet  approximation  may  be  used  [3].  The  deriva¬ 

tion  of  this  formulation  is  summarized  in  Subsection  3.2  for  use  in  computer 
program  1.  If  the  frequency  is  higher,  so  that  d  is  then  comparable  to 
the  shell  material  may  be  assumed  to  support  traveling  waves.  Here  a  trans¬ 
mission  line  analysis  is  presented  in  Subsection  3.3  for  use  in  computer  pro¬ 
gram  2.  Lastly,  if  the  shell  is  circular  then  an  infinite  series  solution  is 
possible  using  Bessel  functions  [4].  This  is  presetted  in  Section  4  for 
use  in  computer  program  3.  The  desired  quantity  in  all  three  formulations 
is  the  field  at  points  interior  to  the  shell.  This  nay  be  expressed  as  an 
integral  over  electric  and  magnetic  currents  on  C  and  procedures  for  this 
computation  are  given  in  Subsection  3.4.  Descriptions  of  the  computer  pro¬ 
grams  as  well  as  sample  input/output  data  are  given  in  Appendix  B. 
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3.2  IMPEDANCE  SHEET  APPROXIMATION 


The  total  field  everywhere  in  Fig.  3-1  is  the  sum  of  the  incident 
i  i  S  S 

field  (E  ,  H  )  and  a  secondary  field  (E  ,  _H  )  due  to  the  presence  of  the  shell. 
This  secondary  field  may  be  generated  by  an  equivalent  electric  polarization 
current  which  effectively  replaces  the  shell.  This  current  is  given  by  [5] 


J 


[joj(e  -  e  )  +  a]  £  in  S 
0  outside  S 


(2) 


where  E_  is  the  total  electric  field  (E1  +  E^) .  Equation  (2)  may  be  rewritten 
as 


-E  (Jj  + 


j(jj(£ 


£o) 


+  o 


in  S 


(3) 


where  £  is  an  electric  field  operator  defined  by  Equation  (8)  of  [6].  If  the 
shell  thickness  d  is  much  less  than  the  wavelength  in  region  b,  then  one 

may  approximate  Equation  (3)  by  specializing  it  to  the  contour  C  in  Fig.  3-2 
and  replacing  the  volume  current  J_  with  a  surface  current  J<..  One  then  obtains 
the  loaded  body  equation  [7] 


- -  ES  (J  )  +  ZT  J 

n  — t  L  — c 

o 


on  C 


(4) 


where  the  factor  of  k  /r\  has  been  multiplied  through  for  later  convenience. 

o  o 

The  subscript  t  denotes  tangential  component  evaluated  on  C. 


The  normalized  impedance  load  is  given  by 


Z 


L 


(5) 
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which,  if  £ 


£  ,  reduces  to 


Z 


1 


_ o 

at 


(6) 


This  is  the  low  frequency  limit  used  in  [1].  Equation  (4)  is  solved  by  a 
moment  method  procedure  as  in  Section  2  [6]  where  the  matrix  equations  are 
written  as 


for  the  TM  case  and 


(7) 


(8) 


q  h  -+ih 

for  the  TE  case.  The  matrices  [Z  ],  (Z  ]  and  vectors  V  .  V  are  exactly 

the  same  as  those  in  Section  2.  The  vectors  I  and  I  contain  the  coeffi- 

cients  of  expansion  for  J  in  the  TM  and  TE  cases,  respectively,  which  is 

^  e 

the  same  as  that  used  for  J  in  [6].  For  the  TM  case,  the  elements  of  Z 

— c 

are  given  by 


,  k  A 
/  o  m 


(ZL> 


mn 


if  in  1  n 


if  m  ^  n 


where  8  is  defined  by 


8  = 


(9) 


(10) 
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For  the  TE  case,  the  elements  of  2^  are  given  by 
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T  T  (Am-l  +  V  if  m  =  n 


k  A 


«J>  '  r 


mn 


k  A 
q  n 

36 


V. 


if  m  =  n-1 

or  n+1  -  MC 

if  id  =  n+1 

or  n-1  +  NC 


(11) 


Note  that  in  this  formulation  J_  is  assumed  tangential  to  C.  Any 
normal  component  which  the  actual  polarization  current  may  have  has  been  neg¬ 
lected.  This  is  probably  acceptable  for  the  TM  case  since  J  is  j  directed. 

For  the  TE  case,  however,  this  assumption  is  no  good  unless  k^d  <<  1  and  even 
then  depends  upon  the  incident  field.  For  example,  the  configuration  in 
Fig.  3-3  would  produce  erroneous  results  by  the  above  formulation.  A  more 
accurate  solution  could  be  obtained  by  allowing  both  components  of  the  polari¬ 
zation  current  [8,9]. 


1 


Figure  3-3.  Incident  Field  Causing  Normal  Component 
of  Polarization  Current 
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3.3  TRAVELING  WAVE  APPROXIMATION 


The  problem  in  Fig.  3-1  may  also  be  looked  at  as  a  three-region 
problem,  where  equivalent  electric  and  magnetic  currents  are  assumed  to  exist 
on  surfaces  Cq  and  C^.  [2].  This  formulation  will  not  be  presented  here,  but 
for  the  purposes  of  discussion,  let  "a",  "b",  and  "c"  denote  the  regions  out¬ 
side  Cq,  between  Cq  and  C^,  and  inside  C^,  respectively.  The  total  field  in 
region  "a"  now  may  be  expressed  as  the  sum  of  the  incident  field  and  a  sec¬ 
ondary  field  arising  from  electric  and  magnetic  current  sources  on  Cq.  Now 
if  the  shell  is  a  good  conductor,  then  the  magnetic  current  on  Cq  will  be 
negligible.  Secondly,  if  the  shell  surface  has  no  abrupt  changes  in  curvature 
[10]  one  may  assume  that  an  impedance  relationship  exists  between  the  total 
tangential  component  of  the  electric  field  in  region  "a"  and  the  electric  cur¬ 
rent  on  Cq.  Again,  since  kQd  <<  1,  we  replace  Cq  and  bv  C  in  Fig.  3-2.  The 
condition  that  |k^d|<<  1  as  in  Subsection  3.2  need  not  apply  here.  Hence,  we 
write 

Ea  =  Z  J  on  C  (12) 

— t  L— o 

where  _Ea  is  the  total  electric  field  in  region  a.  Equation  (12)  is  rewritten 
as 


-Ea(J  )  +  ZT  (J  )  =  E*  on  C  (13) 

— t  — o  L  — o  — t 

which  is  again  the  loaded  body  equation  of  Subsection  3.2. 

The  load  impedance,  Z^,  this  time  will  be  determined  by  assuming 
that,  inside  region  "b",  each  subsection  of  C  appears  locally  planar.  Travel¬ 
ing  waves  are  then  assumed  to  exist  in  region  "b"  which  reflect  the  impedance 
seen  at  looking  into  region  "c"  back  to  region  "a".  Standard  transmission 
line  techniques  may  thus  be  used  to  obtain  Z 

First  consider  the  infinite  slab  shown  in  Fig.  3-4  where  a  local 
( u , v , w)  coordinate  system  is  used.  The  electric  surface  current  exists 
everywhere  on  the  plane  u  =  0  and  is  constant  over  all  v.  This  gives  rise  to 
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plane  waves  in  all  three  regions  which  have  no  u  components.  The  electric 
field  is  written  as 


jk  u 

a  J  a 

E  =  A  t  e 


(-  u  traveling  plane  waves) 


b  »  -jk^u  ' 

=  B  t_  e  +Ct_e  C+  u  traveling  plane  waves) 


,  ,  -jk  (u-d) 

E  ^  -  D  t  e  c 


(+  u  traveling  plane  waves) 


The  transverse  unit  vector  t_  lies  in  the  v-w  plane.  The  magnetic 
field  is  obtained  from  the  Maxwell  curl  equation 


=  -jk  n  H 


which  has  been  specialized  to  +  u  traveling  plane  waves.  Thus  we  obtain 


>  jk  u 

,,a  1  J  a 

11  *  —  (n  *  t)  e 

a 


ceJV 


-jk  (u-d) 

HC  =■  S-J-t  o  e  c 
-t  n 


ihe  boundarv  conditions  which  must  be  satisfied  are 


h  <  [  HJ  -  H3  ]  =  .1  at  u  =  0 
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Ea  =  Eb  at  u  =  0  (22) 

— t  — t 

at  u  =  d  (?'’) 

at  u  =  d  (24) 

Now,  si'  e  H  =  -n  x  n  x  jtl,  Equation  (20)  is  rewritten  as 

Ha  -  =  n  x  J  at  u  =  0  (25) 

— t  — t  — o 

Solving  Equations  (22)  through  (25)  simultaneously  for  the  coeffi¬ 
cients  A,  B,  C,  and  D,  one  obtains 


where  A  is  given  by 

A  =  2(1  +  r.,^)  cos  k^d  +  2  j  (r^  +  r^)  sin  k^d  (30) 
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p 


and  r  =  n  /n  ,  r„  =  n,  /n  • 
1  a  b  Z  be 


In  Che  above,  J  is  in  the  c  direction  so  that 


the  load  impedance  is  determined  by  the  ratio 


J  =  J  t.  In  Che  actual  problem,  E  and  H  must  be  continuous  across  'J  so 

—  t  —  t  o 

e  ratio 

cos  k^d  +  j  r0  sin  k^d 


_b 

ii  x  -5.J- 


Hb 
— t 


u=0 


j  sin  k^d  +  r,,  cos  k^d 


(31) 


The  following  limiting  cases  of  Equation  (31)  may  be  used  when 
applicable : 


r 


1  +  J  r2  V 
r2  +  j  V 


1_ 

ad 


-j  nb  cot  i^d 


V. 


IS 


d  •*  0 


|  k^d  ]  ■*  0  and  |  j  <<  hc 


>V  <<  \ 


;  k^d  i  -*•  00  and  |  |  <<  hc 


(32) 


The  tangential  components  of  field  at  surface  are  given  by  the 
expressions 


and 


EC 

-t 


2  ”  J  -ik  (u-d) 

a  ot  2  J  c 
- : -  t  e 


(33) 


H 

— t 


'  ~  2  n  J  „ 

(n  x  t)  a  ot 


-jk  (u-d) 
J  c 
e 


(31) 
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<  Cl 


These  fields  may  be  thought  of  as  arising  from  surface  electric  and 
magnetic  currents  at  u  =  d  defined  by 


n 


(n  x  J  ) 


— o  (1  +  r^r0)  cos  k^d  +  j (r^  +  r0)  sin  k^c 


(35) 


and 


j  =  j  — =• 


— 1  — < o  nc  (1  +  r^r0)  cos  k^d  +  j ( r ^  +  r^)  sin  k^c 


;ain,  as  n,  <<  n  or 


(36) 


4 


(n  x  J  ) 

— o 


'b 

j  sin  k^d 


(37) 


4 


J 

— o 


nc  j  sin  l^d 


(38) 


In  the  above,  we  have  assumed  that  the  field  in  region  c  is  a  plane 

wave.  This,  of  course,  is  not  exactly  true  so  the  approximation  will  probably 

fail  unless  region  c  is  electrically  large.  Note  that  the  normal  unit  vector 

used  here  is  opposite  to  that  used  in  Figs.  3-1  and  3-2. 
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3.4  COMPUTATION  OF  FIELDS  DUE  TO  TWO-DIMENSIONAL  CURRENT 
DISTRIBUTIONS 

Expressions  are  presented  here  for  the  field  at  points  interior  to 
tne  shell,  denoted  by  region  c  in  [2].  For  the  formulation  of  Subsection  3.2, 

the  total  field  in  region  c  is  due  to  an  incident  field  plus  a  secondary  field 

which  is  caused  by  a  two-dimensional  electric  current  distribution  on  C.  For 

the  formulation  of  Subsection  3.3,  the  total  field  in  region  c  is  due  to  elec¬ 

tric  and  magnetic  currents  on  C.  In  both  cases  these  current  distributions 
radiate  in  unbounded  space  filled  with  Uc>  £c .  Hence,  we  represent  the  fields 
by  a  potential  integral  formula  [14]. 

The  actual  fields  computed  are  the  z  components  of  electric  field 
in  the  TM  case  and  magnetic  field  in  the  TE  case.  These  are  written  as 

i 
E 

—  -  z*  [jkA  +  VxF]  (39) 

n  —  J  c—  —  — 

c 

and 

HZ  =  Hz  +  i  *  tY.  X  a  —  j  k^F]  (40) 

where  the  electric  and  magnetic  vector  potentials  are  defined  by 

'  (t')|)  dt’  (41) 

F  =  —  /  M(t')  H  (2)  (k  |  r  -  r'  (t’)|)  dt’  (42) 

c  •/c 

In  the  above,  the  electric  field  and  magnetic  current  have  been 
normalized  by  o  for  computational  convenience.  Both  terms  in  Equations  (39) 
and  (40)  are  used  in  the  formulation  of  Subsection  3.2  and  the  last  terms 
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only  are  used  for  the  formulation  of  Subsection  3.3.  Equations  (39)  and  (40) 
at  a  point  £  in  region  c  where  jr'(t')  is  a  point  on  C  where  arc  length  is  para¬ 
metrically  expressed  in  terms  of  t'.  The  electric  and  magnetic  currents  are 
expanded  as 


J 


NC 

2 


ijj  T  (t) 


I  JL<t> 

n  — n 


(43) 


NC  , 

2  Vn  4(t) 
n=l 


+  V_  T  (t) 


(44) 


where  P  and  T  are  defined  by  Equations  (14)  and  (13)  of  [6],  I  and  V  are 
— n  — n  n  n 

6  S 

complex  coefficients  for  the  TE  case  and  I  and  are  complex  coefficients 
for  the  TM  case.  Equations  (39)  and  (40)  may  be  conveniently  rewritten  in 
terms  of  near-field  measurement  vectors  as 


—  +  Q  f  +  p  V 


(45) 


c  i  -*Vi  -*■  -►h 

H  =  H  -  P  I  +  0  V 
z  z 


(46) 


where  the  tilda  (~)  denotes  transpose. 

3.4.1  Formulas  for  Near-Field  Measurement  Vector,  0 

Each  element  of  Q  actually  represents  the  electric  (magnetic) 
field  due  to  a  z  directed  electric  (magnetic)  current  of  amplitude  1/HC  (h  ) 
with  a  pulse  function  distribution  on  subinterval,  AC^.  The  field  point  is 
denoted  by  r_  and  r_'  denotes  a  point  on  AC^.  Thus  one  may  write 


/ 


n+1 


P  (t  ’  )  H 
n  o 


(2) 


(kc  |r  -  r*  (t')|)  dt' 


(47) 
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Let  r  be  a  vector  from  the  origin  to  the  midpoint  of  AC  and  define  R  (t')  as 
— n  n  -n 


R  (t’) 
— n 


r  -  r 


t  t 

2  — n 


(48) 


which  is  shown  in  Fig.  3-5,  where  -1  t  <_  1.  Then  Equation  (47)  may  be  trans¬ 
formed  to 


k  A 
c 


>  A 

-1 


(2) 


(k  R  (t ' )  )  dt' 
c  — n 


(49) 


The  integrand  of  Equation  (49)  becomes  singular  when  ;R  (t')|  -*■  0. 
"  — n 

To  remedy  the  numerical  difficulty  encountered  when  this  happens,  we  rewrite 
Equation  (49)  as 


■  -  ^  [%m  (tc 


|Rn(f )  | )  dt’ 


+  2i 


£n 


Yk  ]  R  ( t  ’ ) 


dt’ 


'sc; 


k  A  r  T 

_£ _ 2.  /  2j_ 

8  J 


Yk  ! R  ( t ’ ) 


dt T 


\ 


whenever  |r  (t')|  <  e  for  some  small  number  c  >  0  and  subinterval  A  .  The 
first  integral  can  be  done  accurately  by  a  quadrature  rule  as  long  as  the 
integrand  is  never  evaluated  exactly  where  (R^(t')|  =  0.  The  second  inte¬ 
gral  can  be  done  analytically  and  the  following  substitutions  are  made: 

_r  =  x  x  +  y  £ 


t  =  t  x  +  t  y 
— n  nx  —  nv 


3-15 


in 


Figure 


Origin 


3-5.  Geometry  Relating  to 
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r  *  x  x  +  y  v 
— n  n  —  yn  — 


a  =  (A  /2) ‘ 
n  n 


b  =  -A  [(x  -  x  )  t  +  (y  -  v  )  t  ] 
n  n  n  nx  w  -  n  nv 


cn  =  -  *  >2  +  (y  -  yJ2 


D  =  b  -  4  a  c  =  2A  ^  (x  -  x  )  (v  -  y  )  t  t 
n  n  n  n  n  n  '  '  n  nx  nv 


Then  |R  (t')|  becomes 
— n 


V‘'>l  ■  |/v’2  +  V  + 


and  the  second  term  in  Equation  (50)  is  written  as 


k  A 
c  n  2 

4  TT 


£n 


(a  +  b  +  c  )  (a  -  b  +  c) 
n  n  n  n  n  n 


4  a 


£n 


a  +  b  +  c 

n _ n _ n 

a  -  b  +  c 

n  n  n 


(51) 


.  b  -2a. 
-In  n 


X/D~  /  .  2  a  +  b 

,  V  n  -1  n  n 
+  - -  I  tan - — — -  -  tan 

*n  \  \/\ 

for  D  <0.  If  D  >  0,  D  is  replaced  bv  -D  and  tan  ^  is  replaced  bv  tanh  ^ . 
n  n  n  r  n 
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« 


3.4.2  Formulas  for  Near-Field  Measurement  Vector,  P 

The  element  P^  represents  the  magnetic  (electric)  field  due  to  a 

-.t  (_t)-directed  electric  (magnetic)  current  of  unit  amplitude  with  a  triangle 

function  distribution  over  the  interval  AC  ,  V  AC  .  Thus  we  have 

n-1  n 


l 


(2) 


-  77  j  •  7  X  i  T  (t')  H  (k  lr  -  r'(t')  )  dt' 
—  f  — n  o  c  —  — 


(52) 


^  f  Tn(t,)  *  '  ^(C>)  H  (2) 

4 j  J  ^R(t’)n  H1 


(k  iR(t’)  )  dt' 
c  — 


where  R(t')  =  r  -  r'(t').  This  may  be  rewritten  as 

.1  /I  .  t' 


1_ 


k  A 
c 


(2  +r)i-i-^-i(t’> 


H  (2)  (k  |R  (t ' ) ! )  dt' 
1  c  —n-1 


k  1 

rL(j- 

f) 

n 

•  R  (t* ) 

c  n 

/  \  2 

2  ) 

— n 

— n 

R  (t  •)  ! 


Hi(2)  (kc  ! ( t) 1 )  dt’ 


-1 


(53) 


The  integrand  of  Equation  (53)  becomes  singular  when  1 -*■  0.  The 
singularity  is  integrable,  however,  and  after  a  similar  manipulation  to  that 
done  in  Subsection  3.4,1  one  obtains 


P 

n 


P  ,  +  S  ,  + 
n-1  n-1 


+  S 


( 54  ) 
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SECTION  4 


SHELLS  OF  CIRCULAR  CROSS  SECTION  AND 
FINITE  CONDUCTIVITY 

When  the  shell  has  a  circular  cross  section,  the  Helmholtz  wave 
equation  is  separable  in  cylindrical  coordinates.  An  expansion  of  the  fields 
in  each  region  is  then  possible  in  terms  of  the  solutions  to  the  homogeneous 
wave  equation.  The  unknown  coefficients  of  expansion  are  found  by  ensuring 
continuity  of  tangential  fields  across  the  inside  and  outside  shell  surfaces. 

4.1  INTRODUCTION 

An  exact-series  solution  [1,  2]  for  the  penetration  fields  inside 
circular  shells  is  presented  here.  This  is  useful  in  testing  the  approximate 
solutions  developed  in  Section  3.  This  solution  is  especially  useful  in  de¬ 
veloping  low-frequency  approximations  to  the  interior  penetration  fields,  i.e., 
when  the  overall  dimension  of  the  shell  cross  section  is  small  with  respect  to 
a  wavelength.  The  interior  fields  are  shown  to  be  uniform  in  this  limit  and 
a  convenient  equivalent  circuit  model  is  valid  for  the  coupling  mechanism. 

4.2  GENERAL  SOLUTION 
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Plane  Wave  Incident  Upon  a  Circular  Shell 


The  z  component  of  magnetic  field  in  each  region  is  expressed 


as 


<  -  2  j"  -r  w 


eJ  +  H 


-  s 


.n 


“lE  W>  *  C  <V> 


Jno 


(2) 


H°  =  V  jn  dTt  J  (k  r)  eJ 
z  J  n  n  o 


where  H  (x)  =  J  ( x )  -  j  Y  (x). 
n  n  n 


The  <>  component  of  electric  field  is  obtained  from 


Hence 


n  3H 

n  _ z 

jk  3r 


(3) 


<  =  i  ^2  Jn 


TF 

a  H 1 (k  r)  +  J'(k  r) 
n  n  o  no 


Jn<|> 


=  j  %2  j 


.  n 


bnE  j;v>  +cnE-Y;(V> 


JntJ) 


(4) 


V  .n  jTE  , ,  ,,  , 

=  1  H  X  j  d  J  (k  r)  e 
J  o  J  n  no 


jn6 


/ 
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Matching  tangential  components  of  E  and  H  at  r  =  a^,  leads  to 


H  (k  a  )  -J  (k,  a  )  -Y  (Ic  a  )  0 

noo  n  d  o  n  d  o 

°  Jn(kbal}  Yn(kbal}  -Jn(koal 

H;(Vo)  -5  j;<v0>  Y;(Vo>  0 


-J  (k  a  ) 
n  oo 


)  b 


0  £  J’  (k.  a. ) 

n  d  1 


Y'(k.a.)  -J '  (k  a, )  d 
n  p  1  n  o  1  n 


•J’  (k  a  ) 
n  o  o 


where  £  =  n,  /h  .  This  svstem  has  the  determinant: 
bo 


where 


H  (k  a  ) 
Ic^a^  In  oo 


J ' (k  a, 
n  o  1 


)  T  ..  +  £  J  (k  a.  )  T  .  ' 
nil  n  o  1  nl2 


H' (k  a  )  f  4  J'(k  a,)  T  +  J  (k  a.)  T  „„ 
n  o  o  ^  n  ol  n21  n  o  1  n^^ 


Tnll 

’Vi  ; 
2  1 

VVi1 

y;(Vo> 

-  j;(kbao) 

Yn(kbal} 

rrk  a.  | 

Tnl2 

V  1 

“  .■) 

j;Vo} 

Y;(kbai> 

-  Jn(kbal} 

W 

T 

n21 

’Vl 

2 

LJn(kbal} 

Yn(kbao) 

-  Jn(kbao} 

Yn  ^kbal> 

T 

n22 

'Vl  1 
2 

[Jn(kbao) 

y;<Vi> 

-  Jn(kbal) 

Wo’ 

Solving  for  a 


TE 

n 


and  d 


TE 
n  ’ 


one  obtains 


TE 


.  .TE 

Vl  An 


-  J  (k  a  ) 
n  o  o 


J'(ka.)  T„  +;  J  (k  a.)T  ,, 
n  o  1  nil  n  o  1  nlz 


(8) 


+  J '  (k  a  ) 
n  oo 


TJ'fka.)  T  +  J  (k  a.)  T  ,, 
c,  n  o  1  n21  n  o  1  n22 


,TE 


,TE  2  . 

2.  t  k  a  k,  a. 
n  o  o  d  1 


(9) 


These  coefficients  allow  one  to  compute  the  total  fields  external  to  and  in¬ 
ternal  to  the  shell. 


In  a  similar  procedure  for  TM  excitation,  one  obtains: 


z 


jk  x 
o 


.n  T 
J  J 


n(knr) 

n  o 


jny 


(10) 


The  z  component  of  electric  field  in  each  region  is  expressed  as 


E 


a 

z 


.n 


TM 

a  H  (k  r)  +  J  (k  r) 
n  no  no 


ein6 


Z 


2  JB 


b™  w> 


TM 

4-  c  V  (k.r) 
n  no 


ejn(J) 


(11) 


E 


c 

z 


2 


TM 

n 


J  (k  r) 
n  o 


ejnd 


The  i  component  of  magnetic  field  is  obtained  from: 


H 

5 


1 

jkn 


U2) 
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Hence : 


$  jn 


^  ?» J”  [■“  ")'V>  *  w>] eJno 
2  j”  [»™  w  * c”  w]  •3no 


b  — » 


1  V  ^  T  I  /,  1 

t—  >  j  d  J  (k  r)  e 
jn  *“L  J  n  n  o 


jna> 


(13) 


Again,  matching  tangential  components  of  _E  and  H  at  r  =  a,  a  leads  to 


Wo>  -  Wo* 

-  '.V.1 

O 

_l 

TM 

a 

n 

-  J  (k  a  ) 
n  oo 

0  Wi> 

wc 

-  J  (k  a.) 
n  o  1 

b™ 

n 

0 

“;<Vo>  - 1 

- 1  K(W 

0 

TM 

c 

n 

-  J’(k  a  ) 
n  oo 

.  °  I  j;<w 

r  KlW 

- 

n  o  i 

,TM 

a 

n 

0 

•—  *• 

(14) 

where  again  £  =  ri,  /ri  •  Comparison  with  Equation  (5)  leads  to  the  result  that 
TM  TM  Til  0 

,  a^  ,  and  d^  are  given  by  Equations  (6),  (8),  and  (9)  with  %  replaced 
by  1 /%. 
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4.3  LOW  FREQUENCY  APPROXIMATION 


Of  particular  interest  is  the  case  when  the  shell  geometrv  of 
Fig.  4-1  is  "quasistatic".  This  is  the  case  when  k  a  <<  1.  After  using 
the  first  few  terms  in  the  series  for  J^(k^r)  and  neglecting  terms  of  order 
(kQr)  and  higher,  one  obtains  the  following  formulas  for  the  fields  interior 
to  the  shell: 


,TE  .  ,TE  ,  . 

d  +  id,  k  r  cos  f 
o  1  o 


,TE  V  ,TE  .  TE  kor 

_Jno  do  T  “  nodl  C0S  *  "  nod2  ~2~  C0S  2? 


n  3H 

O  Z 


jk  r  36 
o 


k  r 

-  nQd^E  sin  4>  +  jn0d2E  ~ f~  sin 


,TM  .  .TM  . 
d  +  jd.  k  r  cos 
o  1  o 


(15) 


k  r  d™  d™  k  r 

j  ,TM  q  1  .  ,  2  o  » . 

n  do  —  +  —  cos  *  +  n - r  cos  2^ 

o  o  o 


dE" 


ik  n  r  3<t> 
o  o 


,TM 


k  r 

.  .  ,TM  o  ... 

sm  6  -  J —  d?  sin  26 

o 
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We  see  rrom  the  above  that,  when  k^r  <<  1,  the  fields  interior  to  the  shell 
are  uniform  and  are  given  by: 


H°  =  dTE 
z  o 


rc  _  „  ,TE 

E  =  -  n  d, 

v  o  1 


E  =  0 

x 


for  the  TE  case  [Fig.  4-2(a)]  and 


HC  -  i-d™ 
v  n  l 
o 


H  =  0 

x 


for  the  TM  case  [Fig.  4-2(b)].  It  is  now  necessary  to  examine  expressions  for 
the  coefficients  dQ,  d^,  and  hence  the  determinants  A.  To  do  this,  we  confine 
ourselves  to  shells  which  are  good  conductors,  i.e., 


*  d-i  wuo 


nb  - 


GOPG 


O  OQ 


* 


<»)  rt  Case 


Figure  4-2. 


<b)  TM  Case 
Interior  Fields 


Thus,  In  general,  ^bl>>  kQ  and  )?[«  1  and  hence  Ik^l  and  will  be 

quite  large.  The  following  limiting  forms  for  the  Bessel  functions  are  then 
quite  useful  [ 3 ] : 


As  x  -*■  0, 


Jo(x)  -  1  ~  (f ) 


Ji(x)  -  f  -  IT  (I) 


J'(x) 

o 


J[(x) 


*-*(!) 


H  (x) 
o 


Hlw  2  ttx 


H’(x) 

o 


_2L.11 

2  ttx 


H[(x) 


1  _  2 1  L  ix  1_ 

2  IT  \  n  2  +x2 


As  |x|  >>  1  and  | x|  >>  n,  we  have 


Jn(x)  "  /fe  C0S  (X  ■  T  -  l) 

Vx)  *  tfh  sin  (x  -  r  “  i) 
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The  T  matrix  elements  become: 
n 


nil 


nl2 


n21 


n22 


(21) 


where  d  =  a^  -  a^  is  the  shell  thickness.  To  simplify  the  results  we  assume 

that  d  <<  a  or  a.  and  hence  a  ~  a.  =  b.  This  can  be  taken  as  the  mean  radius 
o  1  o  1 

of  the  shell.  After  using  Equations  (19)  through  (21)  in  the  formulas  for  the 
TE  TM 

coefficients  d  and  d  we  obtain 


cos 


k  b 

V-  2T 


sin 


V 


2  k  b 
o 


2  k  b  cos  k.  d  +  —  sin 

o  d  t. 


V 


for  the  TE  case  and 


for  the  TM  case. 


(22'' 


(23) 


(24) 


(25) 
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If  the  thickness  d  is  such  that  |k^d|  <<  1,  then  the  following  low- 
frequency  behavior  for  the  interior  field  coefficients  is  obtained: 


,TE 

a 

o 


1 


_1 _ 

.  ,  ,  ad 

ju)Uob  — 


(26) 


1 


1 


ad  1 
2  jwe  b 


,TM 

a 


_ l__ 

1  +  jw  ^UQb  £n 


yk  b 
o 

2 


ad 


(27) 


(28) 


1 


1  + 


Vd 

2uQb 


j 


we  b 
o 

2 


ad 


(29) 


Using  Equations  (16)  and  (17),  the  electric  and  magnetic  shielding 
ratios  for  the  two  polarizations  can  be  written  as 


for  the  TE  case  and 


(30) 


(31) 


(32) 


(33) 
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for  the  TM  case.  The  low-frequency  behavior  of  shielding  effectiveness  can 
then  be  inferred  from  Equations  (26)  through  (29).  Of  particular  interest  is 
the  relationship  between  the  interior  electric  and  magnetic  fields  as  the 
frequency  u>  becomes  small.  For  the  TM  case  we  find  that 


v 


jj-*0 


Similarly,  for  the  TE  case  we  have 


Hc  uW) 


z 


,TE 
n  d, 
o  1 

,TE 


0 


(34) 


(35) 


4.4  tm  LINE  SOURCE  EXCITIATION 

Consider  the  thin  circular  shell  of  mean  radius  a  in  the  presence  of 

an  electric  line  source  as  shown  in  Fig.  4-3.  The  thickness  of  the  shell  d 

is  much  less  than  the  free  space  wavelength  and  also  Ag,  the  wavelength  inside 

the  shell  material.  If  the  latter  is  true  then  we  assume  an  impedance  rela- 

i  s 

tion  between  the  total  electric  field  ji  +  E^  and  the  electric  current  J_  at 
r  =  a  where: 

s 

ji  =  the  secondary  field  due  to  .J 
E  *  field  due  to  the  line  source  in  free  space 


J 


unknown  electric  polarization  current 


and 


2tt 

HS(J)  =  ij  7  x  /  aI(f)  ^  (k0!l-l'i)  ’  (38) 

*'o  ! 

Equation  (36)  will  be  solved  by  expanding  J_(d> 1 )  in  a  Fourier  series 
in  J '  and  specializing  the  result  to  the  quasi-static  case  (u)  -*■  0).  The 
incident  fields  are  given  by 


I  k.  n  ... 

E1  - - “2-  H(2)  (k  I  r  -  R  ') 

2  4  O  O  '  ““  — s 


I  k  v 

H1  = - -2-  H  2)  (k  |r  -  R  I) 

x  R  1  o — s 


I  k  (x  -  R  )  ,,, 

-^-R — -»1 ’  <koh-  S,l> 


(39) 


/ 


2  2 
(x  -  Rg)z  +  yz 


with  respect  to  the  coordinates  of  Fig.  4-3.  Equation  (36)  becomes: 


k  n 


f 


~  |  aJ(f)  H^2)  (kjr  -  r’|)  d<t>' 


(40) 


♦  !«♦>  '  -  t  '  ;lo2>  ».l£  -  s.i> 


f 
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s 


The  following  addition  theorem  is  now  useful  [3]: 


r(2) 


2  H (2) (k  r')  J  (k  r)  ejn(*  ’  r 
^  n  o  no 


<  r 


V  (koli-  I'l)  -  < 


(41) 


n  »(2>  /i,  J“0J>  -  $’> 


X  J  (k  r’)  fT  ;  (k  r)  e 
•“n  o  n  o 


r  >  r  ’ 


The  Fourier  series  representation  for  J_  is  written  as 


J(4>’) 


—  00 


(42) 


Substituting  Equations  (41)  and  (42)  into  Equation  (40)  and  using  the  orthogo¬ 
nality  of  e2n^  we  obtain: 


I  k  m 

-  HU;  (k  R  )  J  (k  a) 

_ 4  n _ os  no 

^  ( r)\ 

-  HU;  (k  a)  J  (k  a)  +  Z* 

2  n  o  n  o 


(43) 


where  Z1  =  Z  /r)  .  The  total  field  inside  the  shell  is  of  interest  and  is 
s  t  o 

s  i 

obtained  by  adding  _E  and  E_  for  r  <  a.  Thus  we  have 


ES  +  E1 
z  z 


£  e  h(2)  (k  R  )  J  (k  r) 


(41) 


k  an 

— r -  J  (k  a)  HU'  (k  a) 

_ 2 _ n  o _ n _ o _ 

k  a  ^ 

- -  H  (k  a)  J  (k  a)  +  Z' 

2  n  o  n  o 


-  1 
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This  gives  the  correct  result  as  1 '  —  0  for  a  perfect  conducting  shell  ( =  i 
for  n  *  0  and  2  for  n  >  0).  i'sing  the  formula 


r 


V  (k  r’) 

n+1  o 


J  (k  r)  e^n<1  r 


(ko^-  }  =  { 


(15; 


i  .ru  , 

k  r  j  e  r  ■  r 


V. 


we  obtain  the  total  H-field  inside  the  Miell 


Ik  (x  -  k 
o _ 

-J  R 


V 

n  = 


M  r. 


-2 1  cos  no 


k  S'  /  a.  \ 

-  H(2)  (k  a) 

2  n  o 


!  ( k  a )  +  ?■ ' 

n  o 


I  k 


°  V 


n»t  i 


n  . 


ko  a  •  , ,  ( 2 ) 

— - -  H  t <  a  '  .  I<  31  '  i. 

2  n  o  n  o 


As  the  frequency  •.»  -  0,  the  current  is  approximately  aivon 


.6  ) 


(12 


1 


J 


z 


I.  a 


a  ^  st 


1-  1  h 


where 


x  =  I /  u  e 
a  »  o  o 


R 

s 


a 
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SECTION  5 


DIFFUSION  COUPLING  MODELS  FOR  LIGHTNING  AND 
NUCLEAR  ELECTROMAGNETIC  PULSE 

In  this  section,  models  for  the  threats  of  nearby  and  direct  light¬ 
ning  strikes  as  well  as  nuclear  electromagnetic  pulse  (NEMP)  radiation  are 
postulated.  The  diffusion  coupling  formulas  derived  bv  Kaden  [1,2]  are  used 
to  obtain  internal  penetration  fields  for  homogeneous  shells.  Modeling  the 
excitation  time  waveform  as  a  double  exponential,  an  inverse  Laplace  trans¬ 
form  is  performed  as  in  [1]  to  obtain  the  interior  fields  as  a  function  of 
time.  These  fields  are  then  used  to  excite  transmission  lines  which  reside 
inside  the  shell  and  form  the  basis  of  the  computer-aided  design  (CAD)  pro¬ 
gram  discussed  in  Appendix  C. 

5.1  INTRODUCTION 

In  the  near  strike  case,  the  lightning  is  modeled  as  a  tube  of 
current  parallel  to  the  axis  of  the  shell.  It  produces  a  transverse  magnetic 
field  which  penetrates  the  shell  by  a  diffusive  coupling  mechanism.  The  time 
derivative  of  this  penetration  magnetic  flux  density  interacts  with  a  circuit 
loop  area  or  an  equivalent  transmission  line  area  to  produce  a  voltage  drop 
across  the  circuit.  At  sufficient  distances  from  the  lightning  current  column, 
the  electric  and  magnetic  fields  are  related  by  the  impedance  of  free  space. 

For  the  case  of  radiation  from  a  nuclear  electromagnetic  pulse  (NEMP), 
the  incident  field  is  taken  to  be  an  incident  plane  wave.  This  is  exactly  what 
the  near  strike  excitation  produces  far  from  its  source.  The  main  difference 
between  the  two  is  the  frequency  content  in  the  spectrum  of  the  incident  fields 
where  that  of  NEMP  is  much  higher.  Thus  the  same  diffusion  coupling  formulas 
apply  for  NEMP  fields  as  for  the  near  lightning  strike  case  as  long  as  the 
shell  cross  section  is  electrically  small  for  all  frequencies  of  significance. 

In  the  direct  strike  case,  the  lightning  current  waveform  is  assumed 
to  distribute  itself  uniformly  about  the  outside  of  the  shell  surface  and  is 
in  a  direction  parallel  to  the  axis.  It  gives  rise  to  an  electric  ;ield  on 
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which  are  nominal  values.  The  physical  length  of  the  bolt  will  be  assumed 
much  longer  than  the  aircraft  and  the  end  effects  of  the  clouds  are  neglected. 
The  bolt,  now  considered  as  an  infinite  tube  of  current  as  far  as  the  aircraft 
is  concerned,  gives  rise  to  $  directed  magnetic  field  lines  (because  of  sym¬ 
metry)  with  respect  to  an  axis  along  the  bolt.  At  a  point  R  in  the  absence 
of  the  aircraft,  the  magnetic  field  is  given  by 


Hext 


(t) 


Kt) 
2m  R 


(2) 


which  is  actually  nonuniform.  If  the  aircraft  shell  is  small,  however,  the 
external  field  is  usually  considered  to  be  uniform  over  the  shell  cross  sec¬ 
tion  and  R  is  some  mean  radius  from  the  lightning  bolt  to  the  shell. 

The  spectrum  of  Equation  (1),  and  hence  Equation  (2),  is  shown  in 
Fig.  5-2.  It  is  flat  out  to  approximately  2.7  kHz  where  it  rapidly  decreases. 
Thus,  a  low-frequency  analysis  (neglecting  the  term  jou  CqE  in  relation  to  aE) 
may  be  applied  to  the  shell  to  find  the  internal  magnetic  field  jl  C(t).  The 
result  is  given  irr  the  frequency  domain  by 


T  (s) 


HlnC  (s)  =  _ 1 _ 

ext  ,  .  cosh  z  +  %z  sinh  z 

H  (s) 


where 


Tint 


(s) 


spectrum  of  the  magnetic  field  inside  shell 


Hext 


(s) 


spectrum  of  the  magnetic  field  in  the  absence 
of  the  shell 


z 


t 


d 


,/tt7 

2 

diffusion  time  =  U  csd“ 


(3) 


=  1/d  x  volume  to  surface  ratio  for  nonmagnetic 

(y  =  'jq)  shells 
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NORMALIZED  AMPLITUDE  SPEC! RUM 


r 


Figure  5-2,  Normalized  Spectrum  of  Double  Exponential 
Lightning  Waveform 
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The  spectrum  of  the  incident  field  is  given  by 


wext  (s)  =  H  3  -  a 

o  (ct  +  s)  (8  +  s) 


where  =  Io/2ttR.  The  spectrum  of  the  internal  field  is 


Hlnt(s)  =  T(s)  HexC  (s) 


and  H  (t)  can  be  found  by  simply  taking  the  inverse  Laplace  transform  of 
Equation  (5).  Thus,  we  have 


Hint(t) 


.  f' 

27Tj  J 


_ (8  -  a)  eSt  ds _ 

(cosh  z  +  %z  sinh  z) (a  +  s)(g  +  s) 


where  K  is  an  arbitrary  constant.  It  is  simpler  in  Equation  (6)  to  use  the 
transformation 


z  *  i/  s  t 


to  obtain 


Uint(t)  (B  -  q) 

H  "  d  2-rrj 
o  J 


Z2  t/t 


2  2 

(cosh  z  +  sinh  z) (at^  +  z  ) (3t^  +  z  ) 


where  22  is  a  closed  contour  in  the  complex  plane  [1].  The  integrand  in  Equa¬ 
tion  (8)  has  simple  poles  at 

z  =  +  j  X  ( 9 ) 

— ■ J  n 


where  cot  X  = 

n  n 


±j 
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and 


z 


(11) 


The  poles  along  the  negative  Im(z)  axis  are  discarded  due  to  physical  reasons 
and  the  contour  of  integration  is  the  same  as  that  in  [1].  Thus  Equation  (8) 
may  be  computed  by  computing  the  residues: 


Hlnt(t)  3 

Ho 


2(6  -  a) 


+  +  R 

D  O 


+ 


oo 


I 

n=l 


(12) 


lim 


(z  -  i  \fa) 


2  . 
z  tf  t 


i  z-j  atd  (Cosh  z  +  £2  sinh  z)(z"  +  at^) (z~  +  Bt^) 


-at 


2td(6  -  a)  (cos  j/ottd  -  £  ^atd  sin  ^at^) 


R 


lim 


(z  -  J  )z  e 


2 

z  t/t . 


z  j/Btd  (cosh  z  +  r'z  sinh  z)(z"  +  at  )(z‘"  +  St  ) 

d  d 


-St 


2td(a  -  6) (cos  -  C  y/Std  sin 


(13) 


(14) 


2  . 
z  t/t 


lim 


(z  -  jA  )  z  e 
o 


z  ^Ao  (cosh  z  +  %z  sinh  z)(z“  +  at,)(z“  +  St.) 

d  d 


(15) 


I 
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Using  L'Hospital's  rule,  Che  limit  in  Equation  (15)  becomes 


-x;c/td 


(at.  -  xJ)(Bt.  -  xj)(l  +  5  +  (5A  )2)  sin  X 
dodo  o  o 


The  same  thing  is  done  for  R^.  Now,  usually,  t,  >>  1  for  the  shapes  to  be 
considered.  In  this  case 


X  =  nir  (n  ^  0) 
n 


This  allows  us  to  simplify  Rq  and  R^,  thus  giving: 


Hlnt(t) 


COS  y/atj  -  £  /°*d  Sin  /**d  C0S  /^d  ’  /^d  Sln  /^d 


~c/Ctd 

(3  -  a)  tdc,  e 

(1  -  td) (1  -  £0  td) 


2(3  -  a)t  » 

- — ? — i  2 


2  9 

-n  tt“  t/ 1 , 
,  nn  d 

(-1)  e  _ 


n=l  (n~TT  -  atd)(n"ir  -  Std) 
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For  voltages  induced  in  loops  or  transmission  lines  located  at 
points  inside  the  shell,  it  is  actually  the  time  derivative  of  the  internal 
magnetic  field  which  is  important.  This  is  given  by 


Hlnt(t) 


cos  /Btd  -  /std  sin  vAetd  cos  sin/^d 


-t/Ct 


(6  -  a)  e 

(1  -  c,a  td)  (1  -  £6  td) 


2(6  -  a) 


2  2  . 

-n  it  t/t 

/  i  \  ri  d 

(-1)  e 


*  ~ -v  (■  -  a#  -  s>) 


Alternatively,  one  may  obtain  Equation  (18)  by  convolving  the  im¬ 
pulse  response  h(t)  [1]  of  the  shell  with  the  excitation  of  Equation  (2). 
This  is  written  as 


Hint(t)  =  h(t)  *  HeXt(t) 


where 


h(t)  .  e  +  2  f  M,«  e-<n’r  t/cd 


5-8 


The  result  is  given  by 


which  is  equivalent  to  Equation  (18).  From  Equation  (22)  the  initial  values 
are  easily  seen  to  be 

HlnC(0)  =  0 


HLnt(0)  =  0 


(23) 


5.3  NUCLEAR  PULSE  EXCITATION 


The  situation  of  an  incident  plane  wave  radiated  by  a  NEMP  is  shown 
in  Fig.  5-3.  The  external  field  time  dependence  is  assumed  to  be  of  the 
form  [ 3 ] : 


..ext  „  ,  -at  -8t. 

H  =  H  (e  -  e  ) 
o 


where  typical  values  are  given  by 


154  amps /me ter 
6.3  x  10b  s_1 


3  =  1.89  x  108  s  1 


(24) 
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Aircraft  Shell  Cross  Section 


Figure  5-3.  NEMP  Excitation 

This  is  essentially  the  same  as  Equation  (1)  except  that  X  and  £  here  are  much 
larger.  This  accounts  for  the  substantial  increase  in  spectral  content  of  a 
NEMP  waveform  compared  to  the  lightning  case.  The  spectrum  of  Equation  (.24) 
is  shown  in  Fig.  5-4.  Although  the  frequency  content  in  H  here  is  much 
higher  for  the  near-strike  lightning  case,  the  formulas  in  Subsection  5.2 
can  still  be  used  if  the  overall  dimension  of  the  shell  is  small. 


NORMALIZED  AMPLITUDE  SPECTRUM 


5.4  DIRECT- STRIKE  LIGHTNING  ATTACHMENT 


We  assume  here,  for  simplicity,  that  the  attachment  current  dis¬ 
tributes  uniformly  around  the  shell  cross  section  shown  in  Fig.  5-5.  The 
actual  distribution  can  be  found  by  the  method  of  Subsection  2.8.  Thus,  the 
surface  current  density  Jg  is  approximated  by 

J  (t)  =  (25) 

s  C 

where  C  is  the  shell  outside  circumference  and  I(t)  is  given  by  Equation  (1). 
Coupling  to  the  inside  is  effected  through  the  surface  transfer  impedance  Z  . 
Thus,  the  internal  electric  field  is  in  the  same  direction  as  Js  and  is  given 
by 

in  the  frequency  domain. 


E.  (s)  =  J  (s)  Z  (s) 

int  s  st 


(26) 


Z  (s)  is  given  bv 
st 


Z 

st 


(s) 


n(s) 

sinh  y(s)d 


(27) 


where 


n(s) 

y(s) 

Now  we  let  t  ,  = 

a 


=  >v/s)ja 
2 

pad  and  the  inverse  transform  of  Equation  (26)  is  written  as 


E.  „(t) 
int 


I 

ad  C  2a j 


ds 


(28) 


where  K  is  an  arbitrarv  real  constant  and  Equation  (25)  has  been  used 


To  find  the  step  response,  lec 


t  >  0 


Kt)  = 


t  <  0 


Then  I(s)  =  I  /s.  Equation  (28)  then  becomes 
o 


E.  (t)  =  — 

int  crd 


_  f  2  e 

C  2t j  J  si 


z~  t/ 1 


inh  z 


dz 


(29) 


where  the  substitution  s  =  z~/t.  has  been  used  and  "  is  the  contour  shown  in 

d 

Fig.  5-6.  The  integral  in  Equation  (29)  has  simple  poles  at  zero  and  at 


2^  =  jn’T  ;  n  =  1,2,  .  .  . 

The  integral  in  Equation  (29)  is  thus  equal  to 


2rr  j 


y  R  -  I  R 


n=  1 


(30) 


where 


-  (ntt)  “  t/t 
2(-l)n  e  d 


Thus  the  step  response  is  given  by 


Estep 

int 


(t) 


1  4-  2 


V 

n=l 


(-1)‘ 


-  (n n ) 4 


e 


c/td 


(31) 
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Figure  5-6.  Contour  in  Complex  z  Plane  for  Computation 
of  Integral  in  Equation  (29) 


The  impulse  response,  obtained  by  differentiation,  is  given  by 


Elmp(t)  = 
int 


21  ^  „  o  - (n*rr> “  t/t 

2  (-1)°  (n-)2  e  d 

CTd  C  Cd  n=l 


Equation  (31)  is  the  same  as  that  abstracted  in  Apoendix  B  of  [4 
which  was  for  the  voltage  measured  between  longitudinally  spaced  points  on 
the  interior  surface  of  the  tube  due  to  a  current  step  function. 

Now  the  lightning  waveform  is  modeled  by 


I(t)  =  I  [e 
o 


-at  -fit. 


,  t  •>  0 


whose  spectrum  is  given  bv 


I(s)  =  I 


O  ( jL  +  s)  1 3  +  si 
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The  integral  in  Equation  (35)  is  given  by  the  sum  of  its  residues: 


where 


2TTj 


[  R  +  Rfl  +  f  R  1 
L  °  8  n-1  nJ 


-at 


2  sin 


/«d 


(8  -  a) 


-8t 


2  sin  y/gtj 


(a  -  0)  K8/td 


,  -  (nir)  t/t, 

-  (-l)n  (mr)2  e _ d_ 

(atd  -  (nir)2)(gtd  -  (nn)2) 


Thus  Equation  (35)  can  be  written  as 


Eint(t)  " 


/atd  e 


-at 


~0t 


sin  y a'td  sin  J 6td 


2  (0  -  a)  t,  2 


.  -(niT)  t/t, 
(-l)n  (nu)2  e _ 

d  n-1  (atd  -  (mr)2)(0td  -  (rnr)2) 


(36) 


From  [5],  we  have  the  sum: 


2  cos  nx 

i  2  2 

n-l  n  -  a 


2a 


n  cos  (x  -  tt)  a 
2a  sin  ita 


so  that 


dS  V 

dx  =  ^ 


n  sin  nx 

i  2  2 

n-1  n  -  a 


air  sin  (x  -  it)  a 
2a  sin  ua 
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and 


2 

n  cos  nx 


a7T  cos  (x  -  tt)  a 
2  sin  ira 


Now  let  x  *  tt  and  we  have 


or 


aft 

2  sin  Tra 


00 


-  2 

n-1 


2 

n 

2 


n 


__u _ 

sin  u 


Thus  we  also  have 


oo  n  ? 

2  2  (-D  (nff) 

n»l  (mr)2  -  u2 


sin 


tt 
tt 


2  2 


(mr)2  (~l)n 
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sin 


2  2  (n1T)‘ 
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a 


Now  let  *  ctt^  -  (nTT )  and  x^g  *  8t^  -  (nTt)  then  Equation  (25)  becomes 
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This  can  also  be  derived  by  convolution  and  is  identical  to  the 
voltage  response  derived  in  [4].  Thus  the  voltage  drop  per  unit  length 
inside  the  tube  wall  along  the  tube  axis  is  given  by 


Eint(t> 


-at 


-3t 


sin 


V*  aya 


sin  y  Bya  d 


(38) 


+  2 (a  -  0)  t,  2 


-(nit)  t/t, 
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SECTION  6 


INTERIOR  FIELD  RESULTS  AND  APPLICATION  OF  LOW 
FREQUENCY  COUPLING  MODELS 

6.1  RESULTS  AND  APPLICATION  OF  MODELS 

Computer  programs  were  written  for  the  E-field  and  combined-field 
formulations  presented  in  Section  2  of  this  report.  These  are  documented  in 
Appendix  A  where  sample  output  is  given  for  a  cylinder  of  circular  cross  sec¬ 
tion  with  a  radius  of  0.3828m.  The  frequency  of  the  incident  plane  wave  is 
300  MHz.  The  TE  and  TM  surface  currents  of  the  two  formulations  are  In  excel¬ 
lent  agreement  at  all  points  on  the  cylinder  contour.  An  additional  check  is 
provided  by  the  exact  series  solution  where  a  computer  program  along  with 
sample  output  is  also  documented  in  Appendix  A.  Electric  surface  currents 
and  scattered  far-field  patterns  for  all  three  methods  are  in  excellent  agree¬ 
ment.  Several  examples  of  surface  current  density  and  scattered  field  patterns 
may  be  found  in  [1],  [2],  and  [3],  so  no  plots  of  these  quantities  are  pre¬ 
sented  here.  The  main  purpose  of  Section  2  is  actually  to  provide  the  various 
E-field  impedance  operators  necessary  for  the  shell  penetration  formulations 
of  Section  3. 


An  excitation  other  than  an  incident  plane  wave  is  considered  at  the 
end  of  Section  2  in  which  a  longitudinal-directed  current  is  assumed  to  exist 
on  the  contour  C.  This  is  used  as  a  first-order  simulation  of  a  direct-strike 
lightning  (DSL)  excitation.  If  the  contour  C  is  not  circular,  the  surface  cur¬ 
rent  density  is  not  constant  and  distributes  itself  around  C  inductively.  A 
computer  program  is  documented  in  Appendix  A  which  is  used  to  compute  this 
distribution.  An  example  of  a  computation  for  this  excitation  is  shown  in 
Fig.  6-1  where  the  cross  section  of  a  fuselage  station  [4]  is  considered.  The 
numbers  near  the  dots  on  the  contour  identify  the  center  of  each  A j ,  j=l,  2, 
...,  28  (see  Fig.  2-2)  and  the  numbers  with  arrows  pointing  to  the  contour 
Indicate  surface  currents  computed  at  that  point.  Numbers  in  parentheses  indi¬ 
cate  currents  computed  in  [4]  by  solving  Laplace's  equation  in  a  finite  region 
surrounding  C  with  the  effects  of  three  symmetrically  placed  return  conductors 
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included.  The  integral  equation  solution  of  Section  2  in  effect  has  the  return 
path  at  infinity.  This  accounts  for  the  difference  in  results  obtained  by  the 
two  methods.  Note  that  the  current  distribution  is  symmetric  about  a  vertical 
centerline  of  the  fuselage  section  so  that  only  half  of  the  data  is  shown  in 
the  figure. 

Computer  programs  using  the  impedance  sheet  approximation  of  Sec¬ 
tion  3  and  the  exact-series  solution  of  Section  4  for  a  lossy  shell  (i.e.,  a 
shell  wall  material  having  finite  conductivity  a)  are  documented  in  Appendix  B. 
These-  programs  compute  an  internal  z-directed  electric  field  for  the  TM  case 
and  magnetic  field  for  the  TE  case.  Some  sample  plots  of  shielding  effective¬ 
ness  (defined  by  Equation  (1)  of  Section  3)  are  shown  in  Figs.  6-2  and  6-3 
for  a  shell  of  circular  cross  section.  As  expected,  the  impedance  sheet 
approximation  gives  excellent  results  in  the  TM  case  when  compared  to  the 
exact-series  solution.  This  comparison  is  shown  in  Fig.  6-2  for  different 
conductivities  and  over  a  large  frequency  interval.  Here  the  exact-series 
solution  is  represented  by  the  solid  lines  and  the  impedance-sheet  integral 
equation  solution  is  represented  by  circles.  The  latter  becomes  suspect  at 

g 

frequencies  much  higher  than  10  Hz  because  the  number  of  subsections  with 
which  the  circular  contour  was  approximated  was  held  constant  (NC  *  16  for 
all  frequencies).  The  TE  comparison  shown  in  Fig.  6-3  is  not  as  good  because 
for  this  case  there  are  transverse-directed  components  of  polarization  current 
which  are  neglected  in  the  impedance  sheet  formulation. 

The  results  for  the  TM  case  indicate  that  the  z-directed  incident 
electric  field  is  effectively  not  shielded  at  all  as  the  frequency  becomes 
low.  This  can  be  seen  from  the  impedance  sheet  integral  equation  where  the 
vector  potential  term  becomes  negligible  as  ur*0  and  hence  J  -*•  E^/Z^.  The  TM 
case  is  somewhat  pathological  at  low  frequencies  because  theoretically  no 
charge  separation  can  occur  to  create  a  scattered  electric  field  which  tends  to 
produce  zero  total  electric  field  inside  the  shell  as  occurs  in  the  TE  case.  A 
physically  realizable  shell  will  always  have  a  quasi-static  charge  separation 
at  low  frequencies  and  thus  will  tend  to  shield  an  incident  electric  field 
(more  than  the  TM  infinite  cylinder  result  would  indicate).  As  expected,  the 
magnetic  field  is  essentially  not  shielded  at  low  frequencies  regardless  of 
the  polarization  considered. 


6-2 


Results  for  the  traveling  wave  approximation  to  shell  penetration 
developed  in  Subsection  3.3  are  not  presented  here.  A  computer  program  was 
written  which  incorporated  this  formulation  but  was  found  to  give  questionable 
results  at  low  frequencies.  This  is  because  the  relationship  between  the  tan¬ 
gential  electric  and  magnetic  fields  just  inside  the  shell  is  not  given  by  the 
characteristic  impedance  of  the  material  inside  the  shell.  This  assumption 
essentially  neglects  the  inductive  reactance  caused  by  the  shell  geometry  which 
is  seen  by  the  electric  current  used  to  excite  the  equivalent  transmission  line 
model.  The  traveling  wave  approximation,  however,  does  give  reasonable  results 
for  higher  frequencies  at  which  the  shell  cross  section  dimension  is  electri¬ 
cally  large. 

The  various  computational  methods  presented  so  far  allow  one  to  ade¬ 
quately  determine  the  interior  fields  which  penetrate  a  shell  at  particular 
frequencies  of  an  external  steady-state  electromagnetic  field.  Thus  if  the 
spectrum  of  the  excitation  is  known,  one  can  characterize  the  spectrum  of  the 
interior  field.  In  Section  5,  some  expressions  are  derived  for  the  interior 
field  spectrum,  given  an  assumed  external  field  or  excitation.  If  a  low  fre¬ 
quency  assumption  can  be  made  concerning  the  excitation  spectrum,  then  the 
interior  field  can  be  obtained  as  a  function  of  time.  It  is  also  assumed  to 
be  uniform  in  space  over  a  region  which  encompasses  the  shell  cross  section. 

The  computation  of  these  interior  fields  which  are  due  to  a  nuclear  electro¬ 
magnetic  pulse  (NEMP) ,  near-strike  lightning  (NSL) ,  or  direct-strike  light¬ 
ning  (DSL)  as  they  are  modeled  in  Section  5  are  performed  by  the  computer  pro¬ 
gram  subroutines  which  are  briefly  described  in  Appendix  C. 

Figure  6-4  shows  plots  of  internal  magnetic  fields  for  various  lossy 
cylinders  when  exposed  to  an  NSL  excitation.  The  lightning  current  is  100m 
from  the  shell  and  is  parallel  to  the  shell  axis.  As  expected,  as  0,  d,  or 
the  volume- to-surf ace  ratio  is  increased,  the  rise  time  of  the  responding 
interior  field  increases.  This  can  be  seen  from  the  equivalent  circuit  ana¬ 
logues  to  the  coupling  mechanism.  Inside  the  shell,  a  circuit  or  transmission 
line  .3  exposed  to  the  interior  field.  A  worst-case  situation  is  assumed  where 
the  effective  area  of  the  transmission  line  or  circuit  is  maximally  coupled  to 
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the  transverse  component  of  magnetic  flux  density.  Time  derivatives  of  in¬ 
terior  and  exterior  NSL  fields  are  shown  in  Fig.  6-5. 

Using  the  transmission  line  excitation  formulas  developed  in  [5], 
the  open-circuit  voltage  and  short-circuit  current  measurements  made  for  a 
given  loaded  transmission  line  configuration  may  be  computed.  An  example  of 
these  results  is  shown  in  Fig.  6-6  where  the  interior  circuit  consists  of  10m 
of  RG8-A  along  the  cylinder  shell  axis  terminated  in  a  20-ohm  load.  The  product 
of  the  above  quantities  as  a  function  of  time  gives  an  upper  bound  on  the  instan¬ 
taneous  transient  power  possible  across  the  terminals  of  the  load.  An  example 
of  this  is  shown  in  Fig.  6-7.  Various  parametric  curves  are  of  interest  from 
a  computer-aided  design  (CAD)  standpoint  so  that  elementary  circuits  may  be 
constructed  which  minimize  the  possibility  of  component  burnout.  Examples 
of  these  are  shown  in  Figs.  6-8  through  6-11.  This  type  of  parametric  repre¬ 
sentation  is  discussed  in  more  detail  in  [5]. 

As  mentioned  previously,  the  main  difference  between  NSL  and  NEMP 
excitation  is  the  spectral  content  of  the  two,  that  of  NEMP  being  much  higher. 
Some  representative  interior  field  coupling  examples  for  NEMP  are  shown  in 
Figs.  6-12  through  6-14.  Note  that  the  time  scales  are  now  different  and  the 
high-frequency  resonance  on  the  transmission  line  is  visible. 
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3740  (3763) 


Figure  6-1.  Free  Space  Longitudinal  Current  Distribution  on  Conducting 

Cylinder  Contour,  1  ,  =  20  kA 

7  ’  total 


Figure  6-2.  Electric  Shielding  Effectiveness  (TM  Case)  at  Center 
Circular  Shell,  Radius  =  0.5m,  Thickness  =  1mm 


EXACT  SERIES  SOLUTION 


Figure  6-3.  Magnetic  Shielding  Effectiveness  (TE  Case)  At  Center  of 
Lossy  Circular  Shell,  Radius  =  0.5m,  Thickness  =  1mm 
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Figure  6-4.  External  and  Internal  NSL  Magnetic  Field  for 
Circular  Shell 


Figure  6-6.  Open-Circuit  Voltage  and  Short-Circuit  Current  Induced  on 
Transmission  Line  Inside  Shell  of  Figure  6-5 


Upper  Bound  on  Induced  Power  Available  in  Transmission  Line 
of  Figure  6-6 


Figure  6-8.  Log  of  V  and  I  Versus  Transmission  Line  Length  for 
Shell  A  of  Figure  6-4 
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Figure  6-9.  Los  of  P  and  Wunsch  Constant  Versus  Transmissi 
&  "  max 

Line  Length  for  Shell  A  of  Figure  6-4 


Figure  6-10.  Log  of  and  I  Versus  Normalized  Load  Impedance 
for  Shell  A  of  Figure  6-4 
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Figure  6-12.  External  and  Internal  NEMP  Magnetic  Field  for 
Circular  Shell 
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Figure  6-13.  Open-Circuit  Voltage  and  Short-Circuit  Current  Induced 
on  Transmission  Lines  Inside  Circular  Shells  of 
Figure  6-L2 
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TIME  (us) 

(a)  Radius  •  la.  Thickness  ■  8  Ply,  3  •  10"* 


TIME  (us) 

(b)  Radius  *  10a,  Thickness  ■  8  Fly,  3  «  2  x  10"* 

Figure  6-14.  Upper  Bounds  on  Induced  Power  Available  in 
Transmission  Line  Cases  of  Figure  6-13 
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APPENDIX  A 


PROGRAMS  FOR  PERFECTLY  CONDUCTING  CYLINDERS  OF 
ARBITRARY  CROSS  SECTION 

The  purpose  of  this  appendix  is  to  define  the  necessary  data  cards 
required  by  the  E-field  and  combined- field  programs.  The  main  program  seg¬ 
ments  are  listed  in  Subsections  A.l  and  A. 2.  A  program  to  calculate  the 
exact  series  solution  is  included  in  Subsection  A. 3.  The  function  subprograms 
and  matrix  element  subroutines  are  not  explained  in  detail. 

The  first  thing  one  must  do  is  approximate  the  contour  C  of  the 
cylinder  by  a  finite  number  of  straight  line  segments.  Best  results  are 
usually  achieved  when  A  <_  0.1X  which  puts  a  limit  on  the  electrical  size  of 
the  objects  considered  since  matrix  methods  are  being  used.  An  example  of 
approximating  C  is  given  in  Fig.  2-2.  Note  that  the  contour  need  not  be 
closed  for  the  E-field  formulation  but  must  be  closed  for  the  combined- field 
formulation.  The  excitation  is  such  that  the  incident  magnetic  field,  H1,  is 
equal  to  unity  at  the  origin.  This  is  done  for  both  polarizations. 

Data  is  read  from  data  cards  in  the  main  program  according  to  the 
format  statements: 

100  FORMAT (615) 

101  FORMAT (2E20. 7) 

102  FORMAT (6E11. A) 

The  data  cards  appear  In  the  sequence  shown  in  Table  A-l  and  are  defined  as 


follows : 

NGQ 

=  Order  of  Gaussian  quadrature  formula 
used  to  approximate  integrals 

A(i),  T(i)  for 
i  -  1,2, . . . ,  NGQ 

=  Weights  and  nodes,  respectively,  of 
Gaussian  quadrature  formula 
(divided  by  2) 

ITM  »  integer  option 

*  1  for  TM  case 

=  0  bypass  TM  case 

A- 1 


ITE  =  integer 
option 

ISC  =  integer 
option 

NX 

NP 

PHIO 

DPHI 

PHII  (i) 

AMU 

EPS 

BETA 

(xi,yi)  for 

i  -  1,2 . NC+1 


1  for  TE  case 
0  bypass  TE  case 

1  for  normalized  scattered  field 
pattern  to  be  computed 

Number  of  angles  at  which  plane 
wave  is  incident 

Number  of  points  at  which  scattered 
far  field  is  computed 

First  angle  at  which  scattered  far 
field  is  computed 

Increment,  in  degrees,  at  which  far 
field  pattern  is  computed 

i  =  1,2,..., NX  =  Angle  of  plane 
wave  incidence  measured  in  degrees 
counter-clockwise  from  x-axis 

Permeability  of  material  in  which 
cylinder  is  imbedded 

Permitivity  of  material  in  which 
cylinder  is  imbedded 

Used  only  for  combined-field 
program 


x,y  coordinates  of  end  points  of 
straight  line  segments 


NFR 

FMC 


Number  of  frequencies  of  incident 
plane  wave  to  be  considered 

Frequency  in  MHz  (read  NFR  times) 


The  three  programs  were  run  for  a  cylinder  of  circular  cross  section  with  a 
radius  of  0.3828m  at  300  MHz.  The  results  follow  each  program  listing. 
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Table  A-l.  Input  Data  Card  Sequence 


Format  Number 

Data  Punched  on  Card 

100 

NGQ 

101 

A(l) ,  T (1) 

l 

l  l 

101 

A(NGO) ,  T (NGQ) 

100 

ITM,  ITE,  ISC,  NX,  NP 

102 

PHIO,  DPHI 

102 

PHII  (1) 

102 

PHI  (NX) 

102 

AMU,  EPS,  BETA  (only  for 
combined-field  program) 

100 

NC 

102 

X1  yl 

102 

x2  y2 

# 

,  , 

• 

•  • 

• 

*  * 

102 

^C+l  yNC+l 

100 

NFR 

102 


FMC  (NFR  times) 


The  various  subroutines  and  function  subprograms  needed  by  each 
program  will  not  be  described  in  detail.  Instead,  they  are  shown  in  Table  A- 2 
corresponding  to  the  computation  for  which  they  are  used.  The  variables  which 
are  stored  in  common  blocks  are  defined  in  Table  A-3. 


Table  A- 2. 

Subroutines  Corresponding 

to  Computation 

Equation  Number 
in  Section  2 

Computation 

Subroutine 

pp  2-5  to  2-7 

x-y  components  of 

t  ,  y  R 

-in  m  — m 

CDATA  for  E-field 

CDATACF  for  combined- 

field 

(22) 

elements  of  [Z  ] 

ZMNE 

(26) 

elements  of  [Z^1] 

SZH 

(27) 

-*-ie 

elements  of  V 

TMX 

(28) 

elements  of  V^'n 

TEX 

(33) 

0 

elements  of  [T  ] 

TMNE 

(36), (37) 

elements  of  [T^1] 

STH 

(46) 

6Vie  +  Iie 

THXCF  (for 
combined- fie Id) 

(47) 

6Vih  +  Iih 

TEXCF  (for 
combined- field) 

IB) 

VcA 

TMS  for  TM  case 

TES  for  TE  case 

— 

»o<2>  <*> 

HANK02 (X) 

— 

H1(2)  (x) 

HANK12 (X) 

(23) 

a(z) 

ALPHA  (z) 

— 

- ►  -f 

Solve  Ax  =  b 

DECOMP  and  SOLVE 

for  x 

Solve  Ax  =  b 

GELS 

— ► 

for  x  when 

A  is  symmetric 

J  ,  Y 
n  n 

BES  for  exact 
series  solution 
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Table  A-3.  Common  Block  Variables 


Block 

Variable 

Meaning 

GQI 

NGQ 

Order  of  G-Q  integration 

A 

Weights/2 

T 

Nodes/2 

CUV 

ULX 

t  •  X 

— m  — 

ULY 

t  •  y 
-m 

NC 

NC 

c 

RCX 

R  •  x 
— m  — 

RCY 

R  •  y 
— m 

DC 

A 

m 

CK 

RKX 

k  •  RCX 

RKY 

k  •  RCY 

DK 


k  •  DC 


A. 1  E- FIELD  PROGRAM 


1 .  . . main  RR»GoaM  s-*?1?  C OM*'j T T M3  INDUCE*'  ELECTPrr  c,.lR9F,','rc 

?.  C  ON  T*ft-DI'-ENSI0NAL  CONDUCTING  SMAPrS, 

3*  C  USES  E-FIrUD  INTEGRAL  EG'J  AT  I  ON 

A.  DIMENSION  PHjniCl 

5.  COMPLEX  Z! 1000) * VM<500) 

6.  COMPLEX  ZMNE,S2H 

7.  COMMON  /CUV/,JLX(60)*,JLYC60)*NC/C/RCX(O>*°CY(aC>  .0C<fD) 

9.  COMMON  /3Q I / a ( 1 0 ) #  T  ( 1 0 )  * NGO 

9*  COMMON  /CK/P*X(60)>R<Y(&0)*J<<6D) 

IQ.  DATA  P I /3 • 1 4i 593/ 

11.  C--------R£AO  IN  AvQ  aR  I  NT  OUT  INPUT  DATA 

12.  100  P9RM«T(6I5) 

13.  101  r0RMAT I 2E2C  *  7 ) 

14.  102  FORMAT! 6E1 1 •* ) 

15*  READ (105*100)  NGQ 

16.  READ(io5#ici)  (A{  I  )#T  {  I  )  *  I  »i  * KIG0 ) 

17.  OUTPUT * NGQ 

19*  READ! 105*100)  I Tm* I Tr * I So NX . No 

19.  OUTPUT* ITM, ItE* ISCiNXiNP 

20*  READ! 105* 102)  PMIo/OPMI 

21.  OUTP'JT*PHIOjOPPI 

22.  DO  A  I  *  1  *  NX 

23.  READ ( 105* 102 )  °HII(I> 

24.  OUTPUT*I*PHlt( I) 

25.  4  CONTINUE 

26.  REAQ(  105*  102)  AM'J.EPS 

27.  OUTPUT* AMU* E»S 

29.  CALL  CDA T A  {  I-'C) 

29.  READ! 105. 100)  NFr 

30*  00  50  INF«i*vFR 

31«  READ(  105*  102)  rv,C 

32.  aRITE! 109* 30p)  F«C 

33.  AK»2.*PI*FVC#S0RT( AMu*EPS)*l.E  6 

34.  DO  5  I“l*Nc 

35.  RKX  (  I  )  «RC,X  (  I  )  *A< 

36.  RKY! I ) »RCY (  I  ) *A< 

37.  D<( I )*0C! I  )*4< 

39.  5  CONTINUE 

39*  IF( ITM.NF. i )  Q9  TO  30 

40 •  C- - ----FORM  UPPER  RT,  TRIANUUE  9r  v  A  TR I  X  , 

4f  <«1 

42.  DO  1  I N« 1  *  NC 

43.  DO  1  I M ■ 1  *  I  N 

44.  Z!<)»ZMNE! T“. ’ N ) 

45.  <*<T1 

46.  1  CONTINUE 

47 .  C- - --FORM  TM  E yC I T AT  I  0m  VECTORS. 

» 5 •  CALL  TMX! vm*dwJ I/NX) 
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*9  • 

c »  •  - 

..---Sc*LVE  F»R  NORMALIZES  ELECTRIC  "  jPPrNTS. 

5C  • 

( Nc*l ) /2 

5  1  • 

CALL  3ELS<V",Z*NC*NX,MR1 

52* 

<*1 

o3 . 

09  2  1»1#NX 

54  • 

wRlTE (103*200)  3H!!(]| 

55  ♦ 

*R I  TP ( 108* 20i  ) 

56  • 

09  3  J*1»NC 

37* 

VMM«CA8S( V«( <) ) 

58. 

«9ITP(  103*205)  JfV'KO  ,Vmm 

59  • 

<»K*1 

60  • 

3 

CONTINUE 

61  • 

2 

C9NTINUE 

6?. 

IF f  ISC.EQ.l  )  CALL  T^S(VN,NX.'awn*JS'-II^,0S|4I*N^) 

63  • 

200 

r9RMAT( <0'*lRX# 'T*  CURRENTS  P9R  °MI  « • , E 1 1 « 4 # sX . ' 0r3RFcS '  ) 

64, 

201 

P9RMAT  (  '  1 , 6Xj  '  °ULSE  NO  ,  >,  1 1 X  ,  •  R£  AL  •  ,  1 1  X  ,  »  I  'I  A  3  1  ,  1  j,  x  *  •  *  *3 . 1  ) 

65* 

202 

F  9  R  M  A  T  {  1  O Ii 5*5X*3El5.6) 

66  • 

3G 

I F  c ITE.NE.l )  39  T9  60 

67  • 

C  — 

--•••FORM  UPPER  RT.  T2IANULE  9p  TE  Z  matrix. 

68  • 

<«1 

69  • 

I B« 1  + 1  DC 

70 

09  31  IN.IB*MC 

71  • 

09  31  IM»I8*TN 

72  • 

Z(<)«SZH(  In-i  /  IN-1*  1*  1  )+SZUt  IM-t  *  IN*i  «-1  )•► 

73. 

1SZH(  IM»  IN.o.l,l)+SZH(  1-1,  IN, -1/-1 ) 

74. 

<»K  +  1 

75  • 

91 

C9NTINUE 

76* 

C--- 

- -FORM  TE  EXCITATION  VECTORS . 

77* 

CALL  TEX(VOOHH*NX,  I  DC  J 

78. 

C*  -  - 

. S9LVE  FOR  NORMALIZED  T£  CURRENTS. 

79. 

NF-NC-I3+1 

■iO 

mb*,\ie*(np  +  i  )  /5 

81 . 

CALL  GELS( VM, Z#nf,nX#mr) 

72* 

<*  1 

83* 

09  32  I  *  1  *  N X 

34. 

AR1TE<  108*300)  BWIH!) 

35* 

WRITE! 108*301  ) 

36. 

09  33  J* 13*  Nr 

87* 

VMM«CABS( VM</)  ) 

38. 

wRI TE( 103*505)  J*vM(K)*vmM 

39  • 

<■<  +  1 

90  « 

33 

C9NT INUE 

91. 

32 

continue 

9  P. 

IF(ISC.EO.I)  CALL  TES(Vri*NX*pMI  I*dHIm»0ow!*N5,  IOC) 

93. 

6C 

CONTINUE 

94  . 

^0 

CONTINUE 

95. 

300 

F8RMAT(  *0'*1RX#  'TE  CURRENTS  p9R  PMI  a  <  ,  E 1 1  •  4 ,  -?X ,  Or'iRrcS  <> 

96. 

301 

FORMAT {»  * #3y# « TRIANqle  N9 . » » 1 2X * • RE A L • # 1 1 X #  » T M A 0 • , 1 1 < , • M A 3 . 

97* 

3C2 

FORMAT (  *  1  ’  * 'FREQUENCY  9H  PLANE  vavf  a • *  E 1 5 • 7 # 9X ,  • MuZ  '  ) 

98* 

STOP 

99* 

END 
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Sample  output  data  is: 


PUN 

NGQ  «  4 
ITM  «  1 
ITE  *  l 
ISC  »  1 
NX  *  1 
N  P  ■  4 

PHIO  ■  .000000 
UPHI  •  90.0000 
I  »  1 

PHIMI)  ■  .000000 
AMU  »  1 .25SOO0E-06 
EPS  ■  8.350000E-1? 

-no,  or  STRAIGHT  t.INE  SEGMENTS  aoprRX I « AT l NQ  C  •  24 


JLX 

ULY 

PCX 

RCY 

oc 

1301*: 

CO 

-•9915E 

00 

•3763E 

00 

-  .4955E* 

01 

•9995E-01 

3833E 

GO 

-•9236E 

00 

•  35Q6E 

00 

-•1452E 

00 

•9993E-01 

6o’35E 

00 

-  *  7936E 

00 

•30UE 

00 

- »23loE 

00 

• 9993E-01 

7936E 

00 

- «60*5B 

00 

•231oE 

00 

-•3011E 

00 

•99932-01 

9236E 

00 

- •  3833E 

00 

» 1>52£ 

00 

•»3506E 

00 

•9993P-01 

9915E 

00 

-•1301E 

00 

•4955£ 

•01 

•»3763E 

00 

• 9995F-C1 

991 5E 

00 

•  1 301 F 

00 

-.4955E 

•01 

• • 3763E 

00 

»  9995F -01 

9P36E 

00 

*  3833E 

00 

-•1452E 

00 

• • 35Q6E 

00 

*  9993E-01 

7936E 

00 

•  f>0^5E 

00 

••231oE 

00 

■ • 301  IE 

00 

^OSSE-C* 

6035E 

00 

. 7936E 

00 

- « 301 1 E 

00 

-.2310E 

00 

• 9993E-nl 

3833E 

00 

• 9236E 

00 

- *3506E 

00 

• » 1 452E 

00 

, qo93r-0l 

1301E 

00 

.9915E 

00 

• « 3763E 

00 

•  *  4955E* 

01 

•a995E-01 

1301E 

00 

•991 5E 

00 

-•3763E 

00 

.4955E- 

01 

• a995E-n1 

3S33E 

00 

♦9236E 

00 

-  *  3506E 

00 

•1452E 

00 

•9993E-01 

6035E 

00 

*  7936E 

00 

-•3011E 

00 

•2310E 

00 

. 9993E" 0  * 

7936E 

oo 

»60*PE 

00 

-»23l0t 

00 

•301  IE 

00 

•99932-01 

9236E 

oc 

» 3333E 

00 

•  *  1 452E 

00 

» 3506E 

00 

•99932-01 

9915E 

00 

•1301E 

00 

-«4955E 

■01 

•3763E 

00 

•99952-01 

9915E 

00 

• • 1 301 E 

00 

•4955E 

•  01 

•3763E 

00 

•99952-01 

9236E 

00 

*  *  38332 

00 

•14SPE 

00 

« 35n6E 

00 

•99932-01 

7936E 

00 

•  • 6Ca5E 

00 

♦231cE 

00 

•  "3  0 11 E 

00 

,9993r.«< 

6o35E 

00 

- »  79962 

00 

•3011E 

00 

•  231 qE 

00 

•9993E-01 

3S33E 

00 

-•92362 

00 

. 3506E 

00 

•1*52E 

00 

• 9993E-01 

1301E 

00 

• *991 SE 

00 

»3763E 

00 

.4955E- 

01 

. 39Q5E-0^ 

•CLOSES  CONTOUR,  IOC  «  0 


FREO'JENCV  92  pl^E  WAVE  • 


330C000E  S3 


TM 

CURRENTS  F39  B*I  ■ 

.OOOOE  00 

OFjorrq 

». 

REAL 

I-A3 

*  AO. 

1 

-*1252552  01 

•1678932  01 

.3095492 

01 

2 

-.8825872  00 

•  1 73334E  M 

.  1^*9792 

01 

3 

••188352E  00 

•1781662  01 

•  1791592 

01 

4 

«6042»7E  00 

•1403262  01 

.1527842 

01 

s 

» 1073952  01 

•608l 34T  on 

,1 333662 

01 

6 

•9164922  00 

••2252172  00 

.9437012 

00 

7 

•3237A7E  00 

-.3937962  00 

•6,07312 

00 

8 

•«1734=;7E  00 

-•4262602  00 

.4621092 

00 

9 

■•29C824E  00 

- • 76156*E-31 

.3006402 

00 

10 

•♦1427502  00 

•  1021 49E  00 

.1755332 

00 

1 1 

•lb7C?9E-Cl 

•  7644182-01 

. 7a03802 • 

01 

12 

•  8602Q8E  »0 1 

•2012912-02 

.8605332- 

01 

13 

♦86q3i 6E*0l 

•2015032-02 

.8605522. 

01 

1* 

•157065E-01 

• 764405E-31 

.7303742- 

01 

15 

•♦142748E  00 

•1021472  00 

.1755312 

00 

16 

•  ♦  29q8^6E  00 

-.7615442-01 

.3306412 

01 

17 

• • 1 784^7E  00 

- *4263622  oo 

•  442 1 1 02 

00 

18 

•3237<aiE  00 

-.5987962  00 

.6a07322 

oo 

19 

•916436E  00 

••2252172  00 

.9436942 

00 

30 

•1073372  01 

•6081352  30 

.1233672 

01 

21 

.6042832  00 

.1403252  01 

.1527832 

01 

22 

••  18?351E  00 

•1781682  01 

.1791612 

01 

23 

-•3825882  00 

•1783332  01 

.  la89752 

01 

24 

•  •  1 25255E  01 

.1679942  01 

.209549F 

01 
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SCATTERED  FIELD  rattcrm  rgg  oujl 


-‘O-'tr  30 


T  R  I ANGLE 


*ST9°* 


dm  I 

• COOOE  00 
.9000E  0? 
. i«00E  03 
•  2700E  03 


SOBTfS/L) 
« 1 1322  01 
.10892  01 
•2600E  01 
•  1 089E  01 


TE  CURRENTS  2:5:3  »*!  •  .QO00r  03  DEGREES 


v«. 

PEAL 

I  M  A  G 

«A3. 

1 

• • 154673E 

01 

.  101632E 

01 

.1851032 

» i 

2 

145044E 

01 

•  1 1 40^3E 

01 

.1H45Q9E 

01 

i 

-•U23S9E 

01 

. 144623E 

01 

•1331402 

01 

4 

-•5087«2E 

00 

. 172797E 

01 

.1801312 

01 

5 

•3198S2E 

00 

• 167835E 

71 

.  170S56E 

01 

6 

•105135E 

01 

.  1036132 

01 

.151162=! 

01 

7 

• 1235S7E 

01 

.  102306E 

00 

.1259632 

01 

8 

•9002s9E 

00 

/65437E 

00 

. 1131632 

01 

9 

• 22o4?5E 

00 

• • 105381 E 

01 

.1076622 

01 

10 

-•27S2R5E 

00 

-./13143E 

00 

.7655202 

00 

11 

• . 3873p5E 

00 

- . 694841 E 

•01 

.3940002 

00 

12 

• • 2&o 1®7E 

00 

•*33124E 

00 

.5437362 

00 

13 

•. 1743a3E 

00 

.6^34572 

00 

.715047=- 

00 

14 

• • 26q2o1E 

00 

• *331252 

00 

.5437392 

00 

15 

-♦3873r7E 

00 

-.694811 E 

•01 

.3943022 

00 

16 

•*2782«8E 

00 

••7131355 

00 

.7655132 

00 

17 

•  22o**?7E 

00 

-.105380E 

01 

.1076612 

01 

13 

*9ooHs8E 

00 

•  •  7  65431 E 

00 

.1151672 

01 

19 

•1285S7E 

01 

.102808E 

00 

•  1289682 

01 

20 

•105135E 

01 

.108612E 

01 

.1511622 

01 

21 

•3193S4E 

00 

.1678362 

01 

.  1  70»572 

01 

22 

• . 50  8  7a5E 

00 

.1727982 

01 

.  18Q1332 

01 

23 

••1123S9E 

01 

.1446252 

01 

.1831422 

01 

24 

-.145045E 

01 

•  1140*42 

01 

. 134510E 

01 

T2  SCATTERED  F I  ELD  PATTERN  2»g  PwT!  *  *000^2  00 


OH  I 

< COOOE  00 
.9000E  0? 
.i«OOE  03 
.2700E  03 


SORT ( S/L ) 
.1026E  01 
•661AE  00 
•  1476E  -'l 
•6614E  00 
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COMBINED- FI ELD  PROGRAM 


A.  2 

1.  sriri.?  -  ■jowtjL  AT  T^'J  F»R  09*o'jTIN3  3i_Avr  *Avr 


3.  r  3CATTc«imt  *n  3ERFECTLv  C9\!0'JCTTV3  CYLINOrRS  »r 

3.  C  ARB  I  TP  A9  Y  S^AB(T, 

3  I  MENS  1 9m  OHf  I  (  io)i  I»S(bn) 

5.  COMPLEX  7  ( 1  J?0  )  ,  VM(500)  *0UM(5D0) *XC(  SO)  *  C’JR  < 5~ ) 

6.  ^?MPLEX  T^'E,Z«N<EiSTwJSZta| 

7.  C9MM9\  /CUV/'!LX(  60)  *  JLYc  60)  .'  C/C/RCX  (  O)  *RC*(  O)  *OC(  69 

8.  C9Mm*n  /c</Ri<X{60)»:?l<y(6o),0<f69) 

9.  09MM8N  /TO?/a(  10)  *  Tf  10)  *N39 

10*  DATA  Pl/3«1.4i  593/ 

11.  AO  I N  A^O  PRINT  9UT  H'P'JT  OAta 

12.  100  P8HMAK6I5) 

13.  101  <t9RMAT(2E2o*7> 

19.  102  =,99maT(6£11.(.) 

15.  REAO< 105* 100)  N32 

16.  REAQ< 105*101 ) < *( I ) *T( I )* I«i*VGQ) 

17.  9UTPUT * NGQ 

1*.  RE AO ( i q5  * 1 00 )  ITM, ITr, IbCjVX.MO 

19.  AUTPUTi 1TM< ItE, ISC*nx,NP 

20*  REAO(  105. 102}  o^lo,apm 

21.  9UTPUT*Pui')J^pm 

22.  08  3  I • 1 . NX 

23.  READ(105,102)  pmJI(I) 

24.  StjTPuT/  J*P*l  r  (I  > 

25.  3  C8NTINUE 

36.  REAO(  105*  1C2<  aM'J*EpS*9£Ta 

27.  9UTPIJT  »  Amu#  E«S*  SETA 

25.  CALL  COATACP 

29.  READ (105*lQ0)  N^R 

30.  9UTPUT*NFR 

31.  09  50  INP»l*.'FP 

32.  REAJ( 105. 102'  rYC 

33.  *RITF(  109*  30o)  F*c 

34.  AK«2.»pI*EvC*SQRT(AMM»£P9)*i,r  f, 

35.  9UTP'JT*a< 

36.  09  7  I»1*NC 

37.  R<X( I )»RCX(  I  )*A< 

38.  R<Y(  I  )«RCY{  I  )*A< 

39.  0<( I )»DC( I )*a< 

40.  7  riNTlNgE 

41.  IF( ITM.NF. 1 >  a?  T9  3q 

4?.  TM  OoYgivjFO  r  j rLO  MATRIX 

<*3*  <»1 

44.  09  1  INM^'C 

45*  L»(In»1)*Nc+1 

46.  09  1  I 1  *  I N 

47.  0UM<K)»5FTA*7MNE( 19* IN) 

43.  Z<L)«TMNF(  I  IN  )  ♦OJMC  <  ) 

49.  <  a  <  4 1 
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5C  • 

_»L*1 

5 1  • 

1 

C9NTI\UE 

52» 

oe  2 

53* 

JMl « IM-l 

54 . 

L"  I M*  IM1  /2  +  1 

55 « 

08  2 

56  • 

Z<  I^*NC*(  !*•'««  )  ).T*\E<  IM, 

57. 

l*l*i 

58* 

3 

08NTINUE 

59* 

CaLU  T~1XCE  .  V*,PMII  /NX, BETA  ) 

6  0* 

CALL  0EC9*iP  (  vC ,  I»S/Zi 

51  . 

C--- 

5?. 

L-l 

63  • 

<»1 

6*. 

08  4  IX«liVX 

65* 

^91  T£  ( 108/206  )  JWIIUX5 

66  • 

«9ITE( 1Q3J30, ) 

67. 

08  5  I ■ 1 #  NC 

6». 

XC<  I  ) «VM(K  ) 

69* 

<■<♦1 

70. 

5 

C8NTINUE 

71  • 

CALL  S8LvE(NCi !?S.Z/XC#C J3) 

72. 

09  6  I«1<MC 

73* 

Cm*CaBS  <  CUP ( l )  ) 

74. 

V"<L) «CUP< I ) 

75. 

*WITE<  103,20?)  I ,  C  J-?  {  I  )  j  Cy 

76. 

'.•L*l 

77. 

6 

C8NTINUE 

7?. 

4 

C8NTINUE 

79. 

I E  (  ISC  »ET .  l  )  Call  fS  (  V" ,  V'X  ,  p~  I!  /  ^  I  ,  03*  I  / ' 

30. 

300 

rsur'un  •  O’  ,  lsx,  '  T«  C'jR9EVTc  f**P  »uj 

31* 

?01 

PBK'lATt'  '/6v<«3jlsE  S9«  '  /  11X.  '  PEaL  '  /  1 1  <•  '  ! r 

33. 

202 

F9P^AT(*  ' , li5#6X/3El5.6) 

33. 

30 

!E(  ITE.MEM  )  39  T9  60 

3*. 

c--- 

.....F9PM  JPOEO  PT,  T2 1  A\3LE  T3  l  «»TPIX, 

35« 

<■1 

36. 

09  31 

37. 

L" ( IN»1 ) *NC*« 

3». 

09  31  IM.i,  i». 

39. 

OUM (<)*3ETA*(SZ^(I8-l#^-l#l*t)f5?H(l--*»l\i 

°c . 

1  *SZ-4  (  IM»  IN-1,-1/1  )*SZ*(I*«»TN#«1,-1)  ) 

91  • 

z  (L> »ou3(<  rqi4(  i  —  i/  in-1  ,  i .  i  i ♦« t-* <  i-*-i ,  r  / : 

92. 

93* 

<■<♦1 

94. 

L  ■  L  ♦  1 

95. 

31 

C8NTINUE 

96. 

09  32  1^-2/Na 

97. 

IilT'-l 

93* 

L»I-i*I^1/2a1 

99* 

08  32  I  N •  i ,  1  >•! 
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1 


1  30. 

Z<  I*+NC*  (!“•-.)  )  «DL'*CL  )+ST“*<  TM-1 .  I W  .  l«  1  )  *S  T-t  (  I -- 1 ,  Im,-:  ) 

-101* 

1  ♦STH  (  IM<  IN-1.-1M  )“STM(  I'*/  TN/-li-l  > 

10?. 

l*l+i 

103. 

3? 

-5NTINUE 

1  J4  . 

c  •  •  • 

.-..-F9RM  TE  E  yC  I  T  A  T  I  3\1  vECT«*9S. 

105. 

CALL  TEXCF( Ii\x#3ETA> 

106. 

CALL  OEC*Me(vC# I^SfZ  ) 

107. 

c + •  - 

.-.•-S5LVE  N9^vALIZE0  Tr  C’.i9Rr\iTS . 

1  08. 

<»1 

109* 

_  « 1 

no. 

59  3“  IX,1,Ny 

Ill* 

*«ITE( iq8*3o*>  °wlI(IX) 

11?. 

wKITE< 108*30!  ) 

113. 

59  35  I  M  *  s< C 

11“* 

XC  (  I  )  •  VVMK  ) 

Il5« 

<■<♦1 

116. 

15 

csntinue 

117. 

CaLL  59Lv£{\r, I»S«Z*XC/CjR> 

11?* 

59  36  I»l>SC 

its* 

C^-CAtJSlCUCU  r  )  ) 

i?o. 

V'ML)  !  ) 

1?1 . 

«WITE(  n*#2rJ?)  I»CJ«<  l  J'C'1' 

!??• 

l»l*i 

’  ?3  • 

16 

CONTINUE 

l  • 

3*. 

C9NT I \UE 

1  ?5  • 

I?(  I  SC  »E  3  •  l  )  CALL  TES<V'%\X,B‘dt,'5H!-.,5:sHt,A:=>,  51 

6  0 

C9NT I\UE 

127* 

50 

C9\T  I  VJE 

1??- 

900 

rSKMATt  '  T  ,  1<5X,  '  TE  CL'99tMT5  s^p  pwj  .i.rjj.m 

129. 

9C  1 

FfJWMAT  {  '  '  ,3x,  '  T9  I  ANGLE  V?  . '  .  1 2*  •  '  ?E •  .  1 1 X  • '  t  m  a~,  .  , .  .  *  ,  .  -a3 . 1  ) 

no* 

90? 

C99~IAT  <  '  1  '  ,  '  r«E3'jEMCv  95  INJCI5EV-T  =>la'mE  *avE  a  1  .  r  l  S  .  7  ,  px  *  «  «-<Z  M 

131  • 

ST9F 

1  ??  * 

END 

Sample  output  is  given  by: 


POM 

\QCi  »  4 
IT-1  .  i 
ITE  »  1 
ISC  »  1 
MX  ■  1 

N?  a  4 

•  ‘OOCOOO 
UPhj  ■  99.3000 
I  *  I 

PmI I ( I )  •  *033033 

X"iU  ■  l  •  »q f 

EOS  »  g.?5o6ooE»12 
SET*  ■  10.3003 


M9. 

p*»TNTS  SO£C 

I  E  Y  I  M 

3  C9MT 

•  24 

Jl_  X 

JL  Y 

PCX 

pry 

DC 

* • 1 33 1  - 

33 

•99t5E 

00 

•376iE 

33 

*  *  4955E' 

’31 

.qgq^r--, 

*■*  J833E 

*02960 

00 

•150fE. 

33 

•  *  1 45?E 

93 

•  1993c-') 

• • 6oS5E 

5c 

*79360 

03 

•301  iE 

33 

••2310E 

93 

.  3903c--  ~ 

•• 7936E 

3  n 

*6j«bo 

00 

•231 qfc 

30 

-•331  IE 

90 

. 9P73C-C 

-.92165 

n  f* 

*  3S93E 

00 

» 1 4b?t 

39 

• • 1596E 

39 

*  3993c  »o 

-•9915E 

r> 

■j  J 

•1331E 

00 

• 49b?fc. 

31 

” *  3763E 

33 

• 997SE-9 

■ • 991 5E 

33 

•  1 331  E 

30 

•49556. 

31 

••3763E 

90 

. TSqKC.' 

• • 9?3ftE 

00 

*  36932 

00 

« 1 45?e 

33 

•  *  35 76E 

9  " 

. laqir.: 

•• 7936E 

'J  . 

«6085E 

00 

•23l3t 

39 

••331  IE 

9  9 

. 19115-9 

* • 6o?5E 

->  /> 

•7996E 

00 

•33UE 

r\ 

-•2310E 

33 

.qoaac.n 

•* J«33E 

•9296E 

30 

•3506E 

**  n 

-o  . 

•  *  1 45pr 

9  3 

, baqpr- - 

-•13315 

33 

•  99<  5E 

00 

•  376?£ 

3  n 

••4955E* 

91 

• 9935E  *  3 

•1331- 

33 

*991 5E 

oc 

•37635 

*> 

*‘*9S9E*''i 

,  TQ-jqr.  -> 

*  3g  33E 

O'' 

•9296E 

00 

•  393ft 

*>  n 

•1452E 

. aoqir-^ 

.6365E 

.  79-J6E 

03 

•  301 1 1 

C  n 

•23’CE 

n>  ^ 

,  qa-ipr-" 

*  ^936E 

33 

«63«5E 

00 

•2913E 

•  331  IE 

^  *> 

.  qq-jpr-C 

.92365 

^  <*9 

V  _ 

•3S93E 

30 

*  1 4b?t 

00 

*  35C  6E 

39 

. iqipr-c 

*  99 1 5E 

n 
.>  . 

•13-15 

00 

•  <*9bKE- 

01 

*  3763E 

0  0 

(  qoqqr  •  ■> 

•  99 1  5E 

J  ^  • 

•13-15 

00 

•  49555- 

01 

• 37f 3E 

■V  •-/ 

•  aaosr. 3 

.92365 

*  36930 

30 

•  1*»52E 

0  n 

*  3536E 

33 

,  lonpr .  - 

• 7936E 

~  —  • 

•  6  0  fl5 o 

00 

•211 -5 

0  ~ 

•  3  3 1  1  E 

93 

,  7onr.' 

• 69455 

^  a  ^ 

.7396" 

30 

•  3- 1 1  E 

"9  ^ 

•  ?  j*  pc 

•>  - 

,  101-3C.' 

•  3?  33E 

,  qp-5*,r 

J  C 

•  JSCft 

00 

•  1  4f?E 

-  - 

•  1  q  1 3  r  -  " 

•  130  IE 

•79*30 

30 

•  376  35 

•>  - 

-  J 

, *3s;r . 

-9  1 

,  laisr." 

SCO  a  1 


FREO’JE'vCY  'jr 
a<  «  fc.28'99 

BIJLSE 


INCIDENT 

PLANE  ■‘•avE  » 

.3000000E  03 

M.l? 

tm  CURRENTS  239  B-II  ■ 

«.  «£A| 

.ocoor  00 

I  MAG 

MAG* 

1 

••124978E  01 

•  1 6372®E  01 

.209973=- 

01 

2 

••381894E  00 

.1733972  01 

.  1 99451 2 

01 

J 

••191649E  00 

.1798902  01 

.179715=- 

01 

•597857E  00 

•  1 4 1 1 962  01 

.  153330=- 

01 

5 

•106931E  01 

•6220162  00 

•  1237062 

01 

6 

•918747E  00 

-•211561E  00 

.942790=: 

00 

7 

•3283S5E  0C 

-.5911 91 E  00 

•  676257=: 

00 

a 

••1784P3E  00 

-  *  422236E  00 

,4584532 

00 

9 

••?951-»7E  00 

-.6967  }7E«0l 

.3-3239=' 

?  ^ 

10 
«  i 

-.1494312  OC 

•1129552  00 

. 1 *43732 

o.o 

• 147271E-01 

.8637572-01 

.8762212. 

Ol 

*  * 

12 

."65343E-01 

.  103M5E-31 

.872C*“2. 

Ol 

i  3 

.36584,22-01 

.  1O379PE-01 

.872061r- 

■01 

*  4 

. 1472a5E“01 

.8&3774E-01 

•  87623  72  < 

•01 

* 

15 

«  4 

.,14A**?oe  00 

. 112*542  00 

.1843772 

00 

-•295153E  00 

•  •  69669*2-  01 

.3032402 

00 

A  J 

•  7 

• • 1 7«4^7E  00 

..422^3"E  00 

•4584562 

0^ 

1  8 

• 3283^62  CO 

-.5911932  00 

.6762602 

09 

1  9 

•9187i*6E  00 

-.2119592  00 

.9427392 

20 

•1069912  Cl 

• 6?29l 1 E  00 

.1237062 

01 

?  1 

.5978382  CO 

.1411992  :i 

. 1 53331 r 

01 

22 

•»l9l6s6E  00 

.1789902  01 

.1797152 

01 

23 

-.«8l836E  03 

. 1 73*972  01 

.  1">945C2 

01 

2“ 

••  1249-»9E  Cl 

.1887292  01 

.3099742 

01 

r“  SC  *  ri-£R?o  EJelo  9  a  t  t  ?  oHt  j 


•  ^  0  0  7  9  7  > 


Oh  f 

.coooe  oo 

•9000£  02 
•  1 a00£  03 
•2700E  03 


S0ST(S/L) 
•H31E  01 
•1288E  oi 
•2621E  01 
•1088E  01 


T9 1 AXGlE  V?. 

1 

? 

3 

4 

5 

6 

7 

8 
9 

10 


TE  C'J99?NTS  PJ7  ai-i [  , 
R£.\L 

••154648?  3 . 
••145Q42?  oi 
••112366?  01 
••508975"  00 

•  31 96^7?  73 
•135 HIE  01 

•  1 2  a  6  ^  8  ?  3 1 


•OOOOE  "'0 
1**3 

•lOl^gflc 

•  lHnSB?  01 
*14444 if  31 
,  I7?93cr  oi 

•l57«0aE  01 

•  1  0  «  7 1 1 ?  01 


*901097? 

n 

-.  765P3t?  02 

• 22c9q5? 

00 

-•105425?  01 

• .278208? 

r»  rs 
»-/ 

•• 7 l 3«66?  00 

• • 38«4i 4? 

CO 

•• 700233?- 01 

••261 154? 

cc 

,45304.  r  -v 7 

4  ^  - 

-*26i 1*2? 

•  *  388479? 

•  *278204? 

•22C9«9E 
•901C39E 
*128679? 
•105141? 
*  31 964 7E 

• ♦ 508978? 
”• 1 1 2  346? 

•  •  1  ‘•5043E 


*693445?  03 
, 495944?  33 

••  7001  45E»"'l 
••713856?  00 
••105424?  oi 
••765197?  73 

•  1 03644?  37 

•108710E  01 
•167908?  01 
•172835?  01 
,14464(7?  01 

•  1 1 4q51 E  0 1 


VA3. 

•  1*5131=-  01 
•134512?  31 
.1*3158=-  j, 
,190174=-  01 
.170924=-  31 
•151237?  0! 

•  1 29o25r  0! 
. 11 8212?  31 

•  1  0771  7=-  01 
.766195?  00 
.3 9 4675=”  00 
.549029?  00 
.715298?  79 
.549032?  00 

•394669?  33 

•  766184=-  33 

•  1  07716=-  01 
.118212?  0! 
•129025?  oi 
.  1  51237="  oi 
.170924="  01 
•130173?  01 
•143153=-  01 
.184513?  01 


T-  <3CATTE9[1'0  PIEL5  o4TT?9\  rgo  pmti  a  9; 


►  ST?o*  0 


=h  1 

•COCO?  GO 
•9nooE  02 
•  1  H 0 2 E  03 
•2700E  03 


SORT  t  S/L 1 
•1026?  01 
•6618E  00 

*  1 475E  01 

•  66 1 3E  00 
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A. 3  EXACT  SERIES  PROGRAM 


...-PP9GPA"  ^zcr^ic  !*IOUCO  •* 

wwsa 

scATT^ea  helo  »*ttw  is  »ls-  c-’oteo. 

CS("PUiX  A  ( 50 '  *  3  f  SO )  *  C  ( 50 )  «  0  ( 5C ) 

31I4ENSI9N  c—l  t SOI. CPUS « SOI 

Ot“ENSnN  3J(130).3y<10Q)  sT«*,u»-'N 

‘O/dO»£C^O0»uHl-wP,-. JZ»-P>S  .i 

'3«m5n  *0 ( 33 1 »P ! 2  * p l*.3 17 
PEAL  , 

0*T»  j/ (  O' »  1  <) /»?  ! ''S*  1*1593/ 

rjpriATISEl*.?' 

r9RMAT(5!Sl 

rORfAT ( *51 1 •* I 

R£*0< 105/100) (V3(l)»r*l* S3) 

=12*2./°! 

..!il*87*‘^8.  **  3ESSEU  FUNCTIONS  Taken 
35*0(105.101'  "3.NC.MI 

•UTPUT.N9.NC.Nt 

READ  (  105.  1035  49.3w0.0P*<t 

5UTPUT.aP.P“O.0pu! 

a<«AR*2.*P! 

31*2«/(°1**<' 
ra.3QPT( 2*/p* > 

CALL  3ES(N9«!.a<.3J.9T! 

*N1*9J<  1  )*U*9',(  1 1 

*0*W**1 

30*9J<1>/WM1- 

4P**9Jt  2 ) •U»9T ( 2 ) 

C0« 1 ./WP 

00**9J(2)/UP 

UN*U 

•1N.«1 

00  1  t*l»N9 

4*3 J ( 1*1  I *U*9y( t*t ! 

A ( I l*uN/« 

3(1) .MlN*3u( t *1 ) /N 

HP.4M1  *1  *N/Al< 

o' :  i*«tN»(9Ju  >  •  i  *9u  ( 1*11  /*<>  /mP 

4M1*w 

jN^U.uN 

*jn**min 

-ONTtNUC 

*R I TE ( 103. iOO 1  *0.30 
«R I TE ( 139. 1 30 '  CO. 03 

•sm  13  f  »  1  /  N*  ^ 

«a:Tei!,2*»20M  :*a(I),s<I) 

*3 1  TE (  139.201)  I.C(t)-OlI) 

-•NTINUE 

.....-•CS*'PUrE  r“  CUPPENTS. 

»9tTE( 108.200' 

0Pm*2.*pI/nC 

om.«0Ph/?' 

39  3  t*l«NC 
j2**0 

39  *  u*l«N0 


SO-  jZ*jZ*2«4A<  U)*C9S<°m.j] 

si-  *  continue 

52.  3m«pm.O°m 

S3-  uZ-Pl^UZ 

s*.  jzi»CAas( jzi 

S5 «  «RITE(13S<20l)  I.UZ.JZN 

66.  3  CONTINUE 

S7  •  C-  — — — CO^UTE  Tr  CURRENTS. 

53.  -RITE! 10S»20?I 

=>9.  3m«0 • 

70.  *9  5  I»1»NC 

71.  .P» CO 

72.  00  S  j«'.jN9 

73.  uP« JP*2 • »C ( J) «C9S  < ®h»J ) 

74.  6  CONTINUE 

75.  3M.PU.0PN 

76.  JP«U»P1«UP 

77.  jPN«Ca8S< JP) 

78.  .RITgl lOSi 20l 1  I.U»jjP“ 

79.  *  CONTINUE 

30.  0..-.-.--C3M(,'JTE  T«  anC  t;  NeR**At_I  7r*  SCAT-»rReD  RJ-i.O  pa*t-raiS 

31.  .RITEt 133j 207) 

32.  «RI T£ 1 108*  20b ) 

33*  sw.Pmo»oi/i3;. 

3*.  0°»0Pu 1  »0 1 /I *0 . 

35.  00  7  I o . 1 j N I 

36.  3T-.90 

37.  3TE.OO 

39 •  00  3  J*liN9 

39.  C°«C3S( J»Pu) 

90.  3TN.ST-*2..3( Jl.CP 

91.  STE*STE*J..O(-)»CP 

92.  3  CONTINUE 

93.  SE«P2*C*9S(STE) 

94.  SN.rs.CAgsiS**! 

95.  nRJTEi 108.20*)  »wo,S“,SE 

96.  bmo.<’nc*OPm  i 

97.  9M.PU.QP 

98.  7  CONTINUE 

99*  200  "5R71AT(  - 1  > , 'tm  currents*..  •) 

100.  201  PORMA T {  •  '/ ItO,4£lS.6) 

131.  202  PORO»T( '0> , --E  CURRENTS... • I 

132*  233  PORf*AT!  '  1 '  ,  '  T-  ANO  T£  N9R-Ai_Izr3  SCATTrorQ  ?»£'.'  oA^rrSN^t 

133.  2C*  PORMATI'  1 » 3c 1 1 , 4 ) 

13*«  205  p5RmaT( 'O' .3si ' .«X. ' T.i ,9x, * TE> ) 

135.  STOP 

136.  ?N0 
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Sample  output  is  given  by: 


9UM 

N9  *  10 
MC  •  24 
M  «  36 
AR  >  *382800 

PMO  a  *000000 

0°HI  «  10.0000 

••7560842E-03  .1961226E  01  ,14862275-06  -  . 3*55 1  A 1  E-o3 

• • 1 853659c  ol  • • 367551 4E  00  *962l7Q8r  Q0  *lq97«33E  00 


1 

-.2S7551E  00 

•185366E  ol 

•  •  9621  7i E 

• .  1 a^78?E  09 

1 

• « i 76372E  01 

•  * 81 5625E  oo 

-  •  1 761 8oE  00 

» 39  ^972E  00 

2 

•  *  1 1 7822E  01 

•  1 1 5772E  ol 

•50S773E  o* 

»,4qq9'33E  -'O 

2 

- « 6858Q2E  00 

•.195346E  01 

•109727E  00 

•3l?549E  39 

3 

••116608E  01 

-.287091E  00 

-*57i91oE»01 

• 23?1 31 E  90 

3 

,  l 55797E  01 

-.489399E  00 

-*a3S13oE-01 

-*2»5914E  90 

4 

. 254945E-01 

•*626?l6E  00 

•  1 65156E-0? 

-.40605°E-9i 

4 

. 569038E-01 

*542o7oE  00 

•345724E«,93 

.  4RK096E-01 

5 

•222761E  00 

.813862E-03 

•  •  1  3348oE-04 

•  3653A7E-0? 

5 

-•129210E  00 

.512326E-03 

••15721 6E *04 

•  •  3Q650lE»  32 

6 

-  •  1 1 &936E-Q4 

.585782E-01 

.098494E-07 

-  •  1 3:>623E«  33 

6 

• ,*4430^E*Q5 

-•262491E-01 

.  43O988E-07 

•  ?9736?E- 93 

7 

-•123943E-01 

-.923732E-07 

- • 555 A54E "1 o 

. 74K3S sr. 05 

7 

•  45925 / E-02 

-.349651E-07 

-•579638E-1 o 

-•  7  6 1  3  4  o  E • 0  5 

M 

•45o55*E*09 

-.220967E-02 

•415758E-1 3 

••293901 r«06 

8 

• 145061E-09 

•  701962E-03 

•AP7046E-1 3 

•  2966” 1 E»  06 

9 

.341316E-03 

.146022E-11 

-*1 33030E-16 

, 497B5^r«oa 

a 

• « a5l 423E-Q* 

*4io734E-12 

-♦1 36369E-16 

•  •  4 7  j 7O4E-08 

19 

-.3313Q4E-14 

.465686E-Q4 

•536134E-20 

.,71  i  4  7 1  r  -  1  o 

10 

-.829195E-15 

*  •  11 58o7E»Q4 

•512673E»;5o 

•7i A01?E»1 o 

Tm  C'J^EMTS  .  ,  , 


1 

••126455E  01 

•167788E  01 

•2101Q4T  oi 

3 

- • S96445E  00 

•177994E  01 

•  1  '59?94E  O'. 

3 

-.207342E  00 

• 177882E  01 

•179087E  01 

4 

• 57942CE  00 

•  1 40687E  01 

•152152E  01 

5 

•  1 Q4805E  01 

•625184E  OC 

•122036E  01 

6 

•901943E  00 

" •  196249E  00 

•q23046E  00 

7 

• 025q69E  00 

* • 569P56E  00 

• 635Q7aE  00 

8 

•« 1 66396E  00 

"•407463E  oc 

•440129E  00 

9 

-.280348E  00 

•  •6S5472E-01 

•288607E  00 

lo 

- • ! 39286E  00 

•106374E  00 

•175260E  00 

11 

. 139867E-01 

•834404E-C1 

•®46044E*01 

l? 

• *2595 1E«Q 1 

.  122661E-C1 

•  b35009E-'m 

1? 

- *25952E»01 

.  122641E-01 

•"35007E-01 

14 

•1399Q4E-01 

•834385E-01 

•"46033E-01 

15 

-.139280E  00 

•10&376E  oo 

•175256E  Oo 

16 

-•28Q346E  00 

-.685377E-01 

•289602E  00 

17 

« • 1 664Q7E  00 

••407451E  00 

•4401 22E  ~o 

18 

•  325948E  00 

-.569257E  oo 

•65597oE  00 

19 

•Q0l928E  00 

•  •  1 9fe?7 1 E  00 

•a:s3036E  00 

20 

•104806E  01 

•625153E  00 

•122035E  01 

21 

•  97944 7E  00 

•140685E  01 

•152151E  01 

22 

-•207312E  00 

•  177882E  01 

•179086E  01 

23 

-.896417E  00 

•177994E  01 

•199299E  oi 

24 

• • 1 26454E  01 

•  1 67788E  01 

•210104E  01 

CU9=?E'JTS. 

•  • 

1 

« 1 54852E  01 

-.101437E  01 

•1  351 1 *E  01 

2 

•144972E  01 

" • 1 1 3488E  01 

•134109E  01 

3 

•111 826E  01 

* • 1 43071 E  01 

•  1 31 58*E  oi 

4 

•536173E  00 

* • 169820E  01 

•177203E  01 

5 

-.3Q5532E  00 

-•164074E  01 

•166894E  ql 

6 

-•101326E  01 

••106239E  01 

•146812E  01 

7 

- • 1238Q4£  01 

"•113352E  00 

•124322E  01 

g 

» • 869S56E  00 

•724165E  00 

•113134E  01 

9 

■ • 31 8888E  00 

•100826E  01 

•103179E  0! 

10 

•261561E  00 

•6S8452E  00 

•736465E  0^ 

11 

•36995QE  00 

•721736E-01 

• 376a?4E  OO 

12 

•249369E  00 

-.460950E  00 

♦524080E  OO 

13 

•167391E  00 

••6645Q1E  00 

•  6 3526qE  00 

14 

•249365E  00 

-.460963E  00 

•524090E  00 

15 

•369948E  00 

•721528E-01 

•  3769 1 sr  00 

1  6 

. 26 1 57 IE  00 

•688433E  00 

•  73645 1 E  00 

17 

- «  ?1 8865E  00 

•100«26E  01 

•  1  03 1  7<*r  01 

18 

-.S69833E  00 

•7?41«5E  00 

•1131*42  01 

19 

-•123804E  01 

••113319E  CO 

•  1 24 32 IE  01 

20 

-•101328E  01 

-.106236E  01 

•146*1  IE  01 

21 

• • 30556 IE  00 

-.164Q73E  01 

•16639ur  oi 

22 

•506146E  00 

" «  1 698 21 E  01 

*1772032  01 

23 

•111824E  01 

-.143074E  01 

•  1 3 1 58QE  01 

34 

•144971E  01 

-.113489E  01 

•  1 34 1 1 OE  01 
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T‘l  avo  tc  \i»RM*L!ZE3  SCATTr^cO  c J  EL3  aATTC9vS 


am 

TM 

TE 

ooooE 

00 

•  1 1  3  aE 

01 

•1C20E 

01 

130CE 

02 

. 1133E 

01 

•  1330E 

01 

2000E 

02 

•  1 1 29E 

01 

•  1057” 

01 

30CCE 

02 

•1122E 

Cl 

•1082P 

01 

4000E 

02 

•  1 !  12E 

01 

•1082p 

01 

5000E 

02 

. 1099E 

01 

•1033c 

01 

£>000E 

02 

.1037E 

01 

• °214p 

00 

7000E 

02 

• 10<33E 

01 

•  76392 

00 

8000E 

02 

•  1037E 

01 

•6332E 

00 

9000E 

02 

•  1092E 

01 

•6539c 

CO 

1300E 

03 

•  1078E 

01 

.8261E 

00 

1100E 

03 

•  1034E 

01 

•101*E 

01 

1200E 

03 

•  9878E 

00 

•  1 120E 

01 

1300E 

03 

•  10*?E 

01 

•11222 

Cl 

l^OOE 

03 

• 1315E 

01 

•1080E 

01 

150CE 

03 

. 1737E 

01 

•  1 1 1  OP 

01 

16C3E 

03 

•  21  74E 

01 

•1255c 

01 

17Q3E 

03 

•2496E 

01 

•1421E 

01 

1300E 

03 

•2614E 

01 

•  1 491 2 

01 

1900E 

03 

•  2496E 

01 

•  14212 

01 

2300E 

03 

•2174E 

01 

•  1255E 

01 

2100E 

03 

♦ 1737E 

01 

•11102 

01 

2203E 

03 

•  1315E 

01 

•1080E 

01 

2300E 

03 

•10*8E 

01 

•  U22E 

01 

2400E 

03 

•  9*7SE 

00 

•  1 120E 

01 

2500E 

03 

•  103^E 

01 

•  101^2 

01 

2600E 

03 

•  1078E 

01 

•  8261E 

00 

2700E 

03 

.1092E 

01 

•6539c 

00 

2300E 

03 

•  10S7E 

01 

•63322 

03 

2900E 

03 

•  1 033E 

01 

•7639E 

00 

3000E 

03 

•  1337E 

01 

•9214c 

00 

3100E 

03 

•13 99E 

01 

•  10332 

01 

3200E 

03 

•1112E 

01 

* 1082E 

01 

3300E 

03 

•  1 122E 

01 

• 1082E 

01 

3*00E 

03 

♦1129E 

01 

•1057c 

01 

3500E 

03 

•11 33E 

01 

•1030E 

01 

•STQP*  0 


A. 4  PROGRAM  TO  COMPUTE  CURRENT  DISTRIBUTION  DUE  TO 
IMPRESSED  LONGITUDINAL  CURRENT 

A  computer  program  is  presented  here  for  the  special  excitation 
considered  in  Subsection  2.8.  Subroutine  CDATA  is  used  to  specify  the  con¬ 
tour  C  as  was  done  in  the  E-field  program.  The  data  cards  needed  for  this  are 

exactlv  the  same.  The  total  axial  current  I  is  stored  in  the  FORTRAN  vari- 
'  z 

able  ZI.  DO  loop  1  forms  the  matrix  given  by  Equation  (105)  of  Section  2. 

The  system  of  linear  equations  is  solved  by  RDECOMP  and  RSOLVE. 

The  example  given  here  is  for  the  cross  section  of  a  fuselage  sta¬ 
tion  as  shown  in  Fig.  6-1.  A  total  of  28  subsections  were  used  to  approximate 
the  contour.  The  subsection  number  and  the  total  current  flowing  axially  on 
that  subsection  is  printed  out. 


I 


l,  C--------PP9GRAM  CSMOyyjsjQ  CLIENT  DIST"?!  TUT  1 8"  9N 

?.  C  PERFECTLY  CONDUCTING  CYL  IvDRIC4L  S^r|_L  C  ARRY  f  *<3 

3.  C  TOTAL  LONGITUDINAL  CURRENT  ZI. 

4.  COMMON  /CUV/ULX<60>*ULVC60)*NC/C/PCX(O>iPCYf 

5.  DIMENSION  CM(25D0) *V<{50) . ALPM(50W IPSI5C) 

6.  100  FORMAT (615) 

7.  101  F8R*14T<  2E20«7) 

«.  102  C9RMAT(6E11«4) 

9.  103  C9RMATI I10*2e15*7) 

10.  104  r8R‘lAT(2T5,El5.7»2I5»El5.7) 

11.  CE »  » 222158 

12.  READ( 105* 102)  21 

13.  9UTPUT #  Z  I 

14.  CALL  CDATA ( I  DC ) 

15.  D9  10  I»1*NC 

16.  RCX( I ) »RCX( I )*Cr 

17.  RCY( I )*RCY( I )*Cr 

18.  DC(I)»DCm**CE 

19.  10  CONTINUE 

20.  <»1 

21.  DP  1  J-1/NC 

22.  D2»DC(J)/2. 

23.  UX»02«UL*(J) 

24.  JY»D2*ULV<J> 
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25 

26 
27 
29 

29 

30 

31 

32 

33 
3* 

35 

36 

37 
98 
39 
“0 
»1 
*2 
*3 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


08  1  !«1»NC 

IE< I.EQ.J)  3*  T8  7 

RX-RCX(  I  )-RCx<  J) 

RY.RCV,(  !  )-RCv( J) 

RU- <  RX-UX ) #«5+ <  RV.yy ) **2 
RL,(°X+UX)*404.{  R  Y*UY  )  **2 
0TJ-ULX( J) *Ry*ULY( J) *RY 
DNJ-ABS('JLV(  j)*RX.JLX<  J)*RY) 
THl«ATAN(D2«f'TJ.DNJ) 

TM2-ATAN(  -02-0TJ,0NJ> 

C'1{<)«03*(  AL8G(RU*RL5/2*-2.  ) -0T  J* 4L9G  (  RU/RL  > /?. 
1*DNJ*< T41-TH?) 

<■<♦1 
39  T9  8 

7  C9(0-OC(J)*<AL8G(D25-1*  ) 

*  ■  <  ♦  i 

a  CONTINUE 

1  CONTINUE 

09  3  I-l/N C 

v<n  )*io» 

3  CONTINUE 

CAUL  RDEC9*P{NC# IpS#C*) 

CALL  RS9LVF(NC, IpS#CM,VKf ALpM) 

Sl-O* 

s«o* 

00  4  I«1|NC 
«  «S4ALPut(  I  j  *^C  (  t  ) 

si«si*Dcn ) 

4  Continue 
SUTPUTjS 
output /SI 
C-S/ZI 

09  5  I-l.NC 
ALpH( I ) »AL°Hf I ) /C 
wRITE( 108# 109)  I #  AL°m ( I  ) 

5  CONTINUE 


61 

62 

63 

64 

65 

66 
67 


6 


SC-O* 

09  6  I « 1 #  N c 

sc»sc4al°h( i )-oc( n 

C9NTINUE 

OUTPUT  #  SC 

ST9P 

END 


f 
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Sample  output  is  given  by: 


SUN 

l\  «  20000.0 

-NB.  9F  STRAIGHT  UI^E  SEGMENTS  AP°R9X I “AT  I  NO  C  - 


UL* 

UL* 

RCX 

flCY 

oc 

• 0000E 

00 

•OOOGE 

00 

• 10Q0E 

01 

•OQOOE 

00 

• 1Q00E 

01 

•3000E 

00 

• «500oE 

00 

. OOOOE 

00 

*  1000^ 

01 

.2000E 

01 

.5000E 

•01 

*  9988E 

CO 

♦4994E 

■01 

• .  1500E 

01 

• 25002" 

01 

*  100 1 E 

01 

»  3000E 

01 

•  1 OOOE 

00 

• 9988E 

00 

•49P4E-01 

*  *  250qE 

01 

•  75002* 

01 

•10012 

01 

•4000E 

Cl 

*  2000E 

00 

.9950E 

00 

•99502 

■01 

•  *  35002 

01 

•1500E 

00 

. 10052 

01 

.5C00E 

01 

•3500E 

oc 

• 9889E 

00 

.14*32 

00 

" *  450q2 

01 

.27502 

00 

*10112 

01 

•6000E 

01 

•  4500- 

00 

• 9950E 

00 

•99502 

■01 

• *55002 

01 

•40002 

00 

•  10052 

01 

•63Q0E 

01 

• 5000E 

00 

• 9864E 

00 

• 1644E 

00 

- «615qE 

01 

•4750E 

00 

•30412 

00 

•6000E 

01 

•600CE 

00 

.9487E 

00 

•31622 

00 

- »615oE 

01 

» 5500E 

00 

•31622 

00 

•5000E 

01 

. IOOOE 

01 

.9285 E 

CO 

•  371 4£ 

00 

• »55 OqE 

01 

.80002 

00 

•13772 

01 

•4000E 

01 

•  1 600E 

01 

• 8575E 

00 

•5145E 

00 

• «450o2 

01 

*  1 3qoE 

01 

•11662 

01 

•3000E 

01 

. 2700E 

01 

•6727E 

00 

. 73a9E 

00 

■  *  350oE 

01 

•  21 502 

01 

•  14*72 

01 

• 2000E 

01 

.  4000E 

01 

•6097E 

00 

•  7926E 

00 

• »  250qE 

01 

*  33502 

01 

•16402 

01 

•100CE 

01 

. 4600E 

01 

. 8575E 

00 

•5145E 

00 

" • 1 50qE 

01 

•  43002 

01 

•11662 

*1 

•0000E 

00 

• 475QE 

01 

. 9839E 

00 

•14832 

00 

•  *  5000E 

00 

• 4675E 

01 

•  ionr 

01 

* 1000E 

01 

•46002 

01 

.98S9E 

CO 

- . 1 4*32 

00 

•5000E 

00 

» 4675E 

01 

•  ionr 

H  1 

.2000E 

01 

•4Q00E 

01 

•8575E 

00 

••5145E 

00 

•  15002 

01 

• 430CE 

01 

• 1166T 

01 

•3000E 

01 

• 270CE 

31 

•&C37E 

00 

••79*62 

00 

•250cE 

01 

•33502 

01 

•  164C2 

01 

•4C30E 

01 

*  1 600E 

01 

• 6727E 

00 

• . 73q9E 

00 

•  35002 

01 

.21502 

01 

•  14 A7r 

01 

.50Q0E 

01 

. IOOOE 

01 

.8575E 

00 

• . 5 1 45E 

00 

*450C2 

01 

.13002 

01 

•11662 

01 

•6C30E 

01 

« 6000E 

00 

» 9285E 

00 

* • 371 4E 

00 

•  55002 

01 

•5000E 

00 

•  10772 

01 

.6300E 

01 

•5000E 

00 
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•  9437£ 

00 

•  • 3162E  00 

•  6150E 

01 

• 5500E 

00 

•  3162E 

00 

•6000E 

01 

» 450QE  00 

••9864E 

00 

-•1644E  00 

•  M5cE 

01 

•4750E 

00 

♦  304 1 E 

00 

.5000E 

01 

*  3500E  00 

• « 9950E 

00 

-•99<50E-qi 

•5500E 

01 

•4000E 

00 

•1005E 

01 

•4000E 

01 

.2000E  00 

••98S9E 

00 

-•1483E  00 

• 4500E 

01 

•2750E 

00 

•  131  lE 

01 

•3000E 

01 

•  1 OOOE  00 

•  *  9950E 

00 

-•995oe»oi 

•3500E 

01 

• 1 500E 

00 

•1005E 

01 

» 2000E 

01 

«5000E»01 

* » 9938E 

00 

• • 49a4E»01 

•2500E 

01 

•7500E« 

'01 

•  1001E 

01 

» 1000E 

01 

•OOOOE  00 

•  *  9938E 

00 

•  .4994E-01 

•  1 500E 

01 

•2500E- 

-01 

•  1 001 E 

01 

• OOOOE 

00 

•OOOOE  00 

•• iOOOE 

01 

•OOOOE  00 

•  5000E 

00 

• OOOOE 

00 

•  1000P 

01 

-C18SED  C8NT91R,  IOC  »  0 

S  ■  «69 • 7375 
SI  .  6.30543 


1 

•2200409E 

04 

2 

• 2247948E 

04 

3 

•2354470E 

04 

4 

•25Q4063E 

04 

5 

•2S57402E 

04 

6 

•4112C20E 

04 

7 

•  1 157018E 

05 

3 

•1Q33439E 

05 

9 

• 35n2833E 

04 

10 

•2170148E 

04 

11 

•2087022E 

04 

12 

•2673734E 

04 

13 

•  353751 7E 

04 

14 

•3763301E 

04 

15 

•3763308E 

04 

16 

•35P7514E 

04 

17 

•2678731E 

04 

13 

•2087030E 

04 

19 

•2170130E 

04 

20 

•3502849E 

04 

21 

•1033443E 

05 

22 

•  U57015E 

05 

23 

•  4U2012E 

04 

24 

•28^7416E 

04 

25 

•2594Q52E 

04 

26 

•33R4471E 

04 

27 

•2247945E 

04 

23 

•2200424E 

04 

sc  •  20000.0 

•ST0D*  0 
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APPENDIX  B 


PROGRAMS  FOR  THIN  SHELLS  OF  ARBITRARY  CROSS 
SECTION  AND  FINITE  CONDUCTIVITY 

Two  computer  programs  are  presented  here  for  the  determination  of 
the  field  interior  to  electrically  thin,  lossy  shells.  The  formulations  used 
are  the  impedance  sheet  approximation  derived  in  Section  3  and  the  exact  series 
solution  of  Section  4.  The  required  data  cards  are  described  and  a  sample  of 
the  program  output  is  given. 

B.l  IMPEDANCE  SHEET  APPROXIMATION  PROGRAM 

The  input  data  card  sequence  needed  for  this  program  is  almost 
exactly  the  same  as  that  used  in  the  E-field  program  of  Appendix  A.  In  fact. 
Table  A-l  may  be  used  here  where  the  card  containing  AMU  and  EPS  is  replaced 
by  the  card  sequence: 

ST  *  shell  thickness  in  meters 

SEPS  =  normalized  (to  free  space)  permittivity  of 
shell  material 

SIG  =  conductivity  of  shell  material  in  mhos/meter 

VCO  =  velocity  of  light  in  medium  outside  and  inside 
shell 

NMP  »  number  of  near-field  measurement  points  inside 
shell  at  which  shielding  effectiveness  calcula¬ 
tions  are  to  be  made 

XM(i)  =  X-Y  coordinates  of  interior  field  measurement 
YM(i)  points  in  meters,  i  =  1,2,...,  NMP 

Also,  FHZ  is  the  frequency  of  the  incident  plane  wave  in  hertz. 

The  main  program  is  much  the  same  as  in  the  E-field  case.  The  two 
polarizations  are  handled  separately  and  the  upper  right  triangle  of  the 
E-field  impedance  matrix  is  created  first.  To  this  an  impedance  load  is 
added  using  Equation  (9)  or  (11)  of  Section  3.  The  E-field  excitation  vector 


is  then  formed  and  the  resulting  system  of  equations  is  solved.  The  interior 
field  due  to  the  equivalent  polarization  current  alone  is  computed  bv  sub¬ 
routines  QNFMS  for  the  TM  case  or  PNFMS  for  the  TE  case.  These  computations 
are  defined  bv  Equations  (50)  and  (52)  of  Section  3,  respectively.  To  find 
the  total  interior  field,  then,  the  incident  field  must  be  added  to  these 
"scattered"  fields  at  each  interior  point.  The  shielding  effectiveness  cai- 
culation  may  then  be  done  according  to  Equation  (1)  of  Section  3. 
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«,  c--------t-i!n  Q«r_i  * 


p . 

a  I  N  pr*'joam  ro^  r.j-5  rO^1  /» 

3. 

r 

C^9«r'T'f*9  a  T"t\  S~f__  5!\3  T-“ 

4  « 

/.JV/  _x  !  60  )  i  „i_ y  :  6  0  '  » N  0/0 /pc <  •  -• :  »  5C v '  -  : '  ,  0'  ’  ^ 

5. 

/•.S/X-l  1  00  5  *  yvl  (  1  O*'  !  ,  (  1  O'  5  i  y  '<  '  1  ^0 

6  • 

/C</9<X(A1-)»9<V(fe^)10<;A!'' 

7. 

C9M“9\  /33 I / a ! 1 0 ! , t I  1 0 ! «* 00 

»  . 

0  !*ENSI  3‘  p«»  !  (  10  ) 

9. 

C9mPLE X  Z(lC~C).v*(5O'0>'Cj(-'')#tZ,(5'> 

<  0 . 

09MpU£X  .’«Z“*'EA3ZM#0S09r/?fT*A02XP/2?':#  -Z  0 

1 1  • 

0*Ta  p  I /^. i 4. 593// J/ ( C* ' 1 • 1 / 

12. 

1  00 

r^SMATt  615) 

13. 

101 

r?9*iAT  (  222C  * "» ) 

14. 

1  J2 

r49^AT(<,rij.„) 

15. 

103 

P4C?^aT  (  5215*6) 

16. 

PEAD< 105. ICC )  V 13 

17. 

3FA3( 105 A lcl ' ! * ( 1  )  • T ( n *  !  «1 >K is ) 

18. 

•UTP'jTaMO 

1 9  * 

PEAD< ic5a ICO )  l7** ItEa JSCaNX.NP 

2C. 

•JTP'JT,  !▼*,  I’Ea  I  SC. NX, Np 

21  • 

2EA3(  1Q5a  102^  o^IOaOp-i! 

22. 

5UTpUTa3hI0a''p‘4I 

23* 

09  4  I  •  1  .Nx 

2*  • 

9EAD(105.102)  04IIC) 

25. 

9JTPUT.  I  .pi~!  T  (  I  ) 

26. 

4 

CONTINUE 

27. 

9EAD( 105. 1C2'  ST,?EpS,S1Sa VCO 

28. 

•UTP'JT.STa  SEpS.  SI  3a  V00 

29* 

9EA0<105a100)  n^p 

30. 

9UTPUTaN*P 

31* 

09  6  I»1a\mP 

32. 

9EAD ( 105. 1C2 )  X  M  ( I ).^v( I 5 

33* 

6 

Continue 

3*. 

call  CDATA(  i -»c ) 

35. 

PEAO ( 105, 100 )  nfp 

36  • 

0°  50  I  NC » i .  'T9 

37. 

PE  AD ( 1 05  a  1 02 i  pwZ 

35* 

wRITEt 108a30p)  r«Z 

39. 

aF«2,*P!*PmZ 

<*c. 

AKQ«wF/'/C0 

STK»A<0*ST 

42. 

AMU0«4.*«I*1 .E-07 

43. 

fpso«1 »/ ( Amu3*vcc*vc0 ) 

44  . 

3E7A.ST<*(u#fS2pS-l  .  )  4SI5/(  -*F*EBSc)  ) 

45  * 

9JTPUT. a«ug>2pS0a  wf 

46  . 

09  5  I  ■ 1  a  %C 

47. 

P<X (  I  ) »PCX (  I  1  * A<c 

-8. 

P<Y ( I ) »PCY ( I j *A<0 
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-  9  . 

■'<:  I  )  *  C  '  :  :  ;*»<C 

"2  * 

cq'»t :  \v#* 

^  1  . 

^ A  7  !  *  1  #  N 

52. 

53. 

*M<  (  »  (  ■  (  *  )  # A  <  7 

^  *  • 

y 

C*’sT  !  NLP 

<■1 

S  6  • 

Z*  3  ! • l > \ X 

5P  • 

zstc^st3-: : <  * )*3i/i ac,  > 

pfl  . 

«m.s : s ( ->-i :  1 1 1  *o ; /:  s^. ) 

-9. 

c*  a  _, » *  i  ni-o 

*3  C  • 

tZ !  !<:«C5!xp:  •  , )  »Sv  :  ) 

S  1  . 

<  •  <  •  i 

5?. 

a 

C^\  T j \  jr 

S3. 

1"  (  !T***\p.l  1  3*  ’*  3? 

*  *  . 

' . .  - 

S  5  • 

<•1 

6  6* 

C*>  1  I  N  •  1  a  *  C 

A  7  . 

Z 9  l  !“»!<!  n 

S  ?  • 

z ( < )  »zms-7  ( i v,#  rv i 

<  ■<  ♦  i 

iZ  . 

i 

C9NTI\or 

71  . 

r  »  »  - 

7?  . 

<•1 

P3. 

0^  1C  1 • 1  *  S C 

7<*  . 

Z  <  < ) «Z  <<!*?<(  :  )  /3t’i 

75  • 

I ♦  1 

7 1  • 

1  - 

C?NT  I  'jjr 

77  . 

. f?i)M  ™  EyC!T4’,!?k  VtCTo^*:, 

7  ®  . 

C  *  UL  TMt ( V-|Sh[  I , \ X  ) 

79. 

...-.59Lvr  \*pva'_  I  zrC  t*-  3  ’  C 

sc  • 
si . 

a  p . 
8?. 
a  <* . 
S5  * 
a*. . 
*7  . 
a  a . 
so . 

?C  • 
Pl  . 

3  ?  • 
o  p . 

9  •  . 

3^  * 
J6  . 

p  ’  . 


-P  ! 


-9«sc* { VC*! )/2 

C*Lu  jfc.St  vM,Z#‘ C#'-*, 

'  »1 

,x»  1 

C?  2  !"1»nx 

.91 Tr ( i:«, ?c« ) 

Cl  3  j* i < sr 

V'a'-a^ABS  (  V'-1  (  <  )  ) 

C  J  (  -  )  •  V  '•  f  <  ! 

*9 I  Tr  :  1  z*,t  ?'J5 
<  ■  <  ♦  i 
ISS’JN^r 

C  •  9  !  ~  \  “  d 

,ZIZ> 

'zc*pzc*rz:(  < ! 

•j  TS/,-’; 

1  :  S  #  1  P  ’  <  v.  f  ;  -  ,  #  V  - 


J|VMi^);Vuw 


/- 


Sf 


3  ft 

4  o 


"1 


130  • 

_  X  *  L  X  4  1 

101 . 

q 

CONTINUE 

132. 

p 

C9NTINUE 

133. 

IP  *  ISC.EQ.l  >  CALL  T^5{ 

104. 

230 

R9RMAT  (  *  0  '  ,  15X# '  T*"  CURRENTS;  e-49  Puj  «  .  ,  r  .  4 ,  1  3c  ~9rrq  ' 

105. 

231 

=*0pv,  A  T  (  •  '  ,6X> '  R'JLSE  N9« '  1 1 1  X.  '»EaL  '  t  1  1  *•  '  P‘5*  .  1  1X»  '  mA3. 

136. 

202 

e-9RMAT(*  '  ,  1 1 5*  5x >  3E  1  5 . 6  ) 

107. 

203 

r^RMATt  '0*  t  '*M!EL0ING  EprECT T VEVESS  r9R  T-1  CAqri) 

138. 

30 

I F ( ITE.NE.t )  39  T9  60 

139. 

C  — 

-----F9RM  LDpER  RT.  T  =?  I  ANJLii_r  =*r  TE  1  MATRIX. 

no. 

<■1 

Ill* 

I3«l*I0C 

112. 

09  31  IM» IB/MC 

113. 

09  31  n-IB/T* 

114. 

Z(< ) »SZH ( IM*1# J^-1>1/ 1 )*SZM( T*«-1a IM.1,-1 ) 4 

115. 

1?ZH( IM,  IM-i# -1* l )  4SZW(  IM,  IV# -li  -1 ) 

116. 

<»<4l 

117. 

31 

CONTINUE 

113. 

-.-..ado  L9AD  impedance  matrix 

119. 

nu«nc-ioc 

120. 

1E< IOC.EO. 1 )  G?  79  34 

121* 

0L*O< ( NC ) 

122. 

OU“D< ( 1 ) 

123. 

G9  T9  35 

124. 

34 

OL»OK ( 1 > 

125. 

0U»0K(2) 

126. 

35 

<•1 

127. 

09  36  I ■ 1 #  NU 

128. 

Z  (  K  )  «Z  ( K  )  4  (  Dl.  +  O'J  )  /3 » /BE  T  A 

129. 

OL-OU 

130. 

OU»D<< I + I 3 ) 

131* 

<»K4l+l 

132* 

36 

C9NTINUE 

133. 

<«2 

134. 

NU1 »NU*1 

135. 

09  37  I  -  1 #  NU1 

136. 

Z(O«Z(K)40K(  I4I0C)/(6.*9ETA  ) 

137. 

<»K4l*2 

138. 

37 

C9NTINUE 

139. 

lr< I0C.E3. 1 )  G9  T9  38 

140* 

<«NC*(NC-1 )/P4l 

141* 

Z(K)»Z(<)40<(NC)/6./9ETA 

142. 

3« 

C0NTINUE 

l*'i. 

c--- 

- -E9RM  T£  EyCl7ATI9‘i  VEC79R5. 

144. 

CALL  TEX(V*,PHI I,\X, IOC) 

145* 

09  40  I ■ 1 / NU 

146. 

40 

CONTINUE 

147  . 

c*-- 

_.---59LVE  FPR  N90,iAl I ZEO  7E  C'jRRrNTS. 

148. 

MR»NU* ( NU4 1 ) /? 

149. 

CALL  GELS ( VM, Z/ VU# NX# VR ) 

1 50 . 

LX*  1 
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151  . 

<»1 

152. 

09  32  I «  1 '  X 

153. 

*RITEU08j30«)  pwiI(I) 

15<*. 

*R I TE ( 108*30! 5 

155. 

09  33  J«!a,\o 

156. 

Vm*i»ca0S  (  v«  (  * )  ) 

157. 

CJ(-jT»-vm(k  ) 

15  3. 

*RITE(  108*20?) 

159  * 

<■<★1 

1 6C  • 

33 

C9N7JNJUE 

161  • 

w  1  T  E  <  l05i3CT) 

152. 

09  39 

16  3. 

CALL  PMEvis(CJjM2C<xm<!  p)(y^(IM]  , 

16<*. 

*ZC**ZC*EZ  I  ( i  X  ) 

165. 

9UTPUTi^jC 

16  6. 

SE  *?0**4L93iO(CA9S{EZ!tLX)/uZC) > 

167. 

*RI  TE  ( 10a<  103)  X^(  l^j, yr  (  !«) ,  ?r 

16  8* 

L  X  »  L  X  + 1 

169- 

39 

C9NT I NUE 

170. 

32 

C9NTINUE 

171. 

Jp  (  ISC  «E0. 1  )  CALL  TES(VM#\XiP,-iII#omo.OP^!#'J!S.!0'‘) 

172. 

60 

C9NTI\UE 

173. 

*50 

C9NTJNUE 

17*. 

300 

P9R*AT(  «  ■><  ,  i  =  Xj  ClJRREvTR  r<»3  puJ  ,  .  ,  c  i  j  ,  4  ,  -JX  .  '  0C  ",5:::  '  ; 

175. 

301 

E-^RMAT  (  '  1  t  3y/  *  TR  I  A\3lE  vo  .  '  >  12X#  1  RPaL  '  .  1 1  x>  '  t^aV  ,  ll  x,  '  -"a  3. 

176. 

302 

C9RMAT  (  '  1  '  ,  '  FREQUENCY  «e-  pL4me:  *a  vE  »  *  iE15»7#  3X,  '  ^2  '  5 

177. 

303 

f9RviaT  (  '  0  ' ,  '  I  ELD  I NG  EPPE^T  T  VE‘ EES  r^R  Tr  CASE  1  ) 

175. 

ST9° 

179. 

END 

Sample  output  is  given  by: 


PlIM 

NQG  «  a 
IT*  »  1 
ITE  •  1 
ISC  •  c 
MX  ■  1 
IMP  ■  l 

PHI3  •  .acoooo 

DPMI  »  5*00000 
I  ■  1 

PWIKI)  .  .000000 
ST  .  9*999999E-C4 
SEPS  •  1.00000 
SIG  ■  10000*0 
VCO  •  3.000000E  03 
MMP  •  i 

-N8.  8F  STpAlGMT  SE3MEMTS  APOp?X  I  mATI  MG  C  »  *6 


JLX 

JIT 

pCx 

9CY 

oc 

•5000E 

00 

•OOOOE 

00 

•4620E 

00 

-.1910E 

00 

-» 1951E 

00 

-•9308E 

00 

•  431 o£ 

00 

-•955CE 

•01 

• 1947E 

00 

.3540E 

00 

•  « 354QE 

00 

• . 5523E 

00 

-•8336E 

00 

»408oE 

00 

-•2725E 

00 

•195SE 

00 

•1910E 

00 

-.  A620E 

00 

-.8336E 

00 

* . 55?3E 

00 

*  2725E 

00 

-•*0®0E 

00 

• 1955F 

00 

• OOOOE 

00 

• «500QE 

00 

*.9808E 

00 

*  *  1 951  £ 

00 

•9550E 

-01 

*  *  48 1 OE 

00 

* 1947E 

00 

• • 1 91 OE 

00 

••*62QE 

00 

-.9808E 

00 

•  1 951 E 

00 

-•955oE 

*01 

-•4S10E 

00 

•1947E 

00 

-.3540E 

00 

-•3540E 

00 

• . 8336E 

00 

•552 JE 

00 

-•2725E 

00 

-•4C3CE 

00 

•1955E 

00 

••*620E 

00 

•*  1910E 

00 

-.5523E 

oc 

•3336E 

00 

•  *  4980E 

00 

-.2725E 

00 

. 1955E 

/■*  A 
V 

-.5000E 

00 

•OOOOE 

00 

•« 1951E 

00 

•98C8E 

00 

••481oE 

00 

-  *  9550E 

•01 

•  1947E 

00 

-.'*620E 

00 

• 1 91 OE 

00 

*  195  IE 

00 

•9803E 

00 

- » 4? 1 OE 

00 

.  955CE 

•01 

*  1947E 

00 

*.35*0E 

00 

•3540E 

00 

*  5523E 

00 

.  3336E 

00 

■  *  408  QE 

00 

•2725E 

00 

.  1 955E 

00 

" « l 91 OE 

00 

•4620E 

00 

• 8336E 

00 

. 55P3E 

00 

-•2725E 

00 

•4080E 

00 

•1955E 

CO 

.OOOOE 

00 

. 5000E 

00 

•9808E 

00 

•  1 951E 

00 

• »  955qE 

-01 

« 43 1 OE 

00 

•  1 9  4  7  E 

->n 

•*  V 

. 1 91 OE 

00 

• 4620E 

00 

•9808E 

00 

•• 1951E 

00 

•955oE 

-01 

• 43i oE 

00 

*  19**7E 

00 

.35^0E 

00 

•354QE 

00 

. 8336E 

oc 

• • 55?3E 

00 

•2725E 

00 

• 4080E 

oc 

*  1955E 

00 

. 4620E 

00 

•  1910E 

00 

. 5523E 

00 

- • 3336E 

00 

•408CE 

00 

•2725E 

00 

•  1955E 

A  A 

. 500CE 

00 

OOOOE 

oc 

•  1 95 1 E 

00 

••9303E 

oc 

«49iCE 

00 

•9550E 

•01 

• 1947E 

,-3  A 

'J  j 

-CLOSED  C 

^MTSUPi  IOC  « 

0 

B-7 


FREGUENC*  qr  3LA'JE  aavE  »  »1 OOOOOOE  0®  uZ 

AMUO  •  1.256637E-06 
E°SO  •  8*84194BE«12 
*F  «  6.283186E  07 


TH 

CURRENTS  F9d  owj  ■ 

•ooocf  00 

■'EGRcitg 

3ULSE 

N9  • 

Rr  Al 

I«A3 

HAG. 

1 

•378555E  Cl 

-.2691 16E 

01 

,  4/S4440E 

01 

2 

•35C129E  Cl 

-.275624F 

01 

•  445599E 

01 

3 

♦295902E  01 

-.284401E 

01 

•  4 1 041 6E 

01 

4 

«?242q9E  oi 

• . 29qk25E 

01 

.367035E 

01 

5 

•1475f9E  oi 

• . 29 l 94 1 F 

01 

•  32676 1 E 

01 

6 

•7652S2E  00 

-.237334T 

01 

.  R97350E 

0 ' 

7 

«?212a,5E  CO 

• . 28^01 4E 

01 

.?»0®87E 

01 

8 

-•706900E-C1 

« . 274?5°E 

-1 

.27435CE 

01 

9 

■ . 7Q69b 1 E  *0 1 

».274?59E 

M 

•274350E 

01 

10 

•221276E  00 

*  *  ?8000qE 

01 

.  2®08S2E 

31 

1 1 

•7652S4E  00 

-.2873349 

01 

• ?q735CE 

01 

12 

•1475&8E  01 

-.291S41E 

01 

.3*6760E 

01 

13 

.2243ooE  01 

-.29oS25E 

>  * 

.  3  6  7  0  3  6  f 

01 

14 

•2959o2E  01 

-.284402E 

01 

.4’ Q418E 

01 

15 

•35Q1F8E  Q1 

••275626E 

01 

.445600E 

01 

16 

•3785?5E  oi 

•  *  2691 1 ® E 

01 

.464441E 

01 

ELDING  FrrrCTI vE^ESS  F9R  T*  CASE 

:  «  (4*00781 6E« 

04* -8.124560E-04  ) 

.OOOOOOE 

00 

.OOOOOOE  00  .608581E  0? 

TE 

CURRENTS  fdo  P»tl  « 

.OOOOF  CO 

OEORFr^ 

TRIANGLE 

NO. 

REAL 

I M  A  3 

MA3. 

1 

.981393E  00 

•2?493?E 

00 

.  1006S4F 

01 

2 

» 982880E  00 

.20900QE 

00 

. 1004869 

01 

3 

•9864S3E  00 

• 16315®E 

00 

«9q9865E 

00 

4 

.  989852E  00 

•943096E- 

01 

.994335E 

0^ 

5 

*99o7a7E  no 

. 134543E- 

01 

.9qC8S9E 

00 

6 

•988479E  QO 

-.673963E- 

01 

.99Q774F 

00 

7 

•  98391 7E  CO 

• • 1 36?43E 

00 

» 9933Q4E 

oo 

8 

•979551E  00 

-. 182094E 

00 

.9P6333E 

00 

9 

. 9777R1E  QO 

-.198Q21E 

00 

.qq764lE 

0: 

10 

•979552E  00 

-.  1820<3i5E 

00 

.996334F 

00 

1 1 

•R83916E  00 

-.136243E 

00 

.9q3304E 

3^ 

2 

. 988479E  00 

• . 673°63E" 

01 

• 9q0774r 

3*1 

,  3 

»q9o799E  CO 

. 134S55E- 

ol 

•9q0890E 

o*i 

14 

•9898S4E  00 

.9431 OqE • 

0 1 

.Q043369 

00 

15 

•986462E  00 

. 163159E 

00 

.  9q9864E 

0“ 

16 

•982880E  00 

•209003E 

00 

.  1 0C486E 

01 

SHlrLOTNG  EFFECTIVENESS  F9R  Tr  CASE 

HZC  ■  (-4,7lU51E-04#-5.40l96?E-33i 

•  OOOOOOE  CO  .OOOOOOE  QO  .453160E  0? 
*ST?o*  0 
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B. 2  EXACT-SERIES  PROGRAM 


The  computer  program  listed  here  computes  the  longitudinal  component 
of  an  interior  field  at  the  center  of  a  lossy  shell  of  circular  cross  section. 
The  formulas  of  Section  4  are  used.  The  input  data  is  defined  bv: 

AO  =  outer  radius  of  shell  in  meters 

A1  =  inner  radius  of  shell  in  meters 

EPSB  =  normalized  (to  free  space)  permittivity  of 
shell  material 

SIG  =  conductivity  of  shell  material  in  mhos/meter 

FO  =  initial  frequency  in  hertz 

DF  =  frequency-run  increment  (hertz) 

NF  =  number  of  frequencies  at  which  computation  is 
desired 

The  required  Bessel  functions  are  computed  by  the  following 
subroutines: 


J  (z).  Y  (z) 

CBESO 

O  0 

J1(z),  Y1(z) 

CBES1 

H(2)(x) 

HANK02 

o 

h{2) (x) 

HANK12 

In  the  above,  z  is  a  complex  number  and  x  is  a  real  number.  The  variables 
Til,  T12,  T21,  and  T22  are  computed  according  to  Equation  (7)  if  | k_j_>aQ  |  <  RM 
Otherwise,  Equation  (21)  is  used.  Equations  (6)  and  (9)  are  thus  evaluated 
for  n  =  0  and  the  shielding  effectiveness  according  to  Equation  (1)  in  Sec¬ 
tion  3  may  be  determined. 


r«----«--PR0GRAM  To  c^P'JTP  SHIELDING  EFFECT !  \ZE"iESS  AT 

r  CENTER  0F  L9SST  CYLINDRICAL  pwrLL*  B0TW 

r  p9LARtZATr9VS  AT  \>*RHal  !VCJQENCE« 

C9^pLEX  BJ00,BJ0lABJll*»J13A«V03A®vni,3Vir'#avii 

Complex  han<*2*'-,*n,i<i2>E tb#<3#ho»hi>st< 

C9mPLEX  TnjTi?#T2i/T?2/ETq#'»#CSQoT 
COMPLEX  CC0S#CSlNiCl*C2'C3#">D,  A.3.DC#  TE»  Tw 
PEAL  <0 

DATA  U/(0*a1.  )/aP!/3»1A1593/*C/3»E  q/377  •  / 

DATA  £PS0/8*9M 94E-l?//P^/3» / 
no  F0RMAT  {  6E15*7  ) 

102  F0RHATP  O'  !6X>,!rREQ«,i8x''AB?;(<9Ai)',«Xj'K0Al'#I?X,,TE5;E'.llX 
**'TmSEM 

103  F0RMA T (  *  •  ,  1 «/ 4E1 5 • 7 ) 

C  A l  ■  INMEP  RADIUS  (IN  ^ETEPS) 

C  AO  ■  9UTE®  RADIUS  (IN  ^FTPPS) 

c  kq  •  WAVE vUMBER  0F  SHELL  ®E3I9m 

C  A<  ■  WAvEA'UHrtEP  eF  E«EE  SPACE 

0  ETB  •  RELATIVE  IMPeOAnCp  0F  SmE'-L  MATERIAL 

R£AD< 105/ 100?  AO^ Al>EPSB^S!G 
9uTPuTiA0i A1,EPSB>SIG 
DaAQ-Al 

READ  1 105* lOO )  FOjDF 
R£AD( 105/ 1C1 1  NE 
9UTPUT<P0aCP»mP 
SIGN-SIG/EPSo 

c--------pepp9rm  frequency  run 

WR I TE ( 108  a 1  OP  J 
00  1  I X  »  1  /  N‘F 
fhZ«FO 

wpa2,#PJ*FWZ 

<  0  *  W  P  /  c 

KB«CS3PT(EPSq-U*S!G''l/WF)*K0 

ET3»<0/<B 

ST<«K8*D 

AKO»KO*AO 

A<1 »<0*A 1 

C l *p I *kb*a 1 /? • 

C2*2  • *U/ ( P 1  *  A<0 ) 

C3*SQRT( A1/A0) 

ho»hanko?<  A<0) 

MlaHAN<12( A<0) 

HOP*PEAL(HANk'02(  A<1  )  ) 
u1R*REAL(h*n^12( a<1 ) ) 
oabsbCabs i <y«Ai > 

!F(DABS.GT,R“)  39  T9  2 

CALL  CBES0(<B*a0*9J0Cj9Y00) 

CALL  CBESO(<B*Al#BJOl/BTOn 
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50 ♦  CALL  CSES1 (<a**Oi8JlO^BViO) 

51.  CALL  CBESl(<R*Ai,BJU,3Y11 , 

52.  T11»C1*< -3J01*BY10*3J13*BY01 ) 

53.  Tl2.cl*<BJ10#9Yll«BJll#BV10) 

54.  T21»C1*(BJ01*3V00*3J00*B'<01  ) 

55.  T22»C1*<  -BJC0*B'*'ll+3Jll*BY00) 

56.  38  T8  3 

57.  2  CONTINUE 

58.  Tll»C3*CC8S(?TiO 

59.  T12»*C3#CSIN(ST<) 

60*  T21--T12 

61*  T22«T11 

62.  3  CONTINUE 

63.  C--------C8MPUTE  DETERMINANT  FOR  TM  CASE 

64.  DO«HO*( -wiR*Tll+H0R*T12/ET3  )+Wl* ( .Hi R*T21 *ET3*M0R*T2? j 

65.  D0»C2/D0 

66.  OUTPUT*  DO 

67.  SE»-20**ALOG10(CABS(00) ) 

68.  c--------ceMPUTE  determinant  for  te  case 

69.  DD«  B0*(-wiR#Tll4ETS*M0R*Ti2)+m*{.mR*T2l/ET3*-w3R*TR2) 

70.  D0-C2/D0 

71.  OUTPUT  *  DO 

72.  SH«-20.*ALOGtO(CABS(DO> ) 

73*  WRITE<108*100>  rHZ* DABS* AKj *  SB, SE 

74.  EHZ»FHZ+DF 

75.  1  CONTINUE 

76.  STOP 

77.  END 
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Sample  output  is  given  by: 


PUN 

AO  •  *500000 

a 1  a  **99000 

EPSB  «  1.00000 
SIG  ■  10000*0 
FO  •  1 « OOOOOOE  07 
OF  •  1 .OOOOOOE  07 

NF  a  1 

TES£  T"SE 


•A6Q4617E  32  • 6 j98?3%2 


FSEO.  A8S { <9 A i )  <0A1 

CO  •  {3*921350£-0***»a*02064ce-CA} 

DC  a  <•5*^6381 *E*0*i"4*950687f-335 

•ICOOOCCE  08  .4433989E  ^3  *10*5l03£  on 

♦ST0O*  0 
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APPENDIX  C 


CAD  HOMOGENEOUS  SHELL  COUPLING  ANALYSIS  PROGRAM 

An  Interactive  computer  program  was  written  to  determine  the  coupling 
of  an  exterior  electromagnetic  disturbance  to  a  circuit  situated  inside  a 
closed  homogeneous  thick  shell.  The  low-frequency  formulas  of  Section  5  are 
used  for  the  interior  penetration  fields.  The  purpose  of  this  appendix  is  to 
briefly  describe  the  CAD  program  logic.  This  consists  of  descriptions  of  the 
following: 

0  User-specified  input  data  for  a  given  problem 

°  Computations  performed 

°  CAD  results 

A  typical  problem  is  illustrated  in  Fig.  C-l  where  the  interior  cir¬ 
cuit  consists  of  a  two-wire  transmission  line  of  length  L  and  spacing  h.  The 
line  is  loaded  at  one  end  with  an  impedance  Z^.  The  excitation  shown  could 
result  from  a  NEMP  or  NSL  threat.  The  quantities  of  interest  for  computation 
are  the  open  circuit  voltage  and  short  circuit  current  at  the  terminals  shown 
in  Fig.  C-l.  From  this,  the  maximum  power  coupled  to  the  circuit  may  be  deter¬ 
mined  and  thus  various  parametric  design  curves  may  be  plotted  to  facilitate 
a  CAD  solution  based  on  given  burnout  data  for  the  load. 

A  block  diagram  of  the  computer  program  logic  is  shown  in  Fig.  C-2. 
The  various  subroutines  required  are  listed  with  capital  letters  inside 
blocks  which  identify  the  various  program  segments.  Subroutine  THREAT  allows 
the  user  to  define  the  threat  as  one  of  the  following: 

HEMP  -  Nuclear  electromagnetic  pulse,  plane  wave  excita¬ 
tion  given  by  Equation  (24)  of  Section  5. 

Input  H  . 

o 

NSL  -  Near-strike  lightning,  low-frequency  line  source 
excitation  whose  magnetic  field  is  given  by 
Equation  (2)  of  Section  5.  Input  R. 

DSL  -  Direct-strike  lightning,  impressed  longitudinal 
current  density. 
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given  by  Equation  (25)  of  Section  5.  Once  the  threat  is  specified,  subroutines 
THKMAT  and  SHPMAT  are  used  to  define  the  shell  enclosure  part  of  the  problem. 
Data  required  here  is 

a  =  shell  wall  conductivity 

d  =  shell  wall  thickness  (several  units 

acceptable) 

VSR  =  volume  to  emface  ratio  of  enclosure  (meters) 
Conductivities  of  various  materials  are  tabulated  in  a  data  file. 


The  interior  fields  may  now  be  computed  for  the  empty  homogeneous 
shell  according  to  the  formulas  of  Section  5.  This  is  done  in  subroutines 
DSLTG  and  NSLTG  as  follows: 

NEMP  excitation  -  Use  NSLTG  where  internal  field  is  given 
by  Equation  (18)  with  a  and  8  defined 
following  Equation  (24) 

NSL  excitation  -  Use  NSLTG  where  internal  field  is  given 
by  Equation  (18)  with  a  and  3  defined 
following  Equation  (1) 

DSL  excitation  -  Use  DSLTG  where  internal  field  is  given 

by  Equation  (36)  with  a  and  £  defined 
following  Equation  (1). 

To  complete  the  specification  of  the  sample  problem  illustrated  in 
Fig.  C-l,  the  user  must  input  data  which  defines  the  circuit  under  considera¬ 
tion.  This  is  done  in  subroutines  HEIGHT  and  TRANS.  Data  required  here  is: 

h  *  spacing  of  transmission  line  or  effective 

area  of  standard  cable  (tabulated  according 
to  cable  RG  number) 

L  =  length  of  line  in  meters 

Z  «  characteristic  impedance  of  line 

o 

Z^  =  load  impedance 


C-4 


Finally,  the  open  circuit  voltage  and  short  circuit  current  at  the 
transmission  line  terminals  shown  in  Fig.  C-l  may  be  computed.  This  is 
accomplished  by  the  following  subroutines: 

VOC  -  Computes  VQc  according  to  Equation  (15)  of 
Reference  [4],  p  6-6,  in  Section  5 

SCCUR  -  Computes  I  according  to  Equation  (16)  of 
Reference  [4],  p  6-6,  in  Section  5. 

Once  these  computations  have  been  made,  various  optional  CAD  curves  for  a 
given  problem  may  be  plotted.  A  sequence  of  typical  plots  for  a  problem  is 
given  by  Figs.  6-5  through  6-11  of  Section  6. 
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